1. 程式人生 > >歐拉函數線性篩法

歐拉函數線性篩法

style getc string nbsp rim log getchar 線性篩 tchar

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<cmath>
 5 #define lli long long int 
 6 using namespace std;
 7 const int MAXN=10000001;
 8 void read(int &n)
 9 {
10     char c=+;int x=0;bool flag=0;
11     while(c<0||c>9){c=getchar();if
(c==-)flag=1;} 12 while(c>=0&&c<=9) 13 x=(x<<1)+(x<<3)+c-48,c=getchar(); 14 flag==1?n=-x:n=x; 15 } 16 int n,m; 17 bool check[MAXN]; 18 int prime[MAXN]; 19 int phi[MAXN]; 20 int tot=0; 21 int main() 22 { 23 cin>>n; 24 phi[1
]=1; 25 for(int i=2;i<=n;i++) 26 { 27 if(!check[i]) 28 prime[++tot]=i,phi[i]=i-1;// 只有i與它互質 29 for(int j=1;j<=tot;j++) 30 { 31 if(i*prime[j]>n) 32 break; 33 check[i*prime[j]]=1
; 34 if(i%prime[j]==0) 35 { 36 phi[i*prime[j]]=phi[i]*prime[j]; 37 break; 38 } 39 else 40 phi[i*prime[j]]=phi[i]*(prime[j]-1); 41 } 42 } 43 printf("%d\n",phi[n]); 44 return 0; 45 }

註意 if(i%prime[j]==0)

不要寫成if(!i%prime[j])

歐拉函數線性篩法