构造一个大素数条件下的本原元（JAVA实现）

原理

设p为安全素数，即使p=2q+1，且q为素数。由Fermat定理知g^(p-1) mod p=1，
即g^(2q) mod p=1，因而如果
w=min{t>0: g^t mod p=1}
则有w整除p-1=2q，因而由q是素数知，w只能是2或q，此时g是本原元等价于
g^2 mod p !=1 且 g ^q mod p !=1

编程步骤

1.随机生成一个素数q
2.令p=2q+1，并判断p是否是一个素数
3.如果是，那么p就是一个安全素数，否则转1
4.在安全素数p下，随机选取g，1<g<p-1.
5.如果g^2 mod p!=1 且 g ^ q mod p!=1，那么g就是本原元

代码

public static BigInteger generator() {
		Random rnd=new Random();
		BigInteger p,q,g,x,y;
		while(true) {
			q=BigInteger.probablePrime(1024, rnd);
			p=q.multiply(BigInteger.valueOf(2)).add(BigInteger.valueOf(1));
			if(p.isProbablePrime(20)) break;
		}
		
		while(true) {
			g=BigInteger.probablePrime(512, rnd);
			x=g.multiply(g).mod(p);
			y=pow(g,q,p);
			if(x.equals(BigInteger.valueOf(1))==false && y.equals(BigInteger.valueOf(1))==false) {
				break;
			}
		}
		return g;
	}