222 lines
4.4 KiB
C#
222 lines
4.4 KiB
C#
using System;
|
|
|
|
namespace _2DGAMELIB;
|
|
|
|
public class MyRandom : Random
|
|
{
|
|
protected uint x;
|
|
|
|
protected uint y;
|
|
|
|
protected uint z;
|
|
|
|
protected uint w;
|
|
|
|
private const double d_4294967296 = 2.3283064365386963E-10;
|
|
|
|
private const double d_9007199254740992 = 1.1102230246251565E-16;
|
|
|
|
protected uint InitMtSub(uint s, uint i)
|
|
{
|
|
return 1812433253 * (s ^ (s >> 30)) + i + 1;
|
|
}
|
|
|
|
public void Initialize(uint s)
|
|
{
|
|
x = InitMtSub(s, 0u);
|
|
y = InitMtSub(x, 1u);
|
|
z = InitMtSub(y, 2u);
|
|
w = InitMtSub(z, 3u);
|
|
}
|
|
|
|
public MyRandom(uint s)
|
|
{
|
|
Initialize(s);
|
|
}
|
|
|
|
public MyRandom()
|
|
: this(4357u)
|
|
{
|
|
}
|
|
|
|
public void Initialize(uint[] key)
|
|
{
|
|
uint num = (uint)key.Length;
|
|
uint[] array = new uint[4];
|
|
Initialize(1u);
|
|
array[0] = x;
|
|
array[1] = y;
|
|
array[2] = z;
|
|
array[3] = w;
|
|
uint s = w;
|
|
uint num2;
|
|
for (num2 = 0u; num2 < num; num2++)
|
|
{
|
|
array[num2 & 3] ^= (s = InitMtSub(s, key[num2] + num2));
|
|
}
|
|
uint num3 = 0u;
|
|
while (num3 < 3)
|
|
{
|
|
array[num2 & 3] ^= (s = InitMtSub(s, num3));
|
|
num3++;
|
|
num2++;
|
|
}
|
|
x = array[0];
|
|
y = array[1];
|
|
z = array[2];
|
|
w = array[3];
|
|
if (x == 0 && y == 0 && z == 0 && w == 0)
|
|
{
|
|
x = 1u;
|
|
}
|
|
}
|
|
|
|
public MyRandom(uint[] key)
|
|
{
|
|
Initialize(key);
|
|
}
|
|
|
|
public override int Next()
|
|
{
|
|
uint num = x ^ (x << 11);
|
|
x = y;
|
|
y = z;
|
|
z = w;
|
|
return (int)(0.5 * (double)(w = w ^ (w >> 19) ^ num ^ (num >> 8)));
|
|
}
|
|
|
|
public override int Next(int les)
|
|
{
|
|
uint num = x ^ (x << 11);
|
|
x = y;
|
|
y = z;
|
|
z = w;
|
|
return (int)((double)les * 2.3283064365386963E-10 * (double)(w = w ^ (w >> 19) ^ num ^ (num >> 8)));
|
|
}
|
|
|
|
public void Next(int les, out int o)
|
|
{
|
|
uint num = x ^ (x << 11);
|
|
x = y;
|
|
y = z;
|
|
z = w;
|
|
o = (int)((double)les * 2.3283064365386963E-10 * (double)(w = w ^ (w >> 19) ^ num ^ (num >> 8)));
|
|
}
|
|
|
|
public override int Next(int sta, int les)
|
|
{
|
|
uint num = x ^ (x << 11);
|
|
x = y;
|
|
y = z;
|
|
z = w;
|
|
return (int)((double)(les - sta) * 2.3283064365386963E-10 * (double)(w = w ^ (w >> 19) ^ num ^ (num >> 8))) + sta;
|
|
}
|
|
|
|
public int NextM(int max)
|
|
{
|
|
uint num = x ^ (x << 11);
|
|
x = y;
|
|
y = z;
|
|
z = w;
|
|
return (int)((double)((long)max + 1L) * 2.3283064365386963E-10 * (double)(w = w ^ (w >> 19) ^ num ^ (num >> 8)));
|
|
}
|
|
|
|
public int NextM(int min, int max)
|
|
{
|
|
uint num = x ^ (x << 11);
|
|
x = y;
|
|
y = z;
|
|
z = w;
|
|
return (int)((double)((long)(max - min) + 1L) * 2.3283064365386963E-10 * (double)(w = w ^ (w >> 19) ^ num ^ (num >> 8))) + min;
|
|
}
|
|
|
|
public override double NextDouble()
|
|
{
|
|
uint num = x ^ (x << 11);
|
|
x = y;
|
|
y = z;
|
|
z = w;
|
|
uint num2 = (w = w ^ (w >> 19) ^ num ^ (num >> 8)) >> 11;
|
|
num = x ^ (x << 11);
|
|
x = y;
|
|
y = z;
|
|
z = w;
|
|
return 1.1102230246251565E-16 * ((double)(w = w ^ (w >> 19) ^ num ^ (num >> 8)) * 2097152.0 + (double)num2);
|
|
}
|
|
|
|
public double NextDouble(double les)
|
|
{
|
|
uint num = x ^ (x << 11);
|
|
x = y;
|
|
y = z;
|
|
z = w;
|
|
uint num2 = (w = w ^ (w >> 19) ^ num ^ (num >> 8)) >> 11;
|
|
num = x ^ (x << 11);
|
|
x = y;
|
|
y = z;
|
|
z = w;
|
|
return les * 1.1102230246251565E-16 * ((double)(w = w ^ (w >> 19) ^ num ^ (num >> 8)) * 2097152.0 + (double)num2);
|
|
}
|
|
|
|
public double NextDouble(double sta, double les)
|
|
{
|
|
uint num = x ^ (x << 11);
|
|
x = y;
|
|
y = z;
|
|
z = w;
|
|
uint num2 = (w = w ^ (w >> 19) ^ num ^ (num >> 8)) >> 11;
|
|
num = x ^ (x << 11);
|
|
x = y;
|
|
y = z;
|
|
z = w;
|
|
return (les - sta) * 1.1102230246251565E-16 * ((double)(w = w ^ (w >> 19) ^ num ^ (num >> 8)) * 2097152.0 + (double)num2) + sta;
|
|
}
|
|
|
|
public bool NextBool()
|
|
{
|
|
uint num = x ^ (x << 11);
|
|
x = y;
|
|
y = z;
|
|
z = w;
|
|
return (w = w ^ (w >> 19) ^ num ^ (num >> 8)) % 2 == 0;
|
|
}
|
|
|
|
public int NextSign()
|
|
{
|
|
uint num = x ^ (x << 11);
|
|
x = y;
|
|
y = z;
|
|
z = w;
|
|
if ((w = w ^ (w >> 19) ^ num ^ (num >> 8)) % 2 != 0)
|
|
{
|
|
return -1;
|
|
}
|
|
return 1;
|
|
}
|
|
|
|
public double NextNorCos(double mu = 0.0, double sigma = 1.0)
|
|
{
|
|
uint num = x ^ (x << 11);
|
|
x = y;
|
|
y = z;
|
|
z = w;
|
|
uint num2 = (w = w ^ (w >> 19) ^ num ^ (num >> 8)) >> 11;
|
|
num = x ^ (x << 11);
|
|
x = y;
|
|
y = z;
|
|
z = w;
|
|
double d = 1.1102230246251565E-16 * ((double)(w = w ^ (w >> 19) ^ num ^ (num >> 8)) * 2097152.0 + (double)num2);
|
|
num = x ^ (x << 11);
|
|
x = y;
|
|
y = z;
|
|
z = w;
|
|
num2 = (w = w ^ (w >> 19) ^ num ^ (num >> 8)) >> 11;
|
|
num = x ^ (x << 11);
|
|
x = y;
|
|
y = z;
|
|
z = w;
|
|
double num3 = 1.1102230246251565E-16 * ((double)(w = w ^ (w >> 19) ^ num ^ (num >> 8)) * 2097152.0 + (double)num2);
|
|
return System.Math.Sqrt(-2.0 * System.Math.Log(d)) * System.Math.Cos(System.Math.PI * 2.0 * num3) * sigma + mu;
|
|
}
|
|
}
|