当n为质数时就变成了费马小定理
#include<iostream>
using namespace std;
typedef long long LL;
LL phi(int x){
LL res=x;
for(int i=2;i<=x/i;i++){
if(x%i==0){
res=res/i*(i-1);
while(x%i==0) x=x/i; //欧拉函数的个数与质因数的指数无关
}
}
if(x>1) res=res/x*(x-1);
return res;
}
int main(){
int n;
scanf("%d",&n);
while(n--){
int x;
scanf("%d",&x);
cout<<phi(x)<<endl;
}
return 0;
}
#include<iostream>
using namespace std;
typedef long long LL;
const int N=1000010;
int prime[N],euler[N],cnt;
bool st[N];
int get_euler(int n){
euler[1]=1;
for(int i=2;i<=n;i++){
if(!st[i]){
prime[cnt++]=i;
euler[i]=i-1;
}
for(int j=0;prime[j]<=n/i;i++){
st[prime[j]*i]=true;
if(i%prime[j]==0){
euler[prime[j]*i]=prime[j]*euler[i];
break;
}
euler[prime[j]*i]=(prime[j]-1)*euler[i];
}
}
}
int main(){
int n;
scanf("%d",&n);
get_euler(n);
LL res=0;
for(int i=1;i<=n;i++) res+=euler[i];
cout<<res<<endl;
return 0;
}
当i%pj==0时 例如i=8 pj=2时 i的质因子肯定包括pj 只需用i的质因子就行 在i里已经乘过(1-1/pj) (pj存的时质数)
当i%pj!=0时 例如i=9 pj=2时 i的质因子肯定不包括pj本身 需要再乘一下(1-1/pj) (pj存的是质数)