hdu 1460 完数
注意：num1和num2的大小未知，需比较！
有两种方法：
法一：素数打印+素数分解（求因数和公式）
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #include<string> 5 #include<cmath> 6 #include<algorithm> 7 using namespace std; 8 int p[10005];//记录包括自身的完数个数 9 bool vis[10005]; 10 int prime[1229+1];//如何估算比n小的素数的个数？？ 提前打印即可，小于10000共1229个素数 11 void print_prime(){ 12 memset(vis,false,sizeof(vis)); 13 int i=2; 14 int index=0; 15 for(;i<10000;i++){ 16 if(!vis[i]){ 17 prime[index++]=i; 18 } 19 for(int j=0;j<index&&prime[j]*i<10000;j++){ 20 vis[prime[j]*i]=true; 21 if(i%prime[j]==0){ 22 break; 23 } 24 } 25 } 26 } 27 int work(int n){ 28 int i=0,ret=1,total=1,temp=n; 29 for(;prime[i]*prime[i]<=n;i++){ 30 int sum=1; 31 int num=1; 32 while(n%prime[i]==0){ 33 num*=prime[i]; 34 n/=prime[i]; 35 sum+=num; 36 } 37 total*=sum; 38 } 39 if(n>1){ 40 total*=n+1; 41 } 42 return total-temp; 43 } 44 void get_p(){ 45 memset(p,0,sizeof(p)); 46 p[2]=0; 47 p[3]=0; 48 int i=4; 49 for(;i<10000;i++){ 50 if(i==work(i)){ 51 p[i]=p[i-1]+1; 52 } 53 else{ 54 p[i]=p[i-1]; 55 } 56 } 57 } 58 int main()//10000 59 { 60 int n; 61 cin>>n; 62 print_prime(); 63 get_p(); 64 //cout<<p[6]<<endl; 65 while(n--){ 66 int num1,num2; 67 cin>>num1>>num2; 68 if(num1>num2){ 69 num1=num1+num2; 70 num2=num1-num2; 71 num1=num1-num2; 72 } 73 cout<<p[num2]-p[num1-1]<<endl; 74 } 75 return 0; 76 }
法二：筛法
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #include<string> 5 #include<cmath> 6 #include<algorithm> 7 using namespace std; 8 int main()//10000 9 { 10 int n; 11 cin>>n; 12 while(n--){ 13 int num1,num2; 14 cin>>num1>>num2; 15 if(num1>num2){ 16 num1=num1+num2; 17 num2=num1-num2; 18 num1=num1-num2; 19 } 20 int i=num1; 21 int num=0; 22 for(;i<=num2;i++){ 23 int j=2; 24 int sum=1; 25 for(;j<=i/2;j++){ 26 if(i%j==0){ 27 sum+=j; 28 } 29 } 30 if(sum==i){ 31 num++; 32 } 33 } 34 cout<<num<<endl; 35 } 36 return 0; 37 }
