Linear Congruential Generator

A linear congruential generator (LCG), is an algorithm that produces a sequence of pseudo-random numbers calculated with a piecewise linear equation.

Usamos uma relação de recorrência linear da seguinte forma:

[math]X_i = m.X_{i-1}+p[/math]

onde m e p são dois números inteiros adequadamente escolhidos.

E as variáveis aleatórias uniformes são obtidas após o dimensionamento,

[math]U_i = \frac{X_i}{p}[/math]

Exemplo em Python

commoldura|direita|1 milhão de sorteios usando LCG

def pseudo_uniform_good(mult=16807,mod=(2**31)-1,seed=123456789,size=1):
    """
    A reasoanbly good pseudo random generator
    """
    
    U = np.zeros(size)
    x = (seed * mult + 1) % mod
    U[0] = x / mod

    for i in range(1, size):
        x = (x * mult + 1) % mod
        U[i] = x / mod

    return U