Blum, Blum and Shub Pseudo Random Bit Generator
NewBitGenerator( seed, opt1, opt2, ... )
an integer, the seed for the generator
opt1, opt2, ...
(optional) argument of the form option=value where option is one of primes, numbits, or output
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 x1,⁢x2,⁢… from which cryptographically secure pseudo-random bits z1,⁢z2,⁢… are extracted:
n is a product of two large primes p and q,
x0 is determined from the seed s.
Each iteration of zk+1, the least significant m bits of xk+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⁡x,n=1 and (ii) x=y2modn for some integer y.
Definition: An integer z in Z mod n where n=p⁢q is a pseudo-square if (i) gcd⁡z,n=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=2⁢s+1 and q=2⁢t+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 x0 for any quadratic residue other than 1 is either s−1 or t−1 or lcm⁡s−1,t−1, all of which are large. Random primes p and q, and their product n=p⁢q, 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 x0 used must be a "random" quadratic residue which is not equal to 1 (to avoid a short period). The value used for x0 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.
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.
seed ≔ 1050109365100103751001961108357097849013652340237134723870:
B ≔ NewBitGenerator⁡seed
B ≔ procz ≔ z[1+1..−1];whilenops⁡z<1dox ≔ irem⁡x^2,n;z ≔ op⁡z,T[irem⁡x,1024]end do;op⁡1..1,zend proc
B ≔ NewBitGenerator⁡seed,numbits=10,primes=1024:
B ≔ NewBitGenerator⁡seed,numbits=32,output=integer:
Download Help Document