原根

单求原根个数phi[phi[n]].

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
typedef long long ll;
const int maxn = 1e5;
bool flag[maxn]; //标记数组
ll phi[maxn]; //欧拉函数值
int prime[maxn]; //同时得到素数筛
int cnt = 0;
void Get_phi(int n)
{
    cnt = 0;
    memset(flag,true,sizeof(flag));
    phi[1] = 1;
    for(int i=2;i<=n;i++)
    {
        if(flag[i]) //素数
        {
            prime[cnt++] = i;
            phi[i] = i-1; //素数的欧拉函数值是i-1
        }
        for(int j=0;j<cnt;j++)
        {
            if(i*prime[j]>n)
            {
                break;
            }
            flag[i*prime[j]] = false;//素数的倍数不是素数
            if(i%prime[j]==0) //i%mod prime = 0,那么phi(i*p) = p*phi(i)
            {
                phi[i*prime[j]] = prime[j]*phi[i];
                break;
            }
            else phi[i*prime[j]] = (prime[j]-1)*phi[i];//i mod prime != 0, 那么 phi(i * prime) == phi(i) * (prime-1)
        }
    }
}
int main()
{
    ll n;
    Get_phi(maxn-10);
    while(cin>>n)
    {
        cout<<phi[phi[n]]<<endl;
    }
    return 0;

求具体原根，一步一步来。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
typedef long long ll;
const int maxn = 1e5;
//欧拉函数
ll euler(ll n){ //返回euler(n)
     ll res=n,a=n;
     for(ll i=2;i*i<=a;i++){
         if(a%i==0){
             res=res/i*(i-1);//先进行除法是为了防止中间数据的溢出
             while(a%i==0) a/=i;
         }
     }
     if(a>1) res=res/a*(a-1);
     return res;
}
//线性素数筛
int prime[maxn] = {0},num_prime = 0;
int isNotPrime[maxn] = {1,1};
void is_prime(int N)
{
    for(int i=2;i<N;i++)
    {
        if(!isNotPrime[i])
        {
            prime[num_prime++] = i;
        }
        for(int j=0;j<num_prime&&i*prime[j]<N;j++)
        {
            isNotPrime[i*prime[j]] = 1;
            if(!(i%prime[j]))
            {
                break;
            }
        }
    }
    return;
}
//分解质因数
int curprime[maxn];
int cnt = 0;
void dividep(int n)
{
    int t = sqrt(1.0*n);
    for(int i=0;prime[i]<=t;i++)
    {
        if(n%prime[i]==0)
        {
            curprime[cnt++] = prime[i];
            while(n%prime[i]==0) n /= prime[i];
        }
    }
    if(n>1) curprime[cnt++] = n;
}
int quick_pow(int a,int n,int mod)
{
    int res = 1;
    while(n)
    {
        if(n&1) res = res*a%mod;
        n /= 2;
        a = a*a%mod;
    }
    return res;
}
int main()
{
    ll n;
    is_prime(maxn-10);
    /*for(int i=0;i<num_prime;i++)
    {
        printf("%d ",prime[i]);
        if(num_prime%10==0)
        {
            printf("
");
        }
    }*/
    while(cin>>n)
    {
        cnt = 0;
        dividep(n-1);
        int ans = 0;
        for(int i=2;i<n;i++)
        {
            int flag = 0;
            for(int j=0;j<cnt;j++)
            {
                int t = (n-1)/curprime[j];
                if(quick_pow(i,t,n)==1)
                {
                    flag = 1;
                    break;
                }
            }
            if(!flag)
            {
                ans++;
            }
        }
        cout<<ans<<endl;
    }
    return 0;
}