पार्क-मिलर रैंडम नंबर जनरेशन एल्गोरिथम को लागू करने के लिए C++ प्रोग्राम

पार्क-मिलर रैंडम नंबर जेनरेशन एल्गोरिथम रैंडम नंबर जेनरेट करने का एक और तरीका है।

इस प्रकार के एक यादृच्छिक संख्या जनरेटर (RNG) का एक सामान्य सूत्र है:X_{k+1} =g X(k) mod n

जहां मापांक n एक अभाज्य संख्या या अभाज्य संख्या का घात है, गुणक g उच्च गुणक क्रम मॉड्यूल n का एक तत्व है, और बीज X0 n से सहअभाज्य है।

एल्गोरिदम

Begin
   Declare variables n, a, b, c and seed
   Read variables n, a, b, c and seed
   Uniform()
   Declare variable hi, lo, t
   hi=seed divided by b
   lo = seed - b * hi
   t = a * lo - c * hi
   if (t > 0)
      seed = t;
   else
      seed = t + n;
      return seed;
   Done
   For i =0 to n
      Call the function random
   Done
End

उदाहरण कोड

#include <iostream>
using namespace std;
const long n = 2145678965L;
const long a = 763214L;
const long b = 88844L;
const long c = 7766L; i
static long seed = 12345678L;
double uniform() {
   long hi = seed / b;
   long lo = seed - b * hi;
   long t = a * lo - c * hi;
   if (t > 0)
      seed = t;
   else
      seed = t + n;
   return seed;
}
int main(int argc, char **argv) {
   double A[10];
   for (int i = 0; i < 10; i++)
      A[i] = uniform();
      cout << "Random numbers are:\n";
   for (int i = 0; i < 10; i++)
      cout << A[i] << endl;
}

आउटपुट

Random numbers are:
6.50293e+10
4.27187e+10
2.1539e+10
4.62058e+10
1.70792e+10
8.24569e+09
5.93381e+10
3.63839e+10
4.81931e+10
8.91007e+09