欧拉函数(转载)

这位博主写的很好，我只是搬运工

然后是板子

//欧拉筛求phi 
void make_phi()
{
	phi[1]=1;
	for(int i=2;i<=n;i++)
	{
		if(!vis[i])
		{
			prime[++cnt]=i;
			phi[i]=i-1;
		}
		for(int j=1;j<=num&&prime[j]*i<=n;j++)
		{
			vis[i*prime[j]]=1;
			if(i%prime[j]==0)
			{
				phi[i*prime[j]]=phi[i]*prime[j]; 
				break;
			}
			else	phi[i*prime[j]]=phi[i]*phi[prime[j]];//欧拉函数性质(i和prime[j]互质) 
                       //另一种写法: else phi[ i*prime[j] ]=phi[i]*(prime[j]-1);
		}
	}
}
void euler()
{
	for(int i=1;i<=n;i++)	phi[i]=i;
	for(int i=2;i<=n;i++)
	{
		if(phi[i]==i)//说明i是素数 
		{
			for(int j=i;j<=n;j+=i)
				phi[j]=phi[j]/i*(i-1);
		}
	}
}