【算法讲17：原根】阶 和 原根

$\mathfrak{S}\mathfrak{o}\mathfrak{u}\mathfrak{r}\mathfrak{c}\mathfrak{e}$

• 《初等数论及其应用》第六版第九章
  省略了很多帮助理解的例子和证明。浓缩了一下密集的知识点。

前置

• 乘性函数
  可以见博客的算法讲4、5、7

整数的 $\lceil$阶$\rfloor$ 和 $\lceil$原根$\rfloor$

整数的阶

• 定义：
  设 $a、n$ 为互素的整数，$a\ne 0$，$n>1$，使得 $a^x\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}n)$ 的最小正整数 $x$ 称为 $a$ 模 $n$ 的阶。记作 $or{d}_{n}a$
• 定理 $9.1$
  如果 $a\perp n$ 且 $a\ne 0$，$n>0$，那么正整数 $x$ 是 $a^x\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}n)$ 的一个解当且仅当 $or{d}_{n}a|x$
• 推论 $9.1.1$
  如果 $a\perp n$ 且 $a\ne 0$，$n>0$，那么 $or{d}_{n}a|\varphi (n)$
  证明：根据欧拉定理得到 $a^{\varphi (n)}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}n)$，根据定理 $9.1$得到。
  根据该推论，我们可以从 $\varphi (n)$的所有因子中，去更快枚举阶。
• 定理 $9.2$
  如果 $a\perp n$ 且 $a\ne 0$，$n>0$,那么 $a^i\equiv a^j(\mathrm{m}\mathrm{o}\mathrm{d}n)$ 当且仅当 $i\equiv j(\mathrm{m}\mathrm{o}\mathrm{d}or{d}_{n}a)$,其中 $i,j$ 为非负整数。

原根

• 定义：
  如果 $r\perp n$ ，$n>0$，当 $or{d}_{n}r=\varphi (n)$ 时，称 $r$ 是模 $n$ 的原根 或者称 $r$ 是 $n$ 的原根。
• 定理 $9.3$
  如果 $r\perp n$ ，$n>0$，且 $r$ 是 $n$ 的原根，那么下列整数：
  $$r^1,r^2,\cdots ,r^{\varphi (n)}$$
  构成了模 $n$ 的既约剩余系。
  证明：根据定理 $9.2$易证。
• 定理 $9.4$
  如果 $or{d}_{n}a=t$ 并且 $u$ 是一个正整数，那么有：
  $$or{d}_n(a^u)=t/(t,u)$$
  证明：令 $s=or{d}_n(a^u)$，$v=(t,u)$，$t=t_{1}v$ 且 $u=u_{1}v$，得到 $(t_1,u_1)=1$
  因为 $t_1=t/(t,u)$，所以要证明 $or{d}_n(a^u)=t_1$，所以先证明 $(a^u{)}^{t_1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}n)$。
  首先有
  $$(a^u{)}^{t_1}=(a^{u_{1}v}{)}^{(t/v)}=(a^t{)}^{u_1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}n)$$
  这是因为 $or{d}_{n}a=t$.由定理 $9.1$ 得到 $s|t_1$
  另一方面，由
  $$(a^u{)}^s=a^{us}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}n)$$
  得到 $t|us$，因此 $t_{1}v|u_{1}vs$，于是 $t_1|u_{1}s$
  由于 $(t_1,u_1)=1$ 故 $t_1|s$
  得到 $s=t_1=t/v=t/(t,u)$
• 推论 $9.4.1$
  令 $r$ 是模 $n$ 的原根,$n>1$，那么 $r^u$ 是模 $n$ 的一个原根当且仅当 $(u,\phi (m))=1$
• 定理 $9.5$
  如果正整数 $n$ 有一个原根，那么它一共有 $\phi (\phi (n))$ 个原根。
  证明：有 $\phi (\phi (n))$ 个这样的整数 $u$ 满足 $(u,\phi (n))=1$

原根的存在性

