原根

原根的定义： 

原根，是一个数学符号。设m是正整数，a是整数，若a模m的阶等于φ(m)，则称a为模m的一个原根。

1. 原根的定义
设m是正整数，a是整数，若a模m的阶等于φ(m)，则称a为模m的一个原根。（其中φ(m)表示m的欧拉函数）
假设一个数g对于P来说是原根，那么g^i mod P的结果两两不同，且有 1<g<P, 0<i<P,那么g可以称为是P的一个原根,归根到底就是g^(P-1) = 1 (mod P)当且仅当指数为P-1的时候成立。(这里P是素数)。

简单来说，g^i mod p ≠ g^j mod p （p为素数）
其中i≠j且i，j介於1至(p-1)之间，则g为p的原根。
求原根目前的做法只能是从2开始枚举，然后暴力判断g^(P-1) = 1 (mod P)是否当且当指数为P-1的时候成立，而由于原根一般都不大，所以可以暴力得到。
2. 原根的性质
（1）可以证明，如果正整数(a,m) = 1和正整数 d 满足a^d≡1(mod 7)，则 d 整除 φ(m)。因此Ordm(a)整除φ(m)。在例子中，当a= 3时，我们仅需要验证 3 的 1 、2、3 和 6 次方模 7 的余数即可。
（2）记δ = Ordm(a)，则a^1，……a^(δ-1)模 m 两两不同余。因此当a是模m的原根时，a^0,a^1，……a^(δ-1)构成模 m 的简化剩余系。
（3）模m有原根的充要条件是m= 1,2,4,p,2p,p^n，其中p是奇质数，n是任意正整数。
（4）对正整数(a,m) = 1，如果 a 是模 m 的原根，那么 a 是整数模n乘法群（即加法群 Z/mZ的可逆元，也就是所有与 m 互素的正整数构成的等价类构成的乘法群）Zn的一个生成元。由于Zn有 φ(m)个元素，而它的生成元的个数就是它的可逆元个数，即 φ(φ(m))个，因此当模m有原根时，它有φ(φ(m))个原根。 

来一道题：

Primitive Roots

Time Limit: 1000MS		Memory Limit: 10000K
Total Submissions: 6380		Accepted: 3638

Description

We say that integer x, 0 < x < p, is a primitive root modulo odd prime p if and only if the set { (xi mod p) | 1 <= i <= p-1 } is equal to { 1, ..., p-1 }. For example, the consecutive powers of 3 modulo 7 are 3, 2, 6, 4, 5, 1, and thus 3 is a primitive root modulo 7. 
Write a program which given any odd prime 3 <= p < 65536 outputs the number of primitive roots modulo p. 

Input

Each line of the input contains an odd prime numbers p. Input is terminated by the end-of-file seperator.

Output

For each p, print a single number that gives the number of primitive roots in a single line.

Sample Input

23
31
79

Sample Output

10
8
24

题目大概意思：
给出一个素数 p(3≤p&lt;65536) p(3≤p&lt;65536)p(3≤p<65536) ，求模 p pp 的原根的个数。

 

分析：

根据以下两个定理即可解决此题：

1.模 m mm 有原根的充要条件是 m=2,4,pa,2pa m=2,4,p^a,2p^am=2,4,pa,2pa（a次方），其中 p 是奇素数， a  是任意正整数。
2.当模 m mm 有原根时，它有 ϕ(ϕ(m)) 个原根。

（其中 ϕ(x) \phi(x)ϕ(x) 为欧拉函数）

由于 p pp 为大于 2 22 的素数，因此一定有原根。我们预处理出 [1,65536] [1,65536][1,65536] 的欧拉函数值的表，再根据定理计算答案即可。
 

 

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int get_euler(int n)
{
    int res=n,a=n;
    for (int i=2;i*i<=a;i++)
    {
        if (a%i==0)
        {
            res=res/i*(i-1);
            while(a%i==0)
            a=a/i;
        }
    }
    if(a>1)
    res=res/a*(a-1);
    return res;
}
int main()
{
    int n;
    ios::sync_with_stdio(false);
    while(cin>>n)
    {
        cout<<get_euler(n-1)<<endl;
    }
    return 0;
}