欧拉函数 模板

1. #include <iostream>  
2. using namespace std;  
3. #define N 100000001  
4.   
5. int euler(int n)  
6. {  
7.     int ans=1;  
8.     int i;  
9.     for(i=2;i*i<=n;i++)  
10.     {  
11.         if(n%i==0)  
12.         {  
13.             n/=i;  
14.             ans*=i-1;  
15.             while(n%i==0)  
16.             {  
17.                 n/=i;  
18.                 ans*=i;  
19.             }  
20.         }  
21.     }  
22.     if(n>1)  
23.         ans*=n-1;  
24.     return ans;  
25. }  
26. int main()  
27. {  
28.     int n;  
29.     while(scanf("%d",&n),n)  
30.     {  
31.         printf("%d/n",n-1-euler(n));  
32.     }  
33.     return 0;  
34. }  

 

下面是另外一种（筛选法）

[c-sharp]  view plain copy

1. #include <iostream>  
2. using namespace std;  
3. #define N 100005  
4. #define ll __int64  
5.   
6. ll eule[N];  
7. int p[N][15];  
8. int num[N];  
9.   
10. void init()  
11. {     
12.     int i,j;  
13.     eule[1]=1;    
14.     for(i=2;i<N;i++)//筛选法得到数的素因子及每个数的欧拉函数值         
15.     {  
16.         if(!eule[i])              
17.         {             
18.             for(j=i;j<N;j+=i)                  
19.             {                 
20.                 if(!eule[j])                      
21.                     eule[j]=j;                
22.                 eule[j]=eule[j]*(i-1)/i;                  
23.                 p[j][num[j]++]=i;                 
24.             }             
25.         }         
26.         eule[i]+=eule[i-1];   
27.     }  
28. }