The NewBitGenerator command outputs a Maple procedure, a pseudo-random bit generator, which when called outputs one pseudo-random bit, a 0 or 1. The generator is a Blum-Blum-Shub (BBS) generator. A BBS generator uses the following quadratic recurrence to generate a sequence of integers $x_1,x_2,\dots$ from which cryptographically secure pseudo-random bits $z_1,z_2,\dots$ are extracted:

where

$n$ is a product of two large primes p and q,

$m=⌊{\log }_2\left({\log }_2\left(n\right)\right)⌋$, and

$x_0$ is determined from the seed s.

Each iteration of $z_{k+1}$, the least significant $m$ bits of $x_{k+1}$, generates $m$ cryptographically secure bits. The output of the generator depends on the options described below. The default is to output one bit.

The cryptographic security of the BBS generator assumes that the number theoretic problem of distinguishing a quadratic residue from a pseudo-square in Z mod n is computationally infeasible when n is the product of two primes p and q and the factorization of n is not known. Thus it also assumes that integer factorization is computationally infeasible. Recall the definitions of a quadratic residue and pseudo-square:

Definition: An integer x in Z mod n is a quadratic residue if (i) $\gcd \left(x\,n\right)=1$ and (ii) $x=y^2\mathbf{mod}n$ for some integer y.

Definition: An integer z in Z mod n where $n=pq$ is a pseudo-square if (i) $\gcd \left(z\,n\right)=1$, (ii) z is not a quadratic residue in Z mod p, and (iii) z is not a quadratic residue in Z mod q.

The primes p and q are chosen to be of the form $p=2s+1$ and $q=2t+1$ where s and t are prime and 2 is a primitive element in both Z mod s and Z mod t.  This choice guarantees that p and q are both congruent to 3 mod 4 (a requirement for the security of the generator) and also that the period of the generator for $x_0$ for any quadratic residue other than 1 is either $s-1$ or $t-1$ or $\mathrm{lcm}\left(s-1\,t-1\right)$, all of which are large. Random primes p and q, and their product $n=pq$, satisfying these requirements of lengths 512 bits, 768 bits, and 1024 bits have been precomputed and the prime factorization discarded (Maple does not have the factorizations). Thus there are three choices for $n$ of lengths 308, 462, and 616 decimal digits, providing a range of security for cryptographic applications.

The argument seed=s determines the seed for the generator. The value used for $x_0$ used must be a "random" quadratic residue which is not equal to 1 (to avoid a short period). The value used for $x_0$ is computed from the seed s. For cryptographic applications the seed should be chosen by the user to be a random integer from a large set, e.g., from 0..10^100. It will be automatically reduced modulo n.

The following optional arguments are supported. They are input as equations in any order.

primes = $512$, $768$ or $1024$

The integer l, which must be one of $512$, $768$ or $1024$ specifies the length of the primes p and q in bits.  The default is $512$.

numbits=integer

This specifies how many bits are computed and output by each call to the generator.  The default is 1.

output= bits or integer

This specifies how the bits are output.  If integer is specified, the output is a Maple integer in the range 0 to 2^b-1.  If bits is specified, the output will be a Maple sequence of b bits.  The default is bits.

The RandomTools[BlumBlumShub] module also exports the NewGenerator function.  This function's interface is compatible with the NewGenerator functions of other RandomTools pseudo-random number generator subpackages.

$\mathrm{with}\left(\mathrm{RandomTools}\left[\mathrm{BlumBlumShub}\right]\right)$

$\left[\mathrm{NewBitGenerator}\,\mathrm{NewGenerator}\right]$

$\mathrm{seed}≔1050109365100103751001961108357097849013652340237134723870\:$

$B≔\mathrm{NewBitGenerator}\left(\mathrm{seed}\right)$

$B≔\mathbf{proc}\left(\right)z≔z\&lsqb;1\&plus;1..\&minus;1\&rsqb;\&semi;\mathbf{while}\mathrm{nops}\left(z\right)<1\mathbf{do}x≔\mathrm{irem}\left(x\&Hat;2\&comma;n\right)\&semi;z≔\left[\mathrm{op}\left(z\right)\&comma;T\&lsqb;\mathrm{irem}\left(x\&comma;1024\right)\&rsqb;\right]\mathbf{end\; do}\&semi;\mathrm{op}\left(1..1\&comma;z\right)\mathbf{end\; proc}$

$B\left(\right)$

$1$

$B\left(\right)$

$0$

$\mathrm{seq}\left(B\left(\right)\&comma;i=1..10\right)$

$0,0,1,0,1,1,0,1,1,1$

$B≔\mathrm{NewBitGenerator}\left(\mathrm{seed}\&comma;\mathrm{numbits}=10\&comma;\mathrm{primes}=1024\right)\&colon;$

$B\left(\right)$

$0,1,1,0,0,0,1,1,0,0$

$B\left(\right)$

$0,0,0,1,0,1,0,1,1,1$

$B≔\mathrm{NewBitGenerator}\left(\mathrm{seed}\&comma;\mathrm{numbits}=32\&comma;\mathrm{output}=\mathrm{integer}\right)\&colon;$

$B\left(\right)$

$3026219245$

$\mathrm{seq}\left(B\left(\right)\&comma;i=1..5\right)$

$780542672,3215656114,4009826535,4022932036,3025225627$