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

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

ただいまの
回答率

88.60%

プログラムのエラーの原因がわかりません

解決済

回答 2

投稿

  • 評価
  • クリップ 0
  • VIEW 6,820

jackie687456

score 17

前提・実現したいこと

次のような問題のプログラムで、コンパイルはできるのですが実行するとエラーが出てしまいます。
(cygwinにgccを入れてやっています)

コインを100回投げる試行を行うことを考える。
メトロポリスの方法を使って、二項分布に従うマルコフ過程から10000個分のサンプル列を作成せよ。
ただし、確率変数kとすると、kの次にはk-1かk+1が1/2の確率でそれぞれ選ばれ、0と100はつながっているものとする。

おそらく除算をしているところで不具合が起きているのではないかと考えているのですが・・・。
お手数ですがよろしくお願いいたします。

発生している問題・エラーメッセージ

Floating point exception (コアダンプ)

該当のソースコード

#include<stdio.h>
#include<math.h>
#include<stdlib.h>
#include<stdint.h>

#define K 0

uint32_t xor128(void) {             //乱数の生成(0以上2^32未満)
  static uint32_t x = 123456789;
  static uint32_t y = 362436069;
  static uint32_t z = 521288629;
  static uint32_t w = 88675123; 
  uint32_t t;

  t = x ^ (x << 11);
  x = y; y = z; z = w;
  return w = (w ^ (w >> 19)) ^ (t ^ (t >> 8)); 
}

int fac(int n) {              //階乗の計算
    int result = 1;
    int k;

    for (k = 1; k <= n; k++) {
        result *= k;
    }
    return result;
}

int cmb(int n, int k){        //combinationの計算
    int c;
    c = fac(n) / (fac(k) * fac(n - k));
    return c;
}

int main(void){
    int P[101];
    int W[101] = {0};
    int i, k, l;
    long a;
    double r;

    P[0] = 1;
    P[100] = 1;

    for (i = 1; i <= 99; i++){
        P[i] = cmb(100,i);
    }

    W[K]++;
    k = K;

    for (i = 1; i < 10000; i++){
        a = xor128();
        if (a < pow(2, 31)){
            l = k - 1;
            if (l == -1){
                l = 100;
            }
        }
        else{
            l = k + 1;
            if (l == 101){
                l = 0;
            }
        }

        a = xor128();
        r = (double)P[l] / (double)P[k];
        r = fmin(r, 1);
        if (a < pow(2, 32)*r){
            k = l;
        }

        W[k]++;
    }

    for (i = 0; i <= 100;i++){
        printf("%d\n",W[i]);
    }

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

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

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

    クリップを取り消します

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

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

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

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

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

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

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

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

    質問の評価を下げる

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

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

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

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

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

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

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

    詳細な説明はこちら

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

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

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

回答 2

checkベストアンサー

+2

途中の計算結果が大きすぎて、intで表現できていないのではないでしょうか?

こちらのサイトによると、100C50は100891344545564193334812497256だそうです。
また、combinationの計算方法も、先に分子と分母を求めて除算結果を返しておられるようですが、ここで必要になる100!はさらに大きく93326215443944152681699238856266700490715968264381621468592963895217599993229915608941463976156518286253697920827223758251185210916864000000000000000000000000です。
こちらのサイトで紹介されているような、漸化式による方法がいいかもしれません。

投稿

  • 回答の評価を上げる

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

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

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

  • 回答の評価を下げる

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

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

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

  • 2017/07/22 21:55

    無事解決しました。ありがとうございました。

    キャンセル

0

web関係以外は全くわからないのですが、ちらっとint型は扱える数が限られていると聞いたことがあります。
Bongoさんと同じく、intが使えないものと思います。

投稿

  • 回答の評価を上げる

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

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

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

  • 回答の評価を下げる

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

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

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

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

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

関連した質問

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