阶:al≡1a^l \equiv 1al≡1 (mod(mod(mod m)m)m) 该方程lll的最小正整数解成为aaa在模mmm意义下的阶,记为ordm(a)ord_m(a)ordm(a)
原根:设mmm是正整数,aaa是(正)整数,若aaa模mmm的阶等于φ(m)φ(m)φ(m),则称aaa为模mmm的一个原根。
原根的性质
1.设ggg为模ppp下的原根,原根保证了gl≡1g^{l} \equiv1gl≡1 (mod(mod(mod p)p)p)该方程最小正整数解为φ(p)\varphi(p)φ(p),即g1,g2...gφ(p)−1g^1,g^2...g^{\varphi(p)-1}g1,g2...gφ(p)−1互不同余
2.模mmm有原根的充要条件是m=1,2,4,p,2p,pnm= 1,2,4,p,2p,p^nm=1,2,4,p,2p,pn,其中ppp是奇质数,nnn是任意正整数。
3.当模mmm有原根时它有φ(φ(m))φ(φ(m))φ(φ(m))个原根
原根的求法
假定模数是mmm
首先将φ(m)\varphi(m)φ(m)进行素因数分解
φ(m)=p1c1∗p2c2∗...∗pncn\varphi(m)=p_1^{c_1}*p_2^{c_2}*...*p_n^{c_n}φ(m)=p1c1∗p2c2∗...∗pncn
从222到p−1p-1p−1枚举,如果:∀pi,gφ(p)pi≠1(mod p)\forall p_i,g^{\frac{\varphi(p)}{p_i}}\not=1(mod\ p)∀pi,gpiφ(p)=1(mod p)
(此处不等于就是不同余于),那么ggg就是ppp的一个原根
CodeCodeCode
void euler(int n){ //线性求欧拉函数
phi[1]=1;
for (register int i=2;i<=n;++i){
if (!vis[i]) prime[++tot]=i,phi[i]=i-1,vis[i]=i;
for (register int j=1;j<=tot;++j){
if (prime[j]>vis[i] || prime[j]*i>n) break;
vis[i*prime[j]]=prime[j];
phi[i*prime[j]]=phi[i]*(i%prime[j]? prime[j]-1 : prime[j]);
}
}
}
int PFD(int p,int *prime){ //质因数分解
int tot=0;
for (register int i=2;p!=1;++i){
if (p%i==0) prime[++tot]=i;
while (p%i==0) p/=i;
}
return tot;
}
int find_g(int p){ //求原根
static int prime[MAXM+10];
int tot=PFD(phi[p],prime),tmp;
for (register int g=1;g<p;++g){
tmp=1;
for (register int i=1;i<=tot;++i){
if (fastpow(g,phi[i]/prime[i],p)==1){
tmp=0;
break;
}
}
if (tmp) return g;
}
return -1;
}
离散对数
当模mmm有原根时,设ggg为模mmm的一个原根,则当:gl≡x(mod m)g^l\equiv x(mod\ m)gl≡x(mod m)
有:
Indgx≡l (mod φ(m))Ind_gx\equiv l\ (mod\ \varphi(m))Indgx≡l (mod φ(m))
把IndgxInd_gxIndgx成为底数为ggg,模φ(m)\varphi(m)φ(m)时的离散对数值
离散对数的性质:(和普通对数相似)
Indg(xy)≡Indg(x)+Indg(y) (mod φ(m))Ind_g(xy)\equiv Ind_g(x)+Ind_g(y)\ (mod\ \varphi(m))Indg(xy)≡Indg(x)+Indg(y) (mod φ(m))
Indg(xy)≡y∗Indg(x) (mod φ(m))Ind_g(x^y)\equiv y*Ind_g(x)\ (mod\ \varphi(m))Indg(xy)≡y∗Indg(x) (mod φ(m))
离散对数有什么用呢?
它可以把乘法转换成加法,所以利用它就可以将模mmm同余系下的乘法转化为模φ(m)\varphi(m)φ(m)同余系下的加法问题
只要构造映射,将xxx映射到IndgxInd_gxIndgx即可
由原根的性质一:该映射是满的单射(一一对应)
CodeCodeCode
g=math::find_g(m); //找原根
x=g;
for (register int i=1;i<phi[m]++i){
ind[x]=i; //开个桶建立映射关系
x=x*g%m;
}