梅森旋转产生随机数c语言实现,梅森旋转法产生随机数
直接贴代码了:
#include
#include
typedef unsigned long uint32;
#define N (624) // length of state vector
#define M (397)
#define K (0x9908B0DFU) // a magic constant
#define hiBit(u) ((u) & 0x80000000U) // mask highest bit of u
#define loBit(u) ((u) & 0x00000001U) // mask lowest bit of u
#define loBits(u) ((u) & 0x7FFFFFFFU) // mask all but highest bit of u
#define mixBits(u, v) (hiBit(u) | loBits(v)) // move high bit of u to low bits of v
void seedMT(uint32 seed);
uint32 reloadMT(void);
uint32 randomMT(void);
static uint32 state[N+1];
static uint32 *next;
static int left = -1;
void seedMT(uint32 seed)
{
register uint32 x = (seed | 1U) & 0xFFFFFFFFU, *s = state;
register int j;
for(left = 0, *s++=x, j=N; --j;
*s++ = (x*69069U) & 0xFFFFFFFFFF);
}
uint32 reloadMT(void)
{
register uint32 *p0=state, *p2=state+2, *pM=state+M, s0, s1;
register int j;
if(left < -1)
seedMT(4357U);
left=N-1, next=state+1;
for(s0=state[0], s1=state[1], j=N-M+1; --j; s0=s1, s1=*p2++)
*p0++ = *pM++ ^ (mixBits(s0, s1) >> 1) ^ (loBit(s1) ? K : 0U);
for(pM=state, j=M; --j; s0=s1, s1=*p2++)
*p0++ = *pM++ ^ (mixBits(s0, s1) >> 1) ^ (loBit(s1) ? K : 0U);
s1 = state[0], *p0 = *pM ^ (mixBits(s0, s1) >> 1) ^ (loBit(s1) ? K : 0U);
s1 ^= (s1 >> 11);
s1 ^= (s1 << 7) & 0x9D2C5680U;
s1 ^= (s1 << 15) & 0xEFC60000U;
return (s1 ^ (s1 >> 18));
}
uint32 randomMT(void)
{
uint32 y;
//printf("... left = %d\n", left);
if(--left < 0)
return (reloadMT());
//printf("... *next = %d\n", *next);
y = *next++;
y ^= (y >> 11);
y ^= (y << 7) & 0x9D2C5680U;
y ^= (y << 15) & 0xEFC60000U;
y ^= (y >> 18);
return (y);
}
// test only
int main() {
int seed;
scanf("%d", &seed);
seedMT(seed);
int cnt = 1000000, a = 0, b = 0;
while(cnt--) {
if(randomMT()%2==0) {
a++;
}
else b++;
reloadMT();
}
printf("%d %d\n", a, b);
return 0;
}
梅森旋转产生随机数c语言实现,梅森旋转法产生随机数相关推荐
- JAVA梅森旋转随机算法,C语言实现梅森旋转算法
原理: 梅森旋转算法是一种伪随机数生成算法 1.定义相关宏和结构体 #include #include #include #include typedef unsigned int uint; #de ...
- 线性同余法产生随机数C语言,线性同余生成随机数的一点思考
今天下午 pk 和我讨论了一个问题,他看到在另一个项目组的 lua 代码里有一段使用线性同余产生随机数的代码,但是那个项目组的同事告诉他这个函数生成的随机数是分布不均的.于是他想到了我前两天给他讲的关 ...
- 梅森旋转算法原理c语言,梅森旋转素数算法(MT199937)c语言代码
满意答案 sddzycz 2016.06.23 采纳率:42% 等级:12 已帮助:8913人 #include /*Periodparameters*/ #define N 624 #defi ...
- c语言产生服从正态分布的随机数,C语言产生满足正态分布的随机数
C语言中可以通过rand函数生成满足均匀分布的随机数,但是生成满足正太分布的随机数就没有那么简单了,下面对常用的几种方法进行总结并用C++编程实现. 方法一:由均匀分布的随机数来产生 一个简单可行的并 ...
- JAVA梅森旋转随机算法_梅森旋转算法
概念 梅森旋转算法(Mersenne twister),可以快速产生高质量的伪随机数,修正了古典随机数发生算法的很多缺陷. 常见的两种为基于32位的 MT19937和基于64位的 MT19937-64 ...
- JAVA梅森旋转随机算法_伪随机数生成算法-梅森旋转(Mersenne Twister/MT)
今天主要是来研究梅森旋转算法,它是用来产生伪随机数的,实际上产生伪随机数的方法有很多种,比如线性同余法, 平方取中法等等.但是这些方法产生的随机数质量往往不是很高,而今天介绍的梅森旋转算法可以产生高质 ...
- c语言中什么情况下会出现随机数,C语言中如何产生随机数
今天看到一段小程序 ,里面用到随机数.才发现在C语言中产生随机数不像matlab中那么简单. C中也有rand()函数,但是rand()函数产生的数不是真正意义上的随机数,是一个伪随机数,是根据一个数 ...
- 梅森旋转产生随机数c语言实现,C++生成随机数的实现代码
C风格 C为随机数提供的工具是rand.srand和RAND_MAX,定义在中. srand为rand设置种子,如果不设置,相当于调用过srand(1).rand产生伪随机数,其范围为0到RAND_M ...
- js版梅森旋转生成随机数
用js实现, 梅森旋转生成随机数, 来代替Math.random(); 该文章参考自 梅森旋转算法--伪随机数(加密.身份信息ID号)_Touch_Dream的博客-CSDN博客 <!DOCTY ...
最新文章
- Java---实现运行任意目录下class中加了@MyTest的空参方法(实现图形界面)
- pip 安装指定版本
- LeetCode35.搜索插入位置 JavaScript
- wxWidgets:wxPython 概述
- centos6.2 64安装mysql_centos6.2 64bit安装mysql+php
- php mysql 线程安全_PHP 线程安全与非线程安全版本的区别深入解析
- 论文浅尝 | NumNet: 一种带有数学推理的机器阅读理解模型
- 云云协同解决方案全景图发布 华为云助力科技企业云上创新
- insert into select 优化_数据库优化总结
- [转载]java对cookie的操作
- snap7 python连接_python snap7 windows-找不到snap7库
- VBA InStr 函数
- 巧解快速排序时间复杂度
- Python爬虫入门教程:豆瓣Top电影爬取
- 大数据与个人隐私的平衡:懂你,但不认识你
- 元数据与数据治理|apache-atlas2.0.0 源码编译和安装部署(CHD版本)
- 电磁散射特性与SAR影像之间的关联,典型地物的散射特性
- 半色调技术简介(一)
- 【技术知识】SVAC 2.0安全技术浅析
- Linux定时器alarm()用法讲解