质因数分解 2012年NOIP全国联赛普及组

 1 #include<cstdio>
 2 int n,a,s[6000];
 3 bool v[600000];
 4 bool ps(int x){
 5     if(x<=50000&&v[x]==1) return 0;
 6     for(int i=2;i*i<=x;i++)
 7     if(x%i==0) return 0;
 8     return 1;
 9 }
10 int main(){
11     for(int i=2;i<=50000;i++){
12         if(!v[i]){
13             a=i;s[++s[0]]=a;
14             while(a<=50000){a+=i;v[a]=1;}
15         }
16     }
17     scanf("%d",&n);
18     for(int i=1;i<=s[0];i++){
19         if(n%s[i]==0&&ps(n/s[i])){
20             printf("%d\n",n/s[i]);
21             return 0;
22         }
23     }
24 }