//欧拉函数:求不超过n且与n互素的正整数的个数
long long phi(long long n)
{
long long rea = n;
for (int i=2; i*i<=n; i++)
{
if (n%i==0)
{
rea-=rea/i;
do
n/=i;
while (n%i==0);
}
}
if (n>1)
rea-=rea/n;
return rea;
}
先把素数筛选出来,方便欧拉函数使用
bool boo[50000];
long long p[20000];
void prim()//线性筛素数,在main函数里面先调用函数prim
{
memset(boo, 0, sizeof (boo));
boo[0]=boo[1]=1;
int k=0;
for (int i=2; i<50000; i++)
{
if (!boo[i]) p[k++]=i;
for (int j=0; j<k&&i*p[j]<50000; j++)
{
boo[i*p[j]]=1;
if (!(i%p[j])) break;
}
}
}
long long phi(long long n)
{
long long rea = n;
for (int i=0; p[i]*p[i]<=n; i++)//对于一些不是素数的可以不用遍历
{
if (n%p[i]==0)
{
rea-=rea/p[i];
do
n/=p[i];
while (n%p[i]==0);
}
}
if (n>1)
rea-=rea/n;
return rea;
}