欧拉函数（小于或等于n的数中与n互质的数的数目）&& 欧拉函数线性筛法
【欧拉函数】
在数论，对正整数n，欧拉函数是少于或等于n的数中与n互质的数的数目。此函数以其首名研究者欧拉命名，它又称为Euler’s totient function、φ函数、欧拉商数等。 例如φ(8)=4，因为1,3,5,7均和8互质。 从欧拉函数引伸出来在环论方面的事实和拉格朗日定理构成了欧拉定理的证明。
【证明】：
设A, B, C是跟m, n, mn互质的数的集，据中国剩余定理，A*B和C可建立一一对应的关系。因此φ(n)的值使用算术基本定理便知，
若
n= ∏p^(α（下标p）)p|n
则φ(n)=∏(p-1)p^(α（下标p）-1)=n∏(1-1/p)
p|n p|n
例如φ(72)=φ(2^3×3^2)=(2-1)2^(3-1)×(3-1)3^(2-1)=24，与欧拉定理、费马小定理的关系，对任何两个互质的正整数a, m, m>=2有a^φ(m)≡1(mod m)即欧拉定理：当m是质数p时，此式则为：a^(p-1)≡1(mod m)即费马小定理。（慢慢理解~~）
代码实现：（写一遍欧拉函数，加深印象！）
在线版：
1 #include <bits/stdc++.h> 2 using namespace std; 3 int eular(int n) 4 { 5 int res=1; 6 for(int i=2;i*i<=n;i++){ 7 if(n%i==0){ 8 n/=i,res*=i-1;//保证i一定是素数 9 while(n%i==0) 10 n/=i,res*=i; 11 } 12 } 13 if(n>1) 14 res*=n-1; 15 return res; 16 } 17 int main() 18 { 19 int n; 20 while(scanf("%d",&n)!=EOF){ 21 printf("%d\n",eular(n)); 22 } 23 return 0; 24 }
预处理：
1 #include <bits/stdc++.h> 2 using namespace std; 3 const int N=le5+5; 4 int phi[N]; 5 void pre_eular() 6 { 7 phi[1]=1; 8 for(int i=2; i<N; i++) 9 { 10 if(!phi[i]) 11 { 12 for(int j=i; j<N; j+=i) 13 { 14 if(!phi[j]) phi[j]=j; 15 phi[j]=phi[j]/i*(i-1); 16 } 17 } 18 } 19 }
欧拉函数的和：phi_sum(n) = the sum of phi(i) where gcd(i,n) = 1 and 1 <= i <= n
1）phi_sum(n) = n * phi(n) / 2 (n >= 2)
2）phi_sum(n) = 1 (n == 1)
线性筛：该算法在可在线性时间内筛素数的同时求出所有数的欧拉函数。
需要用到如下性质(p为质数)： 1. phi(p)=p-1 因为质数p除了1以外的因数只有p，故1至p的整数只有p与p不互质 2. 如果i mod p = 0, 那么phi(i * p)=p * phi(i) 证明如下
(上述证明存在bug。。感谢@PrimaryOIer指教) 上面的过程证明了从区间[1,i]->[i+1,i+i]，若整数n不与i互质，n+i依然与i不互质。下面给出另一个证明：若整数n与i互质，n+i与i依然互质
3.若i mod p ≠0, 那么phi(i * p)=phi(i) * (p-1) i mod p 不为0且p为质数, 所以i与p互质, 那么根据欧拉函数的积性phi(i * p)=phi(i) * phi(p) 其中phi(p)=p-1即第一条性质
1 #include<iostream> 2 #include<cstdio> 3 #define N 40000 4 using namespace std; 5 int n; 6 int phi[N+10],prime[N+10],tot,ans; //phi - 欧拉函数的值 ， prime - 素因子的值 7 bool mark[N+10]; 8 void getphi() 9 { 10 int i,j; 11 phi[1]=1; 12 for(i=2;i<=N;i++)//相当于分解质因式的逆过程 13 { 14 if(!mark[i]) 15 { 16 prime[++tot]=i;//筛素数的时候首先会判断i是否是素数。 17 phi[i]=i-1;//当 i 是素数时 phi[i]=i-1 18 } 19 for(j=1;j<=tot;j++) 20 { 21 if(i*prime[j]>N) break; 22 mark[i*prime[j]]=1;//确定i*prime[j]不是素数 23 if(i%prime[j]==0)//接着我们会看prime[j]是否是i的约数 24 { 25 phi[i*prime[j]]=phi[i]*prime[j];break; 26 } 27 else phi[i*prime[j]]=phi[i]*(prime[j]-1);//其实这里prime[j]-1就是phi[prime[j]]，利用了欧拉函数的积性 28 } 29 } 30 } 31 int main() 32 { 33 getphi(); 34 }
