数论之阶与原根

---

一、阶（order）

1.定义

设$a$与$n$是互质的整数且$n>1$，使 $a^x\equiv 1(mod$ $n)$ 成立
的最小正整数$x$称为整数$a$模$n$的阶，记为:$or{d}_{n}a$

2.基本定理

定理一

内容:$gcd(a,n)==1$, $a^m\equiv 1(mod$ $n)$的充要条件为 $(or{d}_{n}a)|m$

证明

必要性：
已知$a^{or{d}_{n}a}\equiv 1(mod$ $n)$，则可以推出：$a^{k\ast or{d}_{n}a}\equiv 1(mod$ $n)$ $k为任意整数$ 。

充分性：
反证法，若存在$m=K\ast or{d}_{n}a+r（0<=r<or{d}_{n}a）$可使得： $a^m\equiv 1(mod$ $n)$，那么就意味着有： $a^r\equiv 1(mod$ $n)$，此时$r<or{d}_{n}a$，根据阶的定义，阶不应该是$or{d}_{n}a$而应该是$r$，故产生矛盾。

---

定理二

内容：$gcd(a,n)==1$,$or{d}_{n}a|\phi (n)$

证明

根据欧拉定理，此时有$a^{\phi (n)}\equiv 1(mod$ $n)$,再利用定理一，就可以证明：$or{d}_{n}a|\phi (n)$

---

定理三

内容:$gcd(a,n)==1$,$a^x\equiv a^y(mod$ $n)$的充要条件为$x\equiv y(mod$ $or{d}_{n}a)$

证明

先令$or{d}_{n}a$为$d$，有$a^d\equiv 1(mod$ $n)$

必要性：
条件转化为$x=k\ast d+y$
$a^x=a^{k\ast d+y}\equiv a^y(mod$ $n)$

充分性：
由于$gcd(a,n)==1$，使用消去律，得：$a^{x-y}\equiv 1(mod$ $n)$
由定理一可以推出$d|(x-y)\Rightarrow x-y=kd\Rightarrow x\equiv y(mod$ $or{d}_{n}a)$

---

定理四

内容:若$u$为正整数，则$or{d}_n{a}^u=\frac{or{d}_{n}a}{gcd(or{d}_{n}a,u)}$

证明

先令$or{d}_{n}a$为$d$，$t$为$gcd(d,u)$
那么就可以写为:$d=pt,u=qt,(p与q互质)$

已知$a^u$与$n$互质，根据定义求$or{d}_n{a}^u$也就是求满足:$(a^u{)}^x\equiv 1(mod$ $n)$的最小正整数$x$
转化一下得：$a^{qtx}\equiv 1(mod$ $n)$,由定理一可知，$d|qtx\Rightarrow pt|qtx\Rightarrow p|qx$
那么满足整除条件的最小的$x$就是p
所以有:$or{d}_n{a}^u=\frac{or{d}_{n}a}{gcd(or{d}_{n}a,u)}$

---

二、原根

1.定义

设$a$与$n$是互质的整数且$n>1$，则当$or{d}_{n}a=\phi (n)$时，称$a$为模$n$的原根。
即对于$a$来说，$\phi (n)$是满足$a^x\equiv 1(mod$ $n)$的最小正整数。

已知$a$为模$n$的原根，若有$a^x\equiv 1(mod$ $n)$，则$x$为$\phi (n)$的倍数。

2.性质

性质一（哪些数有原根）

内容：模$n$有原根的充要条件是n=2、4、$p^k$、$2p^k$，其中$p$为奇素数，$k$为正整数。

证明可见此处

---

性质二（若有原根，有几个不同余的原根）

内容：当正整数$n$有原根时，有$\phi (\phi (n))$个原根。

证明

根据阶中的定理四：若$u$为正整数，则$or{d}_n{a}^u=\frac{or{d}_{n}a}{gcd(or{d}_{n}a,u)}$
若$a$为原根，有$or{d}_{n}a=\phi (n)$

两式合并，得$or{d}_n{a}^u=\frac{\phi (n)}{gcd(\phi (n),u)}$
那么当$gcd(\phi (n),u)=1$,则有$or{d}_n{a}^u=\phi (n)$,

也就是当$gcd(\phi (n),u)=1$,$a^u$为模$n$的原根。

由于1与正整数都互质，所以得到结论：
当正整数$n$有原根时，有$\phi (\phi (n))$个原根。

推广（已知最小原根，求其他原根）

当我们已经得知了最小原根为$a$的情况下，想要求出其他不同余的原根，就只需要找出在$[1,\phi (n))$内与$\phi (n)$互质的数，将其作为指数，再取模求出另外的原根。

例：当$n=9$，已知最小原根为2

φ(9)=6，φ(6)=2，说明有两个原根，与6互质的数分别为1跟5:
$2^1(mod$ $9)=2$
$2^5(mod$ $9)=5$

所以原根分别为2跟5.

---

性质三

内容：若$a$为模$n$的原根，那么$a^1,a^2\dots \dots a^{\phi (n)}$构成了模$n$的既约剩余系

证明

第一步，证明$a^1,a^2\dots \dots a^{\phi (n)}$模$n$的结果两两不同余。

第二步，证明$a^1,a^2\dots \dots a^{\phi (n)}$模$n$的结果与$n$为互质关系。

第一步证明
反证法：若存在$i,j(1\le i<j\le \phi (n))$满足$a^j\equiv a^i(mod$ $n)$
通过消去律可以得到：$a^{j-i}\equiv 1(mod$ $n)$
$(j-i)$小于$\phi (n)$，但题目的原根为$a$，说明满足$a^x\equiv 1(mod$ $n)$的最小正整数为$\phi (n)$,产生矛盾，故不存在。
也就是 $a^1,a^2\dots \dots a^{\phi (n)}$模$n$的结果两两不同余。

