C/C++の宿題を片付けます 86代目
■ このスレッドは過去ログ倉庫に格納されています
0471デフォルトの名無しさん
2007/05/02(水) 20:22:14#include <stdio.h>
#include <math.h>
#define N 200000
#define LCGs(a, c, x) ((a*x)%c) /* 式は適当に変えて */
int main()
{
double p[10], t, sigma;
int a, c, x[N], i, k;
x[0] = 7, a = 1009, c = N; /* 定数と初期値 */
for(i = 2; i <= N; i++)
x[i - 1] = LCGs(a, c, x[i - 2]); /* 乱数を格納 */
for(k = 1; k <= 10; k++){
p[k - 1] = 0;
sigma = 0;
for(i = 1; i <= N - k; i++)
sigma += (double)x[i - 1]/c*x[i - 1 + k]/c;
p[k - 1] = 12.0/N*sigma - 3;
}
for(k = 1; k <= 10; k++){
t = sqrt(N - 2)*p[k - 1]/sqrt(1 - pow(p[k - 1], 2));
printf("P%d = %f\t\tt = %f\t\t", k, p[k - 1], t);
if(abs(t) < 2.326) /* 勝手に絶対値と解釈 */
printf("受容");
printf("\n");
}
return 0;
}
■ このスレッドは過去ログ倉庫に格納されています