pub/bbs-gen.c: Carefully generate numbers of the correct sizes.

diff --git a/pub/bbs-gen.c b/pub/bbs-gen.c
index d2b4b6e5be4f979ac3c9cf8eab91de51904894a4..4796cf5c719942c0bb65f95133a1aa8a05249b48 100644 (file)
--- a/pub/bbs-gen.c
+++ b/pub/bbs-gen.c
@@ -65,7 +65,7 @@ int bbs_gen(bbs_priv *bp, unsigned nbits, grand *r, unsigned n,
   pgen_jumpctx j;
   pgen_gcdstepctx g;
   unsigned nb = nbits/2;
-  mp *x = MP_NEWSEC;
+  mp *x = MP_NEWSEC, *t = MP_NEW;
 
   /* --- Generate @p@ --- */
 
@@ -88,11 +88,15 @@ int bbs_gen(bbs_priv *bp, unsigned nbits, grand *r, unsigned n,
   g.r = mp_lsr(MP_NEW, bp->p, 1);
   g.g = MP_NEW;
   g.max = MP_ONE;
+  t = mp_lsl(t, MP_ONE, nbits - 1);
+  mp_div(&t, 0, t, bp->p);
+  if (MP_CMP(x, <, t)) x = mp_leastcongruent(x, t, x, g.jp.m);
   bp->q = pgen("q", MP_NEWSEC, x, event, ectx, n, pgen_gcdstep, &g,
               rabin_iters(nb), pgen_test, &rb);
   pfilt_destroy(&g.jp);
   mp_drop(g.r);
   mp_drop(g.g);
+  mp_drop(t);
   if (!bp->q) goto fail_q;
 
   /* --- Compute @n@ --- */