听灌多 为何ce
  • 板块灌水区
  • 楼主Harley_Wu
  • 当前回复10
  • 已保存回复10
  • 发布时间2024/7/16 14:33
  • 上次更新2024/7/16 15:23:15
查看原帖
听灌多 为何ce
471331
Harley_Wu楼主2024/7/16 14:33
#include<bits/stdc++.h>
using namespace std;
int n,phi[1000010],prime[1000010],isprime[1000010],cnt,total;
void eular(int n)
{
	for(int i=2;i<=n;i++)
	{
		if(!isprime[i])
		{ 
			prime[++cnt]=i; 
			phi[i]=i-1;
		}
		for(int j=1;j<=cnt;j++)
		{ 
			if(prime[j]*i>n) break;
			isprime[prime[j]*i]=1;
			if(i%prime[j]==0)
			{
				phi[i*prime[j]]=phi[i]*prime[j]; 
				break;
			}
			else phi[i*prime[j]]=phi[i]*(prime[j]-1);
		}
	}
}
void seive(int Maxx)
{ 
	memset(isprime,true,sizeof(isprime));
	isprime[0]=false;
	isprime[1]=false;
	for(int i=2;i<=Maxx;i++)
	{
		if(isprime[i]) prime[++total]=i;
		for(int j=1;j<=total&&i*prime[j]<=Maxx;j++)
		{ 
			isprime[i*prime[j]]=false;
			if(!(i%prime[j])) break;
		}
	}
}
int main()
{
    cin>>n;
    seive(n);
    int ans=0;
    for(int i=2;i<=n;i++)
    {
    	ans+=phi[i];
	}
	cout<<ans;
    return 0;
} 
2024/7/16 14:33
加载中...