DH（Diffie-Hellman）算法本元根(原根)求法

Diffie-Hellman 算法中，比较重要的是大素数的本元根（原根，Primitive Root）的求法。

具体算法原理及证明在《计算机程序设计艺术》卷2中有详细描述，

算法简要流程（伪码）如下：

------------------------------------------------------------------------------------

BIGINT g, p, q;

1.  生成一个大素数 q，直接p是素数，其中 p = 2q + 1

do{
    q = gen_prime();
    p = q * 2 + 1;
}while( ! is_prime(p) );

2.  生成一个随机数 g，1 < g < p - 1，直到g^2 mod p 和 g^q mod p 都不等于 1

do{

    g = random(2, p - 1); 

}while( (pow(g, 2) % p == 1) || (pow(q, 2) % p == 1) );

3. 得到g是p的本元根
print("p的本元根g是:" + g);