1. 程式人生 > >UVa 10791 唯一分解定理

UVa 10791 唯一分解定理

可以比較兩種程式碼:

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<vector>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fod(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
vector<int
>
prime; vector<int>e; typedef long long ll; const int N=sqrt((2<<31)+0.5); int n; void Euler_prime() { bool vis[N]; memset(vis,0,sizeof(vis)); for(int i=2;i<=N;i++) { if(!vis[i]) e.push_back(i); for(int j=0;j<int(prime.size());j++) { if(prime[j]*i>N) break
; vis[prime[j]*i]=1; if(i%prime[j]==0) break; } } } ll solve(int n) { if(n==1) return 2; for(int i=0;i<int(prime.size());i++) { while(n%prime[i]==0) { e[i]++; n/=prime[i]; } if(n==1) break; } ll ans=0,num=0
; for(int i=0;i<int(prime.size());++i) { if(e[i]) num++; ans+=pow(prime[i],e[i]); } if(num==1) ans++; return ans; } int main() { Euler_prime(); int tot=0; while(scanf("%d",&n)&&n) { printf("Case %d: %lld\n",++tot,solve(n)); } return 0; }

第二種:

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<vector>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fod(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
typedef long long ll;
const int N=46340;
ll tot,f[N],Case=0,ans;
void solve(ll n) {
    ll m=(ll)sqrt(n+0.5);
    tot=0;
    for(ll i=2;i<=m&&n>1;i++) {
        if(n%i==0) {
            ll fac=1;
            while(n%i==0&&n>1) {
                fac*=i;
                n/=i;
            }
            f[++tot]=fac; 
        }
    }
    if(n>1) f[++tot]=n;
} 
int main() {
    ll n;
    while(scanf("%lld",&n)&&n) {
         ans=0;
         solve(n);
         if(tot<=1) ans=n+1;
         else for(ll i=1;i<=tot;i++) ans+=f[i];
         printf("Case %lld: %lld\n",++Case,ans);
    } 
}

第二種是AC程式碼

對於此題,我們肯定需要求出其唯一分解後的式子,但是問題在於n是在太大,n=2^31-1;由素數定理可知,小於n的素數最多可以由1e9多個,顯然空間是不夠的.但是我們可以在分解質因數時發現這樣一個事實:如果一直到sqrt(n)之後,其所剩下的質因數必然只會剩下一個或沒有,這樣我們便只需要到sqrt(n)就可以了,這個範圍是46340,完全可以承受