原根和阶 及模板

原根和阶及模板

小记：

摆出定义和算法 $\Rightarrow$ 解释算法

主要参考 算法导论第31章 和传送门1

原根的定义

若 $or{d}_n(g)=\varphi (n)$, 即g对模n的阶等于n的欧拉函数，则g是n的原根。

如果g是n的原根，则对于$Z_n^{\ast }$的任意一个元素，都存在z，使得$g^z\equiv a(modn)$

什么样的N有原根

如果N=2，4，$p^e，2p^e$（p是大于2的质数，e是正整数），则N有原根

原根的意义？

如果$Z_n^{\ast }$包含一个原根，就称$Z_n^{\ast }$是循环的。

求N的原根的算法步骤

step1: 将N的欧拉函数phi进行质因数分解，phi(N)=$\varphi (N)$=$ p_1^{a_1} *p_2^{{a_2}*……p_n}{a_n}$

step2: 枚举i，判断每个i是否有$i^k$ mod N $==1的最小k是不是$$\frac{\varphi (N)}{p_i}$ ,符合条件的就是N的最小原根

step3：利用最小原根得到 其他原根，

​ 其他原根的形式都为 $g^k$形式，且gcd（k,$\varphi (N)$)=1

基础知识 🥛

---

1.模运算

$Z_n$={0,1,……，n-1}

$Z_n^{\ast }$={$a\in Z_n$ 且 gcd(a,n)=1}​

比如$Z_{15}^{\ast }$={1, 2, 4, 7, 8, 11, 13, 14}

2.$\varphi (n)$是什么?

$Z_n^{\ast }$的规模就是欧拉函数，$|Z_n^{\ast }|=\varphi (n)=n\ast {\prod }_{p_i是素数且p_i|n}(1-\frac{1}{p})$

比如 $\varphi （15）=8$

3.ord（）是什么？

a通过（自定义运算 a^k%n，k为正整数）得到有限群。

一个元素的阶等于它所生成子群的规模，即 ord(a)=|<a>|

也是 满足$a^k\%n=1$的最小正整数k。

举例如下

n=7时，$Z_7^{\ast }$={1,2,3,4,5,6}

n=7	0	1	2	3	4	5	6
$2^i\%n$	1	2	4	2	4	1	2
$3^i\%n$	1	3	2	6	4	5	1

如 n=7,i=2时，<a>={1,2,4}, $a^3=1$, |a|=3

​ n=7, i=3时，<a>={1,3,2,6,4,5}, $a^6=1$, |a|=6

​ 发现了吗？这时<a>与$Z_7^{\ast }$同，故3是7的原根

4.欧拉定理

如果n是正整数且$a\in Z_n^{\ast }$, 那么$a^{\varphi （n)}\%n=1$。

算法解释

那么回到原根定义，$or{d}_n(g)=\varphi (n)$

欧拉函数很好求，而g要怎么寻找呢？

我们采用枚举i的方式，对每个i，

看使它为 i^k %n=1的最小正整数k是不是$\varphi （n）$

并且如果$i\in Z_n^{\ast }$ ,那么由 欧拉定理 得到$ i^{\phi（n）}%n=1$，如果不属于的话，我们还要去计算证明一下。

那难道要用$i^k$，k$\in$ [1,$\varphi (n)$) 去检查吗？

不用，只需要枚举$i^{\frac{\varphi (N)}{p_i}}$ %n去就行，

原因是k 分为2部分，不能被$\varphi (n)$整除的，和能被$\varphi (n)$整除的。

对于不能被$\varphi (n)$整除部分

关于不能被$\varphi (n)$整除，

使用结论，若$i^g\%mod=1$,则对于任意不满足c | g的数不存在$i^c\%mod=1$，可以直接pass。

证明：

$使用反证法$
$假设存在不满足c|g使i^c\%mod=1,且c_1是最小的c$
$\therefore i^{c_1}\%mod=1,c_2=g-c_1$
$\because i^g\%mod=1,我们在前面已经确认过$
$\therefore (i^{c_1}\%mod)\ast (i^{c_2}\%mod)=1$
$\therefore i^{c_2}\%mod=1$
如果 c_1|c_2,则$g=c_1+c_2=(1+f)c_1$,与假设条件相矛盾
$\therefore gcd(c1,c2)<c_1$
又$\because$ $i^{c_2-c_1}\equiv i^{c_2}\ast (i^{c_1}{)}^{-1}\equiv 1(\%mod)$ $\therefore$由gcd的更相减损术原理得$i^{gcd(c_2,c_1)}$%mod=1
$又\because gcd(c1,c2)<c1与原假设c_1是最小相矛盾$
$\therefore 假设不成立，故若i^g\%mod=1，不满足c|g，则不能使得i^c\%mod=1$

对于能被$\varphi (n)$整除部分

把$\varphi （n）$进行质因数分解，得到 $\varphi （n)=p_1^{a_1}\ast p_2^{a_2}\ast \dots \dots$

只要 $i^{p_d}\%n\ne 1$,

那么如果 $\underset{\frac{\varphi (n)}{p_d}}{\underbrace{i^{p_d}\ast \dots \dots i^{p_d}}}\ne 1$,即只要证 $i^{\frac{\varphi (N)}{p_d}}$ %n$\neq$1即可

代码

(hdu4492模板）

#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);
            while (x%i==0)x/=i;
        }
    if(x>1)ret=ret/x*(x-1);
    return ret;
}
void getfac(int x) {
    //得到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 is_have_root(int x) {
    //如果N=2，4，p^e ，2*(p^e)（p是大于2的质数，e是正整数），则N有原根
    if(x==2||x==4)return 1;
    if(x%4==0)return 0;
    if(x%2==0)x/=2;
    for(int i=3; i*i<=x; ++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);        //得到x的欧拉函数
    getfac(phix);            //得到phi（x）的质因数
    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() {
//    freopen("in.text","r",stdin);
    while(scanf("%d",&n)==1) {
        if(!is_have_root(n))puts("-1");
        else {
            solve(n);
            for(int i=0; i<ans.size(); ++i)printf("%d%c",ans[i]," \n"[i==ans.size()-1]);
        }
    }
    return 0;
}