#include<iostream>
#include <ctime>
#include <cstdlib>
using namespace std;
double tab[5][5];
double Px[5];
double Fx[5];
double Pywarx[5][5];
double Fywarx[5][5];
double random1,random2;
int x, y;
int ile[5][5];
int main()
{
srand(time(0));
int size = 5;
tab[1][1] = 0.2;
tab[1][3] = 0.4;
tab[2][1] = 0.1;
tab[3][4] = 0.05;
tab[4][1] = 0.2;
tab[4][3] = 0.05;
for (int k = 1; k < 5; k++)
{
for (int z = 1; z < 5; z++)
Px[k] += tab[k][z];
}
for (int k = 1; k < 5; k++)
Fx[k] = Fx[k - 1] + Px[k];
for (int k = 1; k < 5; k++)
for (int z = 1; z < 5; z++)
Pywarx[k][z] = tab[z][k] / Px[z]; //najpierw y
for (int k = 1; k < 5; k++)
for (int z = 1; z < 5; z++)
Fywarx[k][z] = Fywarx[k - 1][z] + Pywarx[k][z];
for (int k = 1; k <= 100000; k++)
{
random1 = ((double)rand() / (RAND_MAX));
for (int z = 1; z < 5; z++)
if (random1 <= Fx[z])
{
x = z;
break;
}
random2 = ((double)rand() / (RAND_MAX));
for (int z = 1; z < 5; z++)
if (random2 <= Fywarx[z][x])
{
y = z;
break;
}
ile[x][y] += 1;
}
for (int k = 1; k < 5; k++)
{
for (int z = 1; z < 5; z++)
cout << ile[k][z] << " ";
cout << endl;
}
return 0;
}
{"html5":"htmlmixed","css":"css","javascript":"javascript","php":"php","python":"python","ruby":"ruby","lua":"text\/x-lua","bash":"text\/x-sh","go":"go","c":"text\/x-csrc","cpp":"text\/x-c++src","diff":"diff","latex":"stex","sql":"sql","xml":"xml","apl":"apl","asterisk":"asterisk","c_loadrunner":"text\/x-csrc","c_mac":"text\/x-csrc","coffeescript":"text\/x-coffeescript","csharp":"text\/x-csharp","d":"d","ecmascript":"javascript","erlang":"erlang","groovy":"text\/x-groovy","haskell":"text\/x-haskell","haxe":"text\/x-haxe","html4strict":"htmlmixed","java":"text\/x-java","java5":"text\/x-java","jquery":"javascript","mirc":"mirc","mysql":"sql","ocaml":"text\/x-ocaml","pascal":"text\/x-pascal","perl":"perl","perl6":"perl","plsql":"sql","properties":"text\/x-properties","q":"text\/x-q","scala":"scala","scheme":"text\/x-scheme","tcl":"text\/x-tcl","vb":"text\/x-vb","verilog":"text\/x-verilog","yaml":"text\/x-yaml","z80":"text\/x-z80"}