C言語　多次元配列に記載の確率を実行

前提・実現したいこと

x＝乱数として、

配列上の確率に応じたaの行動に対してbの行動をする。そのbの行動を新しいaとして格納し、bの行動をする。というのを繰り返したいです。
配列上の２４列２５列の確率が発生した時に、そのループを抜け出したいです。

プログラミング初心者です。お手柔らかにお願いします。

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

１、確率を記載した２５行２５列多次元配列をコード上に載せる

２、配列からデータを取り入れる

３、確率に応じた行動（０〜２３列に記載の確率ではループ、２４、２５列に記載の確率ではループを抜けだす）を行わせる

該当のソースコード

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

int tpm[22][24] = {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.2, 0.5, 0.3, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.7, 0, 0, 0.05, 0.25, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.1, 0.1, 0.1, 0.2, 0.3, 0.1, 0, 0, 0, 0, 0, 0.1},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.25, 0.5, 0.25, 0, 0, 0, 0, 0, 0, 1, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.6, 0.3, 0, 0, 0, 0, 0.1, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.25, 0.5, 0.25, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.1, 0.5, 0.2, 0, 0, 0, 0, 0, 0, 0, 0.2, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.5, 0.3, 0, 0, 0, 0.2},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.2, 0, 0.2, 0, 0, 0, 0, 0, 0, 0, 0.6, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.2, 0.3, 0.5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.25, 0.5, 0.05, 0, 0, 0, 0, 0, 0, 0, 0.1, 0.1},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.5, 0, 0, 0, 0, 0, 0, 0, 0, 0.3, 0.2},
{0, 0, 0.4, 0, 0.2, 0.4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0.5, 0.5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.5, 0.3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.2, 0},
{0, 0, 0, 0, 0.2, 0, 0.2, 0, 0.2, 0, 0, 0.2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.2},
{0, 0, 0, 0, 0, 0, 0.3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.4, 0.3},
{0, 0, 0, 0, 0, 0.1, 0, 0.3, 0, 0.4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.2, 0},
{0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0.1, 0.3, 0.2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0.2},
{0, 0, 0.1, 0.1, 0, 0, 0, 0.7, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.1, 0},
{0, 0, 0.3, 0.1, 0.1, 0.1, 0.2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.2},
{0, 0, 0.3, 0, 0.5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.1, 0.1}};

int main {

    for(i=0;i<2;i++) {
        x = (double)rand()/RAND_MAX;         
        Ma = 0;
        Mb = 0;
        
        for (j = 1; (x <= Ma) || (Mb < x); j++) {
        
            Ma = Mb;
            Mb = ;  /*ここに新しいbの行動*/

補足情報（FW/ツールのバージョンなど）

行列を記載したいのですが、その場合は２次元配列でよろしいのでしょうか？

また、０〜２３列の確率ではaに対してbの行動をとり、そのbを新たなaとして格納するという考えです。２４、２５列では、ループを抜け出すということです。

説明が下手ですみません。行動という表現には語弊がありました。

確率をPとすると
行に記載されている確率は０ <= P <= １で、同じ行を全て足すと１になるとします。
a（０行X列）の確率を実行すると、次は、b（X行Y列）とします。

０行X列の次はX行Y列、その次はY行Z列、その次は…というように確率で繰り返していき、２４、２５列の確率を引き当てるまでループするという考えです。

遅くなってしまい申し訳ありません。
スタートは０行or１行、１２or１３行とします。
よろしくお願いします。

sh_akira

2018/04/26 08:09

プログラミング以前に日本語が良く分かりません。まずは「c 多次元配列初期化」とかで検索して問題１をご自身で試して、間違っていてもいいので自分なりに書いてみてください。一番意味が分からないのは24,25列の確率でループを抜け出すという部分です。これは逆に言えば確率でループを抜け出さない可能性もあるという事ですか？その場合何が起きるんでしょう。ちょっと分からないです。

coco_bauer

2018/04/26 13:14

「２５行２５列多次元配列」って、２５ｘ２５の2次元配列ではないのですか？　３次元以上に多次元配列になるのなら、その理由と構成方法を説明してください。

sh_akira

2018/04/27 03:14

ただ行列を書くだけなら2次元配列で良いでしょう。「bの行動」ってなんでしょうか。何か計算式を入れるんですか？

sh_akira

2018/04/27 05:34

何となくわかりましたが、25x25の定義を書くのは面倒なので、ご自身で書いてみてください。それから、24列目と25列目の時はループから抜けるのであれば、24,25行目はいらないですよね？23行25列で良いのではないでしょうか？

行動規範の内容に同意します

回答2件

ベストアンサー

c
1#include <stdio.h>
2#include <time.h>
3#include <stdlib.h>
4
5double tpm[23][25] = {
6	{ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.2, 0.5, 0.3, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
7	{ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.7, 0, 0, 0.05, 0.25, 0, 0, 0, 0, 0, 0, 0 },
8	{ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.1, 0.1, 0.1, 0.2, 0.3, 0.1, 0, 0, 0, 0, 0, 0.1 },
9	{ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.25, 0.5, 0.25, 0, 0, 0, 0, 0, 0, 1, 0, 0 },
10	{ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.6, 0.3, 0, 0, 0, 0, 0.1/*, 0*/ },
11	{ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.25, 0.5, 0.25, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
12	{ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.1, 0.5, 0.2, 0, 0, 0, 0, 0, 0, 0, 0.2, 0 },
13	{ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.5, 0.3, 0, 0, 0, 0.2 },
14	{ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.2, 0, 0.2, 0, 0, 0, 0, 0, 0, 0, 0.6, 0 },
15	{ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.2, 0.3, 0.5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
16	{ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
17	{ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.25, 0.5, 0.05, 0, 0, 0, 0, 0, 0, 0, 0.1, 0.1 },
18	{ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.5, 0, 0, 0, 0, 0, 0, 0, 0, 0.3, 0.2 },
19	{ 0, 0, 0.4, 0, 0.2, 0.4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
20	{ 0, 0, 0, 0, 0, 0.5, 0.5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
21	{ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.5, 0.3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.2, 0 },
22	{ 0, 0, 0, 0, 0.2, 0, 0.2, 0, 0.2, 0, 0, 0.2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.2 },
23	{ 0, 0, 0, 0, 0, 0, 0.3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.4, 0.3 },
24	{ 0, 0, 0, 0, 0, 0.1, 0, 0.3, 0, 0.4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.2, 0 },
25	{ 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
26	{ 0, 0, 0, 0, 0, 0.1, 0.3, 0.2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0.2 },
27	{ 0, 0, 0.1, 0.1, 0, 0, 0, 0.7, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.1, 0 },
28	{ 0, 0, 0.3, 0.1, 0.1, 0.1, 0.2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.2 },
29	/*{ 0, 0, 0.3, 0, 0.5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.1, 0.1 }*/
30};
31
32int main()
33{
34	int Ma = 0;
35	int Mb = 0;
36	double y = 0;
37
38	srand((unsigned int)time(NULL));
39
40	while (Ma <= 23) {
41		double x = (double)rand() / RAND_MAX;
42		y = 0;
43		printf("Ma:%d x:%lf\n", Ma, x);
44		for (Mb = 0; Mb < 25; Mb++) {
45			y += tpm[Ma][Mb];
46			if (x <= y) {
47				break;
48			}
49		}
50		Ma = Mb;
51		if (Ma >= 24) {
52			break;
53		}
54	}
55	printf("End Ma:%d\n", Ma);
56	return 0;
57}