• 定理 $9.15$
  正整数 $n(n>1)$有原根当且仅当：
  $$n=2,4,p^t,2p^t$$
  其中 $p$ 为奇素数且 $t$ 是正整数。
  略去了一大片证明…

模板题

• 【模板】原根 | 洛谷 P6091
• 代码：
  时间复杂度：$O(n\log n+T(n^{0.25}\log n+\phi (n)\log \phi (n)))$
  比较复杂…

/*
 _            __   __          _          _
| |           \ \ / /         | |        (_)
| |__  _   _   \ V /__ _ _ __ | |     ___ _
| '_ \| | | |   \ // _` | '_ \| |    / _ \ |
| |_) | |_| |   | | (_| | | | | |___|  __/ |
|_.__/ \__, |   \_/\__,_|_| |_\_____/\___|_|
        __/ |
       |___/
*/
const int MAX = 1e6+50;

ll qpow(ll a,ll n){/* */ll res = 1LL;while(n){if(n&1)res=res*a%MOD;a=a*a%MOD;n>>=1;}return res;}
ll qpow(ll a,ll n,ll p){a%=p;ll res = 1LL;while(n){if(n&1)res=res*a%p;a=a*a%p;n>>=1;}return res;}

int cnt;
bool vis[MAX];
int prime[MAX];
int phi[MAX];
bool hasrt[MAX];
void init(int n){
    cnt = 0;
    phi[1] = 1;
    for(int i = 2;i <= n;++i){
        if(!vis[i]){
            prime[++cnt] = i;
            phi[i] = i - 1;
        }
        for(int j = 1;j <= cnt && i * prime[j] <= n;++j){
            vis[i * prime[j]] = 1;
            if(i % prime[j] != 0){
                phi[i * prime[j]] = phi[i] * (prime[j] - 1);
            }else{
                phi[i * prime[j]] = phi[i] * prime[j];
                break;
            }
        }
    }
    hasrt[2] = hasrt[4] = true;
    for(int i = 2;i <= cnt;++i){
        ll tmp = prime[i];
        for(int j = 1;j <= n && tmp <= n;++j){
            hasrt[tmp] = true;
            if(tmp * 2 <= n)hasrt[tmp*2] = true;
            tmp = tmp * prime[i];
        }
    }
}
int prm[MAX];
int get_factors(int ph){
    int shu = 0;
    for(int i = 2;i*i<=ph;++i){
        if(ph % i == 0){
            prm[++shu] = i;
            while(ph % i == 0)ph /= i;
        }
    }
    if(ph > 1)prm[++shu] = ph;
    return shu;
}
bool check(int n,int g,int shu){
    if(qpow(g,phi[n],n) != 1)return false;
    for(int i = 1;i <= shu;++i){
        if(qpow(g,phi[n]/prm[i],n) == 1)return false;
    }
    return true;
}
int get_minrt(int n,int shu){
    for(int i = 1;i < n;++i){
        if(check(n,i,shu))return i;
    }
    return 0;
}
int ans[MAX];		/// 原根存在这里
void get_rt(int n,int rt){
    int tmp = 1;
    int shu = 0;
    for(int i = 1;i <= phi[n];++i){
        tmp = (ll)tmp * rt % n;
        if(__gcd(i,phi[n]) == 1){
            ans[++shu] = tmp;
        }
    }
    return ;
}
int main()
{
    init(1000000);
    int T;scanf("%d",&T);
    while(T--){
        int n,d;scanf("%d%d",&n,&d);
        if(!hasrt[n])printf("0\n\n");
        else{
            int res = phi[phi[n]];			/// 原根个数
            printf("%d\n",res);
            int shu = get_factors(phi[n]);
            int rt = get_minrt(n,shu);
            get_rt(n,rt);
            sort(ans+1,ans+1+res);
            for(int i = 1;i <= res/d;++i)
                printf("%d ",ans[i*d]);
            puts("");
        }
    }
    return 0;
}