#include #include #include #include #include #include //Obliczamy silnie z argumentu a long double silnia (int a) { long double silnia=1; int i; if (a==0) { return silnia; } else { for ( i=1 ; i<=a ; i++) { silnia=silnia*i; } } return silnia; } //Obliczamy n nad k long double n_nad_k (int n, int k) { long double mnon=1,s; int wybor,i; if(n-k>k) { wybor=n-k;s=silnia(k);} else {wybor=k;s=silnia(n-k);} for(i=wybor+1;i<=n;i++) { mnon*=i; } return mnon/s; } //Obliczamy prawdopodobienstwo long double prawd (int k0, int k1) { long double sum=0,nnadk; int n=61,i; for (i=k0;i<=k1;i++) { nnadk=n_nad_k(n,i); sum+=nnadk*(long double)pow((0.33),i)*(long double)pow((0.33),(91-i)); } return sum; } //Obliczamy parametr alfa long double alfa (int k) { long double sum=0,nnadk; int n=91,i; for (i=0;i<=k;i++) { nnadk=n_nad_k(n,i); sum+=nnadk*(long double)pow((0.33),i)*(long double)pow((0.33),(91-i)); } return sum; } //Obliczamy parametr beta long double beta (int k) { long double sum=0,nnadk; int n=91,i; for (i=k;i<=n;i++) { nnadk=n_nad_k(n,i); sum+=nnadk*(long double)pow((0.33),i)*(long double)pow((0.33),(91-i)); } return sum; } int main(int argc, char **argv) { FILE *plik1; int i,j,n=91; long double delta=0.07; plik1 = fopen("dane", "w"); int k0max=0,k1min=91; //Sprzawdzamy nasze warunki for (i=0;i(1-delta))&&(alfa(i-1)<=(delta/2))&&(beta(j+1)<=(delta/2))){ if(k0maxj) {k1min=j;} fprintf(plik1, " TAK k0:%d k1:%d Pr>1-delta (%.10Lf > %Lf) alfa