C语言线性同余法产生随机数

线性同余法产生随机数序列的公式为：
    x[k+1] = (a*x[k] + c) % m

参数 a, c, m 都是由用户来设定的，并和一个种子数（例如 x[0]）一起作为命令行参数传递给生成的程序。

一个简单的例子：a=7, c=1, m=13, and seed=5

一个复杂的例子：a=69069, c=0, m=2^32=4294967296, seed=31

下面的代码将输出一个随机数序列（最多有 m-1 个不同的值），然后继续循环。

#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) {  // 输出 8 个数字以后换行
            printf("\n");
            Sleep(1);  // 暂停 1 秒
        }
    }
    printf("\n");
    return 0;
}

构建成功后，需要到命令行（cmd.exe）运行该程序才能看到效果。这里涉及到 main() 函数的传参问题，请查看：C语言main()函数详解

比如，我的项目放在 E:\cDemo\ 目录下，在命令行中的输入截图如下：


注意：输入的参数个数（cDemo.exe也算一个参数）必须等于 5 程序才能正确执行。

上面的截图中，我给 m 传的值为 56，产生了 55( m-1 ) 个随机数。