HDU 4992 / BC 8D Primitive Roots

题意：对于给定的n，求所有的x，满足x的y次方模n等于1，并且不存在比φ(n)更小的y   其中φ是欧拉函数

如果x1 满足以上性质 我们称x1为n的原根 那么n的原根之间有没有联系呢？

对于x1满足不存在更小的y使得 x1^y%n=1  那么当且仅当x1^(kφ(n))%n才能为1(k为任意正整数) 

对于x1²满足 x1²^φ(n)%n=x1^2φ(n)%n=1  并且2φ(n)之内只有x1^φ(n)是满足性质的 也就是说如果x1²不是原根，当且仅当x1²^y=x1^φ(n)

也就是说 对于x1  x1^k(k与φ(n)互质)均为n的原根

所以我们暴力(不能太暴力)求出一组x1后 再算出其他x即可

#include <iostream>
#include <cstring>
#include <string>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <vector>
#include <queue>
#include <map>
#include<set>
#include<stack>
#define bug puts("bugbugbug");
using namespace std;
typedef long long ll;
bool visit[1000010]= {0};
int prime[100000];
int  nump;
int init_prim(int n)
{
    memset(visit, true, sizeof(visit));
    int num = 0;
    for (int i = 2; i <= n; ++i)
    {
        if (visit[i] == true)
        {
            num++;
            prime[num] = i;
        }
        for (int j = 1; ((j <= num) && (i * prime[j] <= n));  ++j)
        {
            visit[i * prime[j]] = false;
            if (i % prime[j] == 0) break; //点睛之笔
        }
    }
    return num;//返回质数
}
ll pow_mod(ll a, ll n, ll m)
{
    ll res = 1;
    while(n)
    {
        if (n & 1) res = res * a % m;
        n >>= 1;
        a = a * a % m;
    }
    return res;
}
ll euler(ll x)
{
    ll res = x;
    for (ll i = 2; i <= x / i; i++) if (x % i == 0)
        {
            res = res / i * (i - 1);
            while(x % i == 0) x /= i;
        }
    if (x > 1) res = res / x * (x - 1);
    return res;
}
ll gcd(ll a,ll b)
{
    return b==0?a:gcd(b,a%b);
}
vector<int>ans;
vector<int>vec;
void go(ll n)
{
    ans.clear();
    vec.clear();
    ll k=euler(n);
    ll kk=k;
     for(ll i=1; i<nump&&k!=1; i++)
            if(k%prime[i]==0)
            {
                ll l=vec.size();
                ll cc=1;
                while(k%prime[i]==0&&k!=1)
                {
                    k/=prime[i];
                    cc*=prime[i];
                    vec.push_back(cc);
                    for(ll j=0; j<l; j++)
                        vec.push_back(vec[j]*cc);
                }
            }
            if(k!=1){
                int l=vec.size();
                for(int i=0;i<l;i++)
                {
                    vec.push_back(vec[i]*k);
                }
            vec.push_back(k);
            }
    ll a=0;
    for(int i=2; i<=n; i++)
        if(pow_mod(i,kk,n)==1)
        {
            int flag=1;
            for(int j=0; j<vec.size(); j++)
                if(vec[j]!=kk&&pow_mod(i,vec[j],n)==1)
                {
                    flag=0;
                    break;
                }
            if(flag)
            {
                a=i;
                break;
            }
        }
        if(a==0)
        {
            printf("-1\n");
            return;
        }
        memset(visit,0,sizeof(visit));
        ans.push_back(a);
        ll res=a;
        res=res*res%n;
        for(int i=2;i<kk;i++,res=res*a%n)
        if(gcd(i,kk)==1)
        {
            ans.push_back(res);
        }
        sort(ans.begin(),ans.end());
        for(int i=0;i<ans.size();i++)
        {
            printf("%d%c",ans[i],i==ans.size()-1?'\n':' ');
        }
}
int main()
{
    nump=init_prim(1000005);
    ll n;
    while(~scanf("%lld",&n))
    {
        if(n==2)
        {
            puts("1");
        }
        else if(n==4)
        {
            puts("3");
        }
        else
        {
            int nn=n;
            if(nn%2==0)nn/=2;
            for(int i=2; i<nump&&prime[i]<=nn; i++)
                if(nn%prime[i]==0)
                {
                    while(nn%prime[i]==0)
                        nn/=prime[i];
                    break;
                }
            if(nn==1) go(n);
            else puts("-1");
        }
    }
}