pok

#include <iostream>
#include <random>
#include <algorithm>

using namespace std;


constexpr size_t N = 1000000000;

int main(){
    random_device rd;
    mt19937 gen(rd());

    vector <int_fast8_t>  talia(9+38);
    generate(begin(talia),end(talia),[n=0]()mutable{ return (n++)<9?1:0 ;});

    vector <size_t> wins(10,0);
    size_t cel = 1000;

    constexpr double ex =  1.0-38.0*37.0/(46.0*47.0);

    for (size_t i=0; i<=N; i++){

        std::shuffle(begin(talia),end(talia),gen);

        for (int gracz=0; gracz<10; gracz++){
            if (talia[2*gracz]==1 || talia[2*gracz+1]==1)
                wins[gracz]++;
        }

        if (i==cel){
            cout<<i<<" ";
            for (int gracz=0; gracz<10; gracz++){
                cout<<wins[gracz]/(double)i<<" ";
            }
            double err = 3.0*sqrt( ex*(1-ex)/i );
            cout<< "  ["<<ex-err<<" ,"<<ex+err<<"]\n";
            cel*=10;
        }


    }

    return 0;
}