Last active
March 15, 2022 01:27
-
-
Save adamveld12/6c0350d1cfd2da449dc6 to your computer and use it in GitHub Desktop.
Mersenne prime random number generator in C#
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
/// <summary> | |
/// Implementation of Mersenne Twister random number generator | |
/// </summary> | |
public class MersennePrimeRandom | |
{ | |
private readonly uint[] _matrix = new uint[624]; | |
private int _index = 0; | |
public MersennePrimeRandom() : this((uint)(0xFFFFFFFF & DateTime.Now.Ticks)) { } | |
/// <summary> | |
/// Initializes a new instance of the MersennePrimeRNG with a seed | |
/// </summary> | |
/// <param name="seed"></param> | |
public MersennePrimeRandom(uint seed) | |
{ | |
_matrix[0] = seed; | |
for (int i = 1; i < _matrix.Length; i++) | |
_matrix[i] = (1812433253 * (_matrix[i - 1] ^ ((_matrix[i - 1]) >> 30) + 1)); | |
} | |
/// <summary> | |
/// Generates a new matrix table | |
/// </summary> | |
private void Generate() | |
{ | |
for (int i = 0; i < _matrix.Length; i++) | |
{ | |
uint y = (_matrix[i] >> 31) + ((_matrix[(i + 1) & 623]) << 1); | |
_matrix[i] = _matrix[(i + 397) & 623] ^ (y >> 1); | |
if (y % 2 != 0) | |
_matrix[i] = (_matrix[i] ^ (2567483615)); | |
} | |
} | |
/// <summary> | |
/// Generates and returns a random number | |
/// </summary> | |
/// <returns></returns> | |
public int Next() | |
{ | |
if (_index == 0) | |
Generate(); | |
uint y = _matrix[_index]; | |
y = y ^ (y >> 11); | |
y = (y ^ (y << 7) & (2636928640)); | |
y = (y ^ (y << 15) & (4022730752)); | |
y = (y ^ (y >> 18)); | |
_index = (_index + 1) % 623; | |
return (int)(y % int.MaxValue); | |
} | |
/// <summary> | |
/// Generates and returns a random number | |
/// </summary> | |
/// <param name="max">The highest value that can be returned</param> | |
/// <returns></returns> | |
public int Next(int max) | |
{ | |
var randomValue = Next(); | |
return randomValue % max; | |
} | |
/// <summary> | |
/// Generates and returns a random number | |
/// </summary> | |
/// <param name="min">The lowest value returned</param> | |
/// <param name="max">The highest value returned</param> | |
/// <returns></returns> | |
public int Next(int min, int max) | |
{ | |
if(min > max) | |
throw new ArgumentException("min cannot be greater than max", "min"); | |
return min + Next(min - max); | |
} | |
} |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
I'm not sure if this is correct or not, but in my testing
Next(int max)
never returnsmax
; the maximum it returns ismax - 1
. I modified line 63 to bereturn randomValue % (max + 1)
and it began returningmax
as the maximum number generated.