5 /* find a prime p of length n and a generator alpha of Z^*_p */
6 /* Alg 4.86 Menezes et al () Handbook, p.164 */
8 gensafeprime(mpint *p, mpint *alpha, int n, int accuracy)
14 genprime(q, n-1, accuracy);
16 mpadd(p, mpone, p); /* p = 2*q+1 */
17 if(probably_prime(p, accuracy))
20 /* now find a generator alpha of the multiplicative */
21 /* group Z*_p of order p-1=2q */
24 mprand(n, genrandom, alpha);
25 mpmod(alpha, p, alpha);
26 mpmul(alpha, alpha, b);
28 if(mpcmp(b, mpone) == 0)
30 mpexp(alpha, q, p, b);
31 if(mpcmp(b, mpone) != 0)