ID非公開

2022/5/21 10:52

22回答

n個の点(x, y)を乱数で生成させ，モンテカルロ法により円周率を近似 せよ.nの値は100，10000，1000000の場合を出力すること.また，点 の座標は構造体にし，

n個の点(x, y)を乱数で生成させ，モンテカルロ法により円周率を近似 せよ.nの値は100，10000，1000000の場合を出力すること.また，点 の座標は構造体にし， 円周率の計算は関数MonteCarlomethodを作っ て行うこと.円周率の出力は小数点5位まで表示すること. 点の座標は0~1の範囲(double)にすること. (double)rand() / RAND_MAX; で乱数を生成できる. 誰か優しい方教えてください。 全然わかんなくて困っています。 プログラムだけで大丈夫です。

C言語関連・114閲覧

ベストアンサー

0
Windcat

カテゴリマスター

2022/5/21 11:46

(1150911677さん) 例： #include <stdio.h> #include <stdlib.h> #include <time.h> #define N 10000000 typedef struct { double x; double y; } XY; double my_random( void ) { return (double)rand() / RAND_MAX; } XY *makeData(int n) { XY *p; int i; p=(XY*)malloc(sizeof(XY)*n); for(i=0; i<n; i++) { p[i].x = my_random(); p[i].y = my_random(); //printf("(%f, %f)\n") } return p; } double MonteCarlomethod(XY *d, int n) { int i, cnt=0; for(i=0; i<n; i++) if((d[i].x * d[i].x + d[i].y * d[i].y)<=1.0)cnt++; return (double)cnt / n * 4.0; } void ope(int n) { XY *dt; double ret; dt=makeData(n); ret=MonteCarlomethod(dt,n); printf("n=(%7d) pi = %.5f\n",n,ret); free(dt); } int main(void) { srand((unsigned)time(NULL)); rand(); ope(100); ope(10000); ope(1000000); return 0; } ~~~~~~~~~~~~~~~~~~~~~~ n=( 100) pi = 3.24000 n=( 10000) pi = 3.13920 n=(1000000) pi = 3.14105 n=( 100) pi = 3.60000 n=( 10000) pi = 3.11720 n=(1000000) pi = 3.13969 n=( 100) pi = 3.28000 n=( 10000) pi = 3.12760 n=(1000000) pi = 3.14072

ThanksImg質問者からのお礼コメント

ありがとうございます！！

お礼日時：5/22 19:51

その他の回答（1件）

0
匿

カテゴリマスター

2022/5/21 11:29

#include <stdio.h> #include <stdlib.h> #include <time.h> // 点 の座標は構造体 typedef struct {     double x;     double y; } POINT; double MonteCarlomethod(int, POINT[]); int main(void) {     POINT *p;     // nの値は100，10000，1000000     int pat[] = {100, 10000, 1000000};     int i, j;     double r;     srand((unsigned)time(NULL));     for (i = 0; i < sizeof(pat) / sizeof(int); i++) {         p = (POINT *)malloc(sizeof(POINT) * pat[i]);         // n個の点(x, y)を乱数で生成させる         for (j = 0; j < pat[i]; j++) {             // 点の座標は0~1の範囲(double)にすること.              // (double)rand() / RAND_MAX;で乱数を生成できる.             p[j].x = (double)rand() / RAND_MAX;             p[j].y = (double)rand() / RAND_MAX;         }         printf("%d times\n", pat[i]);         r = MonteCarlomethod(pat[i], p);         // 円周率の出力は小数点5位まで表示         printf("π ≒ %.5f\n", r);         free(p);     }     return 0; } // 円周率の計算は関数MonteCarlomethodを作っ て行う double MonteCarlomethod(int n, POINT a[]) {     int c = 0, i;     double r;     for (i = 0; i < n; i++) {         r = a[i].x * a[i].x + a[i].y * a[i].y;         if (r <= 1) {             c++;         }     }     return (double)c / n * 4; } ー 実行結果 ー 100 times π ≒ 2.96000 10000 times π ≒ 3.14680 1000000 times π ≒ 3.14320 100 times π ≒ 3.12000 10000 times π ≒ 3.13160 1000000 times π ≒ 3.14120 100 times π ≒ 3.24000 10000 times π ≒ 3.13000 1000000 times π ≒ 3.14090 100 times π ≒ 3.32000 10000 times π ≒ 3.12640 1000000 times π ≒ 3.14382