您当前的位置:首页 > 计算机 > 编程开发 > C语言

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

时间:12-29来源:作者:点击数:

线性同余法产生随机数序列的公式为:
    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 ) 个随机数。

方便获取更多学习、工作、生活信息请关注本站微信公众号城东书院 微信服务号城东书院 微信订阅号
推荐内容
相关内容
栏目更新
栏目热门