阶&原根

求阶的方法：

根据性质2，直接对ϕ(m)求出因子即可，从小到大依次判断是不是符合a^d = 1（mod m）（d是ϕ(m)的因子）

求最小的原根的方法：

根据性质8，对ϕ(m)求出素因子，从1开始不断测试即可，因为最小的原根很容易暴力得到。

求原根代码：（下面代码是求素数p的原根，如果不是素数，需要求出p的欧拉函数值，由于是素数，所以欧拉函数值为p-1）

 ll pow(ll a, ll b, ll m)
 {
     a %= m;
     ll ans = ;
     while(b)
     {
         if(b & )ans = ans * a % m;
         a = a * a % m;
         b /= ;
     }
     return ans % m;
 }
 int tot;//素因子个数
 int a[];
 void get_fact(ll n)//质因数分解n
 {
     for(ll i = ; i * i <= n; i++)
     {
         if(n % i == )
         {
             a[tot++] = i;
             while(n % i == )n /= i;
         }
     }
     if(n != )a[tot++] = n;
 }
 bool g_test(ll g, ll p)//测试g是不是p的原根 此处p是素数 欧拉函数值为p - 1
 {
     for(ll i = ; i < tot; i++)
     {
         if(pow(g, (p - ) / a[i], p) == )return false;
     }
     return true;
 }
 int proot(ll p)
 //求解p最小原根，本题p为素数
 //返回最小的原根
 {
     get_fact(p - );//素数的欧拉函数值为p - 1
     int g = ;
     while()
     {
         if(g_test(g, p))return g;
         g++;
     }
 }