ssl1232函数 线形筛

题目给的那个公式,是欧拉函数的性质
那么就可以用线形筛求欧拉函数
线形筛是个好东西,时间复杂度 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;
}