第二步证明
由条件可知，$a$与$n$为互质关系，设存在$a^s$，$s$为正整数，问$a^s$模$n$是否与$n$互质？
设$a^s\equiv b(mod$ $n)\Rightarrow$ $a^s=kn+b$
若$b$与$n$不互质，有质因子$t$,$n=ut,b=vt$
上式可改写为:$a^s=(ku+v)t$
那么就是说明在$a$中同样会包含该质因子$t$，若包含，则与条件的$a,n$互质相矛盾，故不包含，也就是 $a^s$(s为正整数)模$n$仍然与$n$保持互质关系。

---

3.如何求原根

方法

根据之前的分析，想要求模$n$的所有原根，就需要
先判断n的形式是否存在原根，
再求出最小原根，
最后逐一推出其余的原根。

原根需要通过互质进行枚举，但对于原根的判断有点麻烦，所以可以直接利用公式来判断（没有证明）
将$\phi (n)$唯一分解为$p_1^{k_1}{p}_2^{k_2}\dots \dots p_m^{k_m}$的形式，
如果对于每一个质因子来说，都不满足：$a^{\frac{\phi (n)}{p_i}}\equiv 1(mod$ $n)$
那么就可以说明a是模n的原根。

bool check_minroot(ll a,vector<ll> tmp,ll n){//判断是否为原根 

	int sz=tmp.size();//tmp存的是将euler[n]唯一分解的质因子
	for(int i=0;i<sz;i++){
		if(Pow(a,euler[n]/tmp[i],n)==1)return 0;
	}
	return 1;
}

例题·Primitive Roots

Primitive Roots

代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll a[100005], len;
const int N=1e6+10;
int primes[N], euler[N], cnt;
bool st[N];

void get_eulers(){ //欧拉筛的扩展，同时求欧拉函数 
	int n=N;
	euler[1] = 1;
	st[0]=st[1]=1;
	for (int i=2; i <= n; i++ ){
		if ( !st[i] ){
			primes[cnt++] = i;
			euler[i] = i-1;
		}
		for ( int j=0; primes[j] <= n/i; j++ ){
			st[ primes[j]*i ] = true;
			if ( i % primes[j] == 0 ) {
				euler[ i*primes[j] ] = euler[i] * primes[j];
				break;
			}
			euler[ i*primes[j] ] = euler[i] * ( primes[j]-1 );
		}
	}
}

ll Pow(ll x,ll y,ll mod){//带模快速幂 
	ll ans=1;
	while(y){
		if(y&1){
			ans=(ans*x)%mod;
		}
		y>>=1;
		x=(x*x)%mod;//更新底数 
	}
	return ans;
}

ll gcd(ll a,ll b){
	return b?gcd(b,a%b):a;
}

vector<ll> divide(ll n){//分解质因数 
    vector<ll> ans;
    ll tmp=sqrt(n+0.5);
    for (ll i=0; primes[i]<=tmp;i++){
        if (n%primes[i]==0){
        	ans.push_back(primes[i]); 
            while(n%primes[i]==0)
                n/=primes[i];
        }
    }
    if (n>1)
        ans.push_back(n);
    return ans;
}

bool check(ll n){
	
	if(n%2==0)n/=2;
    if(!st[n])return true;
    if(n%2==0)return false;
    
    for(ll i=1; i<cnt; i++){//判断奇质数，所以要从3开始 
        if(n%primes[i]==0){
            while(n%primes[i]==0){
                n/=primes[i];
            }
            return n==1;
        }
    }
    return false;
} 

bool check_minroot(ll a,vector<ll> tmp,ll n){//判断是否为最小原根 

	int sz=tmp.size();
	for(int i=0;i<sz;i++){
		if(Pow(a,euler[n]/tmp[i],n)==1)return 0;
	}
	return 1;
}
vector<ll> primitive_root(ll n){
	vector<ll> ans;//存原根 
	vector<ll> tmp;//存对φ(n)唯一分解的质因数 
	//判断有无原根 
	if(n==2){
		ans.push_back(1);
		return ans; 
	}
	else if(n==4){
		ans.push_back(3);
		return ans;
	}
	if(!check(n)){
		return ans;	
	}

	tmp=divide(euler[n]);

	for(ll i=2;i<n;i++){//找最小原根 
		if(gcd(i,n)!=1)continue;
		if(check_minroot(i,tmp,n)){
			ans.push_back(i);
			break;
		}
	} 
	
	//计算其他的原根 
	ll temp=ans[0];
	for(int i=2;i<euler[n];i++){
		temp=temp*ans[0]%n;
		if(gcd(euler[n],i)==1){
			ans.push_back(temp);
		}
	}
	sort(ans.begin(),ans.end());
	return ans; 
}

int main(){
	get_eulers();
	ll n;
	while(~scanf("%lld",&n)){
		vector<ll>ans;
		
		ans=primitive_root(n);
		
		ll sz=ans.size();
		if(sz==0)cout<<"-1";		
		for(int i=0;i<sz;i++){
			cout<<ans[i];
			if(i<sz-1)cout<<" ";
		}
		cout<<endl;
	}	
}

4.简单应用

当n为素数时，求解同余方程$x^5\equiv a(mod$ $n)$

解：
当a与n互质时，若有原根g，则必定可以用$g^y\equiv a(mod$ $n)$,再使用BSGS求出满足条件的y使得$X^5=g^y$有解，这个解就是原方程的解。

具体应用见于数论之指标的讲解：

---