【UVA 10791】最小公倍数的最小和

<span style="font-size:14px;">特殊情况多，看注释</span>#include<cstdio>#include<cstring>#include<iostream>#define maxn 100000#include<cmath>#define LL long longusing namespace std;int cas;LL n;LL mul(LL a,LL b){if(b==0)return 0;LL ans=1;for(;b;b>>=1){if(b&1)ans=ans*a;a=a*a;}return ans;}int  main(){while(scanf("%d",&n)&&n){cas++;if(n==1){//显然printf("Case %d: 2\n",cas);continue;}LL ans=0,tot=0;bool ok=false;for(int i=2;i<=sqrt(n+0.5);i++){int cnt=0;while(n%i==0){cnt++;n/=i;ok=true;//如果本身是素数 }if(cnt)tot++;ans+=mul((LL)i,(LL)cnt);if(n==1)break;}if(!ok)ans=n+1;//本身是素数 else if(n!=1)ans+=n;//唯一分解最后还是剩一个素数 else if(tot==1)ans++;//只有一个质因数（他们的最小公倍数就是自身，所以就是他本身和一） printf("Case %d: %lld\n",cas,ans);}return 0;}