作用:同标题
参考自:https://blog.csdn.net/liuzibujian/article/details/81086324?tdsourcetag=s_pcqq_aiomsg
代码:
void euler(int n)
{
phi[1]=1;///1要特判
for (int i=2;i<=n;i++)
{
if (flag[i]==0)
{
prime[++num]=i;
phi[i]=i-1;
}
for (int j=1;j<=num&&prime[j]*i<=n;j++)
{
flag[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]];///直接利用f(m*n)=f(m)*f(n)
}
}
}
利用到的结论:
①设f(n)为小于n的与n互质的个数,若m,n互质,f(n*m) = f(n)*f(m) (对应代码里else后的求法)
②p是n的质因子,f(n*p) = p * f(n)(对应代码里break上面的求法)