RandomTools[BlumBlumShub] - Maple Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Programming : Random Objects : RandomTools package : BlumBlumShub Subpackage : RandomTools/BlumBlumShub/NewGenerator

RandomTools[BlumBlumShub]

  

NewGenerator

  

Blum, Blum and Shub Pseudo Random Bit Generator

 

Calling Sequence

Parameters

Description

Examples

Calling Sequence

NewBitGenerator( opt1, opt2 )

Parameters

opt1, opt2

-

(optional) argument of the form option=value where option is one of range or seed

Description

• 

The NewGenerator command outputs a Maple procedure, a pseudo-random number generator, which when called outputs a pseudo-random number. 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:

xk+1=modxk2,n

zk+1=modxk+1,2m.

  

where

– 

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

– 

m=floorlog[2]log[2]n, and

– 

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 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) gcdx,n=1 and (ii) x=y2modn for some integer y.

  

Definition: An integer z in Z mod n where n=pq is a pseudo-square if (i) gcdz,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=2s+1 and q=2t+1 where s and t are prime and 2 is a primitive element in both Z mod s and in 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 s1 or t1 or lcms1,t1, 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 following optional arguments are supported. They are input as equations in any order.

  

range=integer..integer or integer

  

The range option specifies the range from which the integer is chosen.  If a range is given, the returned procedure will generate numbers in this range. If an integer is given, the returned procedure will return integers in the range 0 to value-1. By default the range is 0 to 10^12-1.

  

seed=integer

  

The option seed determines the seed for the generator. The value used for x0 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 value of the seed option.  If a seed is not given, a seed is generated by calling RandomTools[MersenneTwister][GenerateInteger] with the argument range = 10^99..10^100.

• 

The RandomTools[BlumBlumShub] module also exports the NewBitGenerator function. This function provides the user with greater control over the properties of the returned procedure.  For cryptographic applications, the NewBitGenerator function should be used.

Examples

withRandomTools[BlumBlumShub]

NewBitGenerator,NewGenerator

(1)

s1050109365100103751001961108357097849013652340237134723870:

BNewGeneratorseed=2

B:=proclocaln&comma;p&semi;p:=proclocali&semi;forifromcby10to40dox:=iremx&Hat;2&comma;n&semi;z:=1024&ast;z&plus;iremx&comma;1024end do&semi;c:=i40&semi;iquoz&comma;2&Hat;c&comma;&apos;z&apos;end proc&semi;n:=p&semi;while1000000000000<=ndon:=pend do&semi;nend proc

(2)

B

493660033144

(3)

B

259218944679

(4)

seqB&comma;i&equals;1..10

17627898170&comma;966628291159&comma;628124827667&comma;692163973638&comma;365738744555&comma;283269255713&comma;255572006637&comma;712645761945&comma;712457622044&comma;426819767701

(5)

BNewGeneratorseed&equals;s&comma;range&equals;1..1000&colon;

B

722

(6)

B

520

(7)

BNewGeneratorseed&equals;s&comma;range&equals;1020..1030&colon;

B

893182270928831182291970575150

(8)

seqB&comma;i&equals;1..5

571028873893626957286791880500&comma;149691015279629477184761484186&comma;893310462378765125162315235384&comma;749992797388668654605785130192&comma;462774729207078150350323115146

(9)

See Also

rand

RandomTools

RandomTools[BlumBlumShub]

RandomTools[BlumBlumShub][NewBitGenerator]

RandomTools[Generate]

 


Download Help Document

Was this information helpful?



Please add your Comment (Optional)
E-mail Address (Optional)
What is ? This question helps us to combat spam