HDU

模板题，可用于求一个数的所有原根。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e6+10,inf=0x3f3f3f3f;
int n,fac[N],nf;
vector<int> ans;
int Pow(int x,int p,int mod) {
    int ret=1;
    for(; p; p>>=1,x=(ll)x*x%mod)if(p&1)ret=(ll)ret*x%mod;
    return ret;
}
int phi(int x) {
    int ret=x;
    for(int i=2; i<=x/i; ++i)if(x%i==0) {
            ret=ret/i*(i-1);
            for(; x%i==0; x/=i);
        }
    if(x>1)ret=ret/x*(x-1);
    return ret;
}
void getfac(int x) {
    nf=0;
    for(int i=2; i<=x/i; ++i)if(x%i==0)
            for(fac[nf++]=i; x%i==0; x/=i);
    if(x>1)fac[nf++]=x;
}
bool check(int x) {
    if(x==2||x==4)return 1;
    if(x%4==0)return 0;
    if(x%2==0)x/=2;
    for(int i=3; i<=x/i; ++i)if(x%i==0) {
            for(; x%i==0; x/=i);
            return x==1;
        }
    return 1;
}
int minRoot(int x,int phix) {
    if(x==2||x==4)return x-1;
    for(int i=2; i<x; ++i)if(Pow(i,phix,x)==1) {
            bool flag=1;
            for(int j=0; j<nf; ++j)if(Pow(i,phix/fac[j],x)==1) {flag=0; break;}
            if(flag)return i;
        }
}
void solve(int x) {
    int phix=phi(x);
    getfac(phix);
    int r=minRoot(x,phix);
    ans.clear(),ans.push_back(r);
    for(int i=2; i<phix; ++i)if(__gcd(i,phix)==1)ans.push_back(Pow(r,i,x));
    sort(ans.begin(),ans.end());
}
int main() {
    while(scanf("%d",&n)==1) {
        if(!check(n))puts("-1");
        else {
            solve(n);
            for(int i=0; i<ans.size(); ++i)printf("%d%c",ans[i]," 
"[i==ans.size()-1]);
        }
    }
    return 0;
}