Update 线性同余法产生随机数.c

ls-Brynn · web-flow · commit aa3fa01c762d · 2021-08-21T15:14:04.000+08:00
diff --git a/&#32463;&#20856;&#39033;&#30446;/&#32447;&#24615;&#21516;&#20313;&#27861;&#20135;&#29983;&#38543;&#26426;&#25968;.c b/&#32463;&#20856;&#39033;&#30446;/&#32447;&#24615;&#21516;&#20313;&#27861;&#20135;&#29983;&#38543;&#26426;&#25968;.c
@@ -1 +1,54 @@
-#
+//&#32447;&#24615;&#21516;&#20313;&#27861;&#20135;&#29983;&#38543;&#26426;&#25968;&#24207;&#21015;&#30340;&#20844;&#24335;&#20026;&#65306;
+//    x[k+1] = (a*x[k] + c) % m
+
+//&#21442;&#25968; a, c, m &#37117;&#26159;&#30001;&#29992;&#25143;&#26469;&#35774;&#23450;&#30340;&#65292;&#24182;&#21644;&#19968;&#20010;&#31181;&#23376;&#25968;&#65288;&#20363;&#22914; x[0]&#65289;&#19968;&#36215;&#20316;&#20026;&#21629;&#20196;&#34892;&#21442;&#25968;&#20256;&#36882;&#32473;&#29983;&#25104;&#30340;&#31243;&#24207;&#12290;
+
+//&#19968;&#20010;&#31616;&#21333;&#30340;&#20363;&#23376;&#65306;a=7, c=1, m=13, and seed=5
+
+//&#19968;&#20010;&#22797;&#26434;&#30340;&#20363;&#23376;&#65306;a=69069, c=0, m=2^32=4294967296, seed=31
+
+//&#19979;&#38754;&#30340;&#20195;&#30721;&#23558;&#36755;&#20986;&#19968;&#20010;&#38543;&#26426;&#25968;&#24207;&#21015;&#65288;&#26368;&#22810;&#26377; m-1 &#20010;&#19981;&#21516;&#30340;&#20540;&#65289;&#65292;&#28982;&#21518;&#32487;&#32493;&#24490;&#29615;&#12290;
+
+
+#include <stdio.h>
+#include <unistd.h>
+#include <stdlib.h>
+#include <Windows.h>
+
+static long seed = 13;
+static long a;
+static long c;
+static long m;
+
+void random_init(long s) {
+    if (s != 0) seed = s;
+}
+
+long random() {
+    seed = (a*seed + c)%m;
+    return seed;
+}
+
+int main(int argc, char * argv[]) {
+    long s;
+    int k;
+    if (argc != 5) {
+        printf("usage: %s a, c, m, seed\n", argv[0]);
+        return 1;
+    }
+    a = atoi(argv[1]);
+    c = atoi(argv[2]);
+    m = atoi(argv[3]);
+    s = atoi(argv[4]);
+    random_init(s);
+  
+    for (k = 0; k < m-1; k++) {
+        printf("%8ld", random());
+        if (k % 8 == 7) {  // &#36755;&#20986; 8 &#20010;&#25968;&#23383;&#20197;&#21518;&#25442;&#34892;
+            printf("\n");
+            Sleep(1);  // &#26242;&#20572; 1 &#31186;
+        }
+    }
+    printf("\n");
+    return 0;
+}
