質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

ただいまの
回答率

88.80%

1000000回繰り返しすべての試行回数を足すプログラム

解決済

回答 4

投稿

  • 評価
  • クリップ 0
  • VIEW 577

Shotayanagi

score 5

C言語 プログラミングに関する質問です。

0~9までの番号が振られたくじがあり、出る確率はすべて等しいものとするとくじ10本をすべて引き切るまでに何回くじを引く必要があるか求め、そのシュミレーションを1000000回行った時の回数の合計と平均値を求めよ。ただしくじは何度でも引けるものとし、rand関数を用いることとする。

という問いに関してなのですが、以下のプログラムを実行すると合計が18億前後、平均は1800ほどになってしまい、本来平均20~40ほどとなるはずなのに大きくずれが生じてしまいます。個人的に不安なのが変数を宣言する位置とそもそもdo文をこのような形でfor文に取り込むことが出来るかという点、また1000000回シュミレーションしたときの合計(sum)を求めるところです。
まだプログラミングを始めて間もなく、わからないことばかりなので教えていただけると幸いです。。。
よろしくお願いします。

#include<stdio.h>
#include<time.h>
#include<stdlib.h>


/*0~9の数字を返す*/
int kuji(){

return (rand()>>2) % 10;
}


int main(){

int i;
double avg=0;        /*1000000回試行したときの平均値*/
unsigned short Flags = 0; /*各値のくじを引いたか下位10ビットで記憶する変数*/
unsigned int N = 0;    /*くじを引いた回数カウンター*/
unsigned sum=0;       /*1000000回引いた時の回数の合計*/

/*for文にdo文を入れdo文を1000000回繰り返す*/
for(i=0; i<1000000; i++){

sum += N;      /*1000000回分のNの合計*/
avg = sum/1000000; /*1000000回引いた時の平均値*/

srand((unsigned) time(NULL));

/*0~9のくじを引き切るまでの回数を求める*/
do{
++N;
Flags |= (1 << kuji() );
}while( Flags != 0x3FF);

}

printf("sum = %u\n",sum);   /*1000000回引いた時の回数の合計を表示*/
printf("avg = %.2f\n", avg ) ;/*1000000回引いた時の平均の回数を表示*/

return 0;
}
  • 気になる質問をクリップする

    クリップした質問は、後からいつでもマイページで確認できます。

    またクリップした質問に回答があった際、通知やメールを受け取ることができます。

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 過去に投稿した質問と同じ内容の質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

質問への追記・修正、ベストアンサー選択の依頼

  • y_waiwai

    2020/07/16 00:28

    どうなるのが正解なんでしょうか

    キャンセル

  • Shotayanagi

    2020/07/16 00:55

    1000000回シュミレーションを繰り返した時、合計で何回くじを引いたか、とその合計を1000000で割った平均値を表示するという問題です。
    分かりずらくてすみません…
    だいたい平均値は20-40前後、合計はその1000000の値が表示されれば問題なしです…

    キャンセル

回答 4

checkベストアンサー

+1

出題通りのコードになっていません。

  • sumの加算タイミングがおかしい。
    試行回数を求める前に加算しています。最後のループ時の試行回数が加算されません。
  • ループごとにN,Flagsを0クリアしていない。
    2回目以降は1度しか振らず、N前ループの値+1になる。
    初回ループで10回の試行ですべての目が出たとしても、それだけでsum1793293664(sum加算タイミング修正前は1792293655)になります。
  • ループごとにsrandで初期化している。
    1秒以内にループが終わった場合次のループで使われる乱数の値が前のループと同じになる。
  • avgの計算が整数演算になっている。
    小数点以下が求まりません。
    srand((unsigned)time(NULL));
    for (i = 0; i < 1000000; i++) {
        N = 0;
        Flags = 0;
        do {
            ++N;
            Flags |= (1 << kuji());
        } while (Flags != 0x3FF);
        sum += N;
    }
    avg = sum / 1000000.0;


結果

sum = 29282182
avg = 29.28

投稿

編集

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2020/07/16 10:33

    ご指摘ありがとうございます!
    丁寧に実行結果までつけてくださって助かりました!

    キャンセル

+1

srand((unsigned) time(NULL));

これは最初の一回だけ、実行するようにしましょう。
1ms以内にループが回ると、同じ乱数列が出てきてしまいます

また、

avg = sum/1000000; /*1000000回引いた時の平均値*/

これがループの中にあるのが意味不明です。

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2020/07/16 01:03

    >1ms以内
    1秒以内ですよね

    キャンセル

  • 2020/07/16 01:06

    ああ、C言語のヤツはそうですね。
    ご指摘ありがとうございます

    キャンセル

  • 2020/07/16 10:34

    なるほど。ご回答ありがとうございます!

    キャンセル

+1

#include <stdio.h>
#include <time.h>
#include <stdlib.h>

int kuji(){
  return (rand() >> 2) % 10;
}

int main(void)
{
  srand((unsigned)time(NULL));

  const int imax = 1000000;
  unsigned int sum = 0;

  for (int i = 0; i < imax; i++) {
    unsigned short Flags = 0;
    unsigned int n;
    for (n = 0; Flags != 0x3ff; n++) { /*# 修正 */
      Flags |= 1 << kuji();
    }
    sum += n;
  }

  double avg = (double)sum / imax;

  printf("sum = %u\n", sum);
  printf("avg = %.2f\n", avg);

  return 0;
}

投稿

編集

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2020/07/16 01:19 編集

    >for (n = 1; Flags != 0x3ff; n++)
    一周目ループ終了時にnが2になるので、初期値を0にしないと一回余分にカウントしますよ。

    キャンセル

  • 2020/07/16 01:21

    おっと、直します。

    キャンセル

  • 2020/07/16 10:36

    ご丁寧に全文つけてくださり、ありがとうございます!
    1000000も文字に代入してしまえばだいぶ見やすくなりまりますね!
    助かりました!

    キャンセル

-3

出る回数が等しい乱数をシミュレーションしたいのであればこの場合は乱数としてi%10を疑似乱数列としてつかってください。そうするとご希望どおりの出現率がおおよそ等しくなります。
シードごとの乱数の分布を知りたいのであれば、シードも保存してください。

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

15分調べてもわからないことは、teratailで質問しよう!

  • ただいまの回答率 88.80%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる

関連した質問

同じタグがついた質問を見る