[欧拉函数]求原根个数

题面描述

【输入格式】
第一行n(1<=n<=10000)
下来n行，每行一个整数p，p保证为质数，表示求p的原根的数目。(1<=p<=2000 0000)
【输出格式】
每次询问输出一行一个整数，即x的原根的数目。
【样例输入】
3
23
31
79
【样例输出】
10
8
24

思路

原根个数是phi[phi[i]]个，我也不知道为什么？
caioj1183代码中的vv数组，其实就是判断是否满足原根条件的数。

代码

caioj1159

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
using namespace std;
const int N=2e7+10;
int a[N];
inline int solve(int x)
{
    int d=x,dd=x;
    for(int i=2;i*i<=x;i++)
    {
        if(x%i==0)
        {
            d=d/i*(i-1);
            while(!(x%i))
            {
                x/=i;
            }
        }
    }
    int last=0;
    if(x>1)d=d/x*(x-1);a[dd]=d;
    return d;
}
int main()
{
    int n;scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        int x;scanf("%d",&x);
        printf("%d\n",!a[x-1]?solve(x-1):a[x-1]);
    }
    return 0;
}

caioj1183

题面描述

给出一个n个正整数，求每个正整数的原根个数。
输入
第一行一个数n

接下来n行，一行一个数xi。

输出
一行一个数，第i行是xi的原根个数（无则输出0）。

样例输入
8
1
2
3
4
8
9
10
18
样例输出
0
1
1
1
0
2
2
2

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
using namespace std;
const int N=2e6+10;
const int M=1e6+10;
const int inf=2e6+10;
bool v[N],vv[N];int phi[N],prime[M],m;
inline void g_p()
{
	m=0;
	prime[++m]=2;phi[2]=1;v[4]=1;phi[4]=2;phi[1]=1;
	for(int i=3;i<=inf;i++)
	{
		if(!v[i])
		{
			prime[++m]=i;
			phi[i]=i-1;vv[i]=1;
		}
		for(int j=1;j<=m&&i*prime[j]<=inf;j++)
		{
			v[i*prime[j]]=1;
			if(i%prime[j]==0)
			{
				vv[i*prime[j]]=vv[i];
				phi[i*prime[j]]=phi[i]*prime[j];
				break;
			}
			else phi[i*prime[j]]=phi[i]*(prime[j]-1); 
		}
	}
	vv[2]=vv[4]=1;
}
int main()
{
	g_p();
	int n;scanf("%d",&n);
	for(int i=1;i<=n;i++)
	{
		int x;scanf("%d",&x);
		if(vv[x]||(!(x&1)&&vv[x>>1]&&x!=8))printf("%d\n",phi[phi[x]]);
		else puts("0");
	}
	return 0;
}