我正在编写一个函数,该函数接受输入n并创建大小为(2n-1)^ 2的一维数组以模拟正方形。即对于n = 1的输入,只有一个点;对于n = 2的输入,它看起来像
0 1 2
3 4 5
6 7 8
对于n = 3,它看起来像
0 1 2 3 4
5 6 7 8 9
10 11 12 13 14
15 16 17 18 19
20 21 22 23 24
每个数字都是一个点。
当在边缘上检测到当前位置并且该点试图移出正方形网格时,该函数终止。
这样做的目的是模拟从n = 2 ^ 0到n = 2 ^ 8的不同大小的正方形访问了多少个点,并返回访问了多少点的总和点在正方形中。
该函数生成一个随机数并将其模数与4进行比较,如果返回0,则将位置向上移动1,如果返回1,则将位置向右移动,向下2,向左3
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
double two_d_random (int n) {
int tot_points = (2 * n - 1)*(2 * n - 1);
int value = (n*n)+((n-1)*(n-1))-1; //center
int length = 2 * n - 1; //length of side
int *array = (int *)malloc (sizeof (int) * tot_points);
int count = 0;
array[value] = 1;
while (1 == 1) {
int r = rand () % 4;
array[value] = 1;
if (r == 0) {//UP
if ((value >= 0) && (value < length)) {
goto a;
}
else {
array[value] = 1;
value -= length;
}
}
else if (r == 1) {//RIGHT
if ((value % length) == (2*n-2)){
goto a;
}
else {
array[value] = 1;
value += 1;
}
}
else if (r == 2) {//DOWN
if ((value < tot_points) && (value >= (tot_points - length))) {
goto a;
}
else {
array[value] = 1;
value += length;
}
}
else if (r == 3) {//LEFT
if (value % length == 0) {
goto a;
}
else {
array[value] = 1;
value -= 1;
}
}
}
a:
for (int i = 0; i < tot_points; i++) {
if (array[i] == 1) {
count += 1;
}
}
free (array);
return 1.0 * count / tot_points;
}
int main ()
{
int trials = 1000;
srand (12345);
for (int n = 1; n <= 256; n *= 2)
{
double sum = 0.;
for (int i = 0; i < trials; i++)
{
double p = two_d_random(n);
sum += p;
}
printf ("%d %.3lf\n", n, sum / trials);
}
return 0;
}
我当前的问题是,当我在机器上运行它时,我得到了一系列我不期望的值:
但是,当一个同事在他们的机器上运行它时,他们得到了我所期望的以下内容:
我意识到这是一个很大的问题。我也意识到我不应该使用goto。我已经花了很多时间,但是我不知道如何解决这个问题。任何帮助将不胜感激。
答案 0 :(得分:1)
您需要初始化数组。调用malloc()
时,它只返回一块未初始化的内存。对其进行初始化,或者使用calloc()
来获取预先清零的内存。
double two_d_random( int n )
{
int tot_points = ( 2 * n - 1 ) * ( 2 * n - 1 );
int value = ( n * n ) + ( ( n - 1 ) * ( n - 1 ) ) - 1; //center
int length = 2 * n - 1; //length of side
int *array = (int *) malloc( sizeof( int ) * tot_points );
int count = 0;
// Initialise the array to zero
for ( int i=0; i<tot_points; i++ )
{
array[i] = 0;
}
array[value] = 1;
while ( 1 == 1 )
{
int r = rand() % 4;
array[value] = 1;
if ( r == 0 )
通过此修改,我得到的结果与您所报告的结果 类似>
1 1.000
2 0.367
4 0.221
8 0.154
16 0.122
32 0.101
64 0.085
128 0.077
256 0.071