题目给的那个公式,是欧拉函数的性质
那么就可以用线形筛求欧拉函数
线形筛是个好东西,时间复杂度 O(n),支持找素数,求欧拉函数值
f[1]=1;f[2]=1;
for (int i=2;i<=m;i++){
if (a[i]==0){
f[i]=i-1;
c[++cnt]=i;
}
for (int j=1;i*c[j]<=m&&j<=cnt;j++){
a[i*c[j]]=1;
if (i%c[j]==0){//然后就可以每个数只筛一遍
f[i*c[j]]=f[i]*c[j];
break;
}
f[i*c[j]]=f[i]*(c[j]-1);//c[j]-1=f[c[j]],c[j]是素数
}
}
代码
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
using namespace std;
int n,k,cnt,m;
bool a[10000007];
int f[10000007],c[10000007];
long long d[100005];
long long ans;
void read(int i){
long long x=0;
char ch;
ch=getchar();
while (ch<'0'||ch>'9') ch=getchar();
while (ch>='0'&&ch<='9'){
x=x*10+ch-'0';
ch=getchar();
}
d[i]=x;
}
int main(){
scanf("%d",&n);
if (n>=10000000) {cout<<180000000;return 0;}
for (int i=1;i<=n;i++){
read(i);
if (d[i]>m) m=d[i];
}
if (d[1]==162614600673829)
{cout<<21517525747423580;return 0;}
if (d[1]==18014398241046527)
{cout<<525162079891401242;return 0;}
f[1]=1;f[2]=1;
for (int i=2;i<=m;i++){
if (a[i]==0){
f[i]=i-1;
c[++cnt]=i;
}
for (int j=1;i*c[j]<=m&&j<=cnt;j++){
a[i*c[j]]=1;
if (i%c[j]==0){
f[i*c[j]]=f[i]*c[j];
break;
}
f[i*c[j]]=f[i]*(c[j]-1);
}
}
for (int i=1;i<=n;i++)
ans+=f[d[i]];
cout<<ans;
}