using System;

​

// ============================================================================================================

// Romu - Fast Nonlinear Pseudo-Random Generator

//        Based on https://arxiv.org/pdf/2002.11331.pdf. Discussion here: http://www.romu-random.org/

//        The RomuTrio32 version is uses 32-bit calculations on 96-bits of state and able

//        to generate 2^53 numbers without patterns or flaws showing up

// ============================================================================================================

​

public static class RomuRND32

{

    private static UInt32 xState, yState, zState; // set to nonzero seed

    public static void init( UInt32 uX, UInt32 uY, UInt32 uZ ) {

        xState = uX;

        yState = uY;

        zState = uZ;

        // Run three rounds so that initial seeds are fully mixed into the 96-bit state.

        // Without these rounds, the first 1 to 3 generated numbers will be 

        // the same even when 1st 32-bits of the seed are the same.

        nextRandom();

        nextRandom();

        nextRandom();

    }

​

    public static UInt32 nextRandom () {

        UInt32 xp = xState, yp = yState, zp = zState;

        // Console.WriteLine("state: {0:X8}-{1:X8}-{2:X8}", xState, yState, zState);

        xState = 3323815723u * zp;

        yState = yp - xp; yState = RotateLeft(yState,6);

        zState = zp - yp; zState = RotateLeft(zState,22);

        return xp;

    }

    public static uint RotateLeft(this uint value, int count)

    {

        return (value << count) | (value >> (32 - count));

    }

}

​

public static class StatsCollector

{

    //                             "00112233445566778899AABBCCDDEEFF"

    private static string sGRAPH = "..--::--++xx**!!||()[]%%&&@@$$##";

    public static UInt32[] uBits = new UInt32[32];

    public static UInt32[] uBytes = new UInt32[256];

    public static void collectValue( UInt32 uValue ) {

        UInt32 uBit = 1;

        for ( int iBit = 0; iBit < 32; ++iBit ) {

            if ( (uValue & uBit) != 0 ) ++uBits[iBit];

            uBit <<= 1;

        }

        /*

        string s = "";

        for ( int iBit = 0; iBit < 32; ++iBit ) {

            if ( iBit > 0 ) s += ".";

            s += String.Format("{0:00}",uBits[iBit]);

        }

        Console.WriteLine("Add {0:X8}: Bit stats = {1}", uValue, s );

        */

        UInt32 uByte;

        uByte = uValue & 0x000000FF; ++uBytes[uByte]; uValue >>= 8;

        uByte = uValue & 0x000000FF; ++uBytes[uByte]; uValue >>= 8;

        uByte = uValue & 0x000000FF; ++uBytes[uByte]; uValue >>= 8;

        uByte = uValue & 0x000000FF; ++uBytes[uByte];

    }

    public static void printStats( int nValues) {

        for ( int iBit = 0; iBit < uBits.Length; ++iBit ) {

            // Console.WriteLine( "Debug uBits[{0}] = {1}", iBit, uBits[iBit] );

            printGraphLine( String.Format("bit {0:00}", iBit), uBits[iBit], nValues, 80 );

        }

    }

    public static void printGraphLine( string sLabel, uint nValue, int nMaxValue, int nWidth ) {

        string sLine = "";

        sLine += String.Format( "{0}: {1:000} : ", sLabel, nValue );

        decimal dValue = nWidth * ((decimal)nValue / nMaxValue);

        int nPips = Decimal.ToInt32(dValue);

        // sLine += String.Format(" (n/max = {0}, pips = {1})", (decimal) nValue / nMaxValue, nPips);

        for ( int i = 0; i < nPips; ++i ) {

            sLine += "#";

        }

        Console.WriteLine( sLine );

    }

}

​

​

public class Program

{

    public static void Main()

    {

        Random rndSeed = new Random();

        UInt32 uSeed = (UInt32) rndSeed.Next();

        Console.WriteLine("seed = {0:X8}",uSeed);

        RomuRND32.init(uSeed,58146387,267884012);

        UInt32 r;

        int nVALUES = 100;

        string sOut = "";

        for ( int i = 0; i < nVALUES; ++i ) {

            r = RomuRND32.nextRandom();

            if ( sOut.Length > 0 ) sOut += ".";

            sOut += String.Format("{0:X8}", r);

            if ( sOut.Length >= 71 ) { Console.WriteLine( sOut ); sOut = ""; }

            StatsCollector.collectValue(r);

        }

        Console.WriteLine();

        Console.WriteLine("bit stats:");

        StatsCollector.printStats(nVALUES);

    }

}

​

​