下面的平方取中生成随机数有没有问题

[采石工]

#include <cmath>
#include <iostream>
using namespace std;

int pow10(int n)
{
	int ret=1;
	for(int i=0;i<n;++i)
	{
		ret*=10;
	}
	return ret;
}

int getplaces(int a)
{
	int n=1;
	while(a/10)
	{
		++n;
		a/=10;
	}
	return n;
}
unsigned int k;
unsigned int next=1;
void srandom(unsigned int seed)
{
	next=seed;
	k=getplaces(next);
}

double random()
{
	double r=(double)next/pow10(k);
	next=((next*next)/pow10(k/2))%pow10(k);
	//next=fmod((double)(next*next)/pow10(k/2),pow10(k));
	return r;
}

void test(const int n)
{
	int a[10]={0};
	for(int i=0;i<n;++i)
	{
		double r=random();
		int j=floor(r*10);
		a[j]++;
	}
	for(int i=0;i<10;++i)
	{
		cout<<a[i]<<endl;
	}
}

int main()
{
	srandom(5466);
	test(5000);
	return 0;
}

测试结果显示，生成的随机值并非近似于均匀分布，不知代码有问题，还是确实为非均匀分布

[5t4rk]

本来就是伪随机数 sleep一下试试。

[采石工]

引用 1 楼 bjtbjt 的回复:

本来就是伪随机数 sleep一下试试。

你好，我知道是随机。我仿照matlab的代码写了上述代码，前者测试结果是：随机序列近似于均匀分布。而上面代码得到的结果却明显不是。不知何故