Lo que sigue está traducido de un post de G. Marsaglia en sci.crypt. Proporcionan varios generadores de números aleatorios que ocupan poco espacio y parecen tener una buena calidad.

    Comienzo del código C:

    Fin del código C.

    Basta colocar las 17 líneas de código al comienzo de tu programa, inicializar la tabla y tener varios RNG muy prometedores.

    Posiblemente quieras cambiar los nombres de las variables z, w, x, y, t, c por algunos más largos para evitar el conflicto con otras con el mismo nombre en tu código.

    Los generadores KISS, MWC, LFIB4, SWB, SHR3 y CONG proporcionan un entero aleatorio de 32 bits. UNI produce una distribución uniforme en (0,1), VNI hace lo mismo en (-1,1).

    Si quieres un generador de periodo astronómicamente grande KISS+SWB tiene un periodo mayor de 2⁷⁷⁰⁰!!

    El generador KISS (keep it simple stupid) combina los dos generadores de multiplicacion con acarreo en MWC con los tres registros de desplazamiento en SHR3 y el generador congruencial CONG usando la suma y XOR. Período aproximado 2¹²³.

    El generador MWC concatena dos generadores de 16 bits del tipo multiplicacion con acarreo, x(n)=36969x(n-1)+carry, y(n)=18000y(n-1)+carry mod 2¹⁶. Tiene un periodo de 2⁶⁰ aproximadamente y al parecer pasa todos los tests de aleatoriedad. Interesante por sí solo aunque sea parte de KISS.

    SHR3 es un generador con tres registros de desplazamiento con periodo 2^32-1 . Usa y(n)=y(n-1)(I+L¹⁷)(I+R¹³)(I+L⁵) con las y's vistas como vectores binarios, L la matriz binaria de orden 32x32 que desplaza un vector a la izquierda 1 y R su transpuesta.
SHR3 pasa todos los tests excepto el de rango binario pues 32 valores sucesivos vistos como vectores binarios deben ser linealmente independientes, mientras 32 enteros aleatorios de 32 bits , serán independientes solo el 29% de las veces.

    CONG es un generador congruencial con el muy popular multiplicador 69069: x(n)=69069x(n-1)+1234567. Tiene periodo 2³². La mitad superior de sus 32 bits pasan bien los test de aleatoriedad pero los 32 inferiores parecen ser demasiado regulares.

    LFIB4 es una extension de la clase anteriormente definida como generadores retrasados (lagged) de Fibonacci: x(n)=x(n-r) op x(n-s) donde las x's pertenecen a un conjunto finito donde hay definida una operación binaria como + o - en los enteros módulo 2³², o como * en los enteros impares citados o la operacion XOR en los vectores binarios. Excepto los que usan multiplicaciones este tipo de generadores fallan varios tests de aleatoriedad excepto si los retrasos (lags) son muy largos. Marsaglia ha desarrollado el generador con 4 retrasos LFIB4, x(n)=x(n-256)+x(n-179)+x(n-119)+x(n-55) mod2³².
Tiene de periodo 2³¹*(2^256-1) o aproximadamente 2²⁸⁷ y parece pasar todos los tests. Combinado con KISS produce un periodo de 2⁴¹⁰, basta usar KISS+LFIB4 en cualquier expresión C.

    SWB es un generador de resta con acarreo desarrollado para mostrar que un método simple puede producir periodos muy largos. x(n)=x(n-222)-x(n-237)-acarreo mod 2³². El acarreo es 0 excepto cuando calcular x(n-1) produjo un overflow en cuyo caso estará a 1. Este generador tiene un periodo MUY largo 2⁷⁰⁹⁸(2^480-1), aproximadamente 2⁷⁵⁷⁸ Parece pasar todos los test de aleatoriedad pero Marsaglia recomienda combinarlo con KISS, MWC, SHR3 o CONG.

    Finalmente, UNI y VNI se emplean cuando se necesitan números reales en el intervalo (0,1) o el (-1,1) respectivamente.

    Los procedimientos MWC, SHR3, CONG, KISS, LFIB4, SWB, UNI y VNI permiten insertarse directamente en cualquier expresión C evitando la carga de tiempo y espacio que lleva una función ordinaria. Las variables globales z, w, jsr y jcong deben tener valores asignados previamente.

    Si se usan LFIB4 o SWB debe haberse inicializado anteriormente la tabla t[256].