生成随机数.cpp

1.书中所有的算法均用C语言描述

#include <iostream>
#include <time.h>
using namespace std;

void rnds(double p[],int n)
{ //生成n个[0 1]间的均匀随机数。
	int i,m;
	unsigned long randSeed;
	randSeed=time(0);//取系统时间为当前种子
    double x0=randSeed+0.0;
    double  *r=&x0; 
	double s,u,v;
	s=65536.0; u=25741.0; v=13849.0;
	for (i=0; i<=n-1; i++)
	{	
		*r=u*(*r)+v;
		m=(int)(*r/s);
		*r=*r-m*s; 
		p[i]=*r/s;
	}
    return;
}
void main()
{ 
	int i,j;
	double p[100];
    rnds(p,100);

	for (i=0; i<=19; i++)
	{ 
		for (j=0; j<=4; j++)
		cout<<p[5*i+j]<<" ";
		cout<<endl;
	}

}