Minimum Sum LCM - UVa 10791

本文探讨了如何将给定的正整数分解为至少两个数的最小公倍数,并寻找元素和最小化的解决方案。通过素因数分解,我们可以找到满足条件的数集,进而计算出其元素和的最小值。

Minimum Sum LCM

\epsfbox{p10791.eps}

LCM (Least Common Multiple) of a set of integers is defined as the minimum number, which is a multiple of all integers of that set. It is interesting to note that any positive integer can be expressed as the LCM of a set of positive integers. For example 12 can be expressed as the LCM of 112 or1212 or 34 or 46 or 1234 etc.

In this problem, you will be given a positive integer N. You have to find out a set of at least two positive integers whose LCM is N. As infinite such sequences are possible, you have to pick the sequence whose summation of elements is minimum. We will be quite happy if you just print the summation of the elements of this set. So, for N = 12, you should print 4+3 = 7 asLCM of 4 and 3 is 12 and 7 is the minimum possible summation.

Input 

The input file contains at most 100 test cases. Each test case consists of a positive integer N ( 1$ \le$N$ \le$231 - 1).

Input is terminated by a case where N = 0. This case should not be processed. There can be at most 100 test cases.

Output 

Output of each test case should consist of a line starting with `Case #' where # is the test case number. It should be followed by the summation as specified in the problem statement. Look at the output for sample input for details.

Sample Input 

 
12
10
5
0

Sample Output 

 
Case 1: 7
Case 2: 7
Case 3: 6

题意:将一个数分解成至少两个数的最小公倍数,求着几个数的和的最小值。

思路:分解素数后,将素数的幂,就是这些数中的一个。

AC代码如下:

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
int vis[100010],prime[100010],num;
int main()
{
    int n,n2,t=0,i,j,k,p;
    ll ans;
    for(i=2;i<=100000;i++)
       if(vis[i]==0)
       {
           prime[++num]=i;
           for(j=i*2;j<=100000;j+=i)
              vis[j]=1;
       }
    while(~scanf("%d",&n) && n>0)
    {
        k=0;ans=0;n2=n;
        for(i=1;i<=9592;i++)
           if(n%prime[i]==0)
           {
               k++;p=1;
               while(n%prime[i]==0)
               {
                   n/=prime[i];
                   p*=prime[i];
               }
               ans+=p;
           }
        if(n!=1)
        {
            ans+=n;
            k++;
        }
        ans+=max(0,2-k);
        printf("Case %d: %lld\n",++t,ans);
    }

}



评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值