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


long double faktorial(int a);
float mocnina(float co,int neco);
float ztraty(int N, float zat);
float zatizeni(int N,float ztr);

int main(int x, char *b[]){

int N,q,o;

float h;

 float *f;


N=atoi(b[1]); // ulozeni poctu vedeni

 f = (float *) malloc((x - 2) * sizeof(float));
    if (f == NULL) {
        printf("chyba\n");
        return 1;
    }


printf("ztraty:");

 for(q=2;q<x;q++){ 
f[q]=atof(b[q]);	// ulozeni parametru do pole
printf("	%9.3lf",f[q]);


}


printf("\nvedeni N\n");


for(q=1;q<=N;q++){	// print vedeni
 printf("%2d",q);
    for(o=2;o<x;o++){
    h=zatizeni(q,f[o]); // vypocet zatizeni do tabulky
    printf("	%8.2f",h);
    
    }    
    printf("\n");
}

free(f);

}



 // mocnina

float mocnina(float co,int neco){
    int l;
    float vysledek=co;
    if(neco==0) return 1;
    if(neco==1) return co;

    for(l=1;l<=neco;l++){
    vysledek=vysledek*co;
    }

return vysledek;

}

// faktorial

long double faktorial(int a){

    long double x;
    int b;

    x=1;

    if(a==0) return 1;
    if(a==1) return 1;


    for(b=1;b<=a;b++){
    x=b*x;

    }
    return x;
}


// ztraty
// funguje do cca N=100

float ztraty(int N, float zat) {
    int i;
    float suma = 0, zt, y = 1;
    
    for (i = 1; i <= N; i++) {
        y =y*zat / i; 
        suma = suma+y;
    }
    
    zt = y;
    zt = zt / suma;
    
    return zt;
}


 // zatizeni

float zatizeni (int N, float ztr){
  int t;
    float prav, lev, z, r;
    lev = N;
    prav = 0;

    for (t = 0; t < 200; t++) {
        r = (prav + lev) / 2;
        z = ztraty(N, r);

//printf(cislo %d r= %f, z= %f\n",t,r,z);
        if (z < ztr)
            prav = r;
        else
            lev = r;
    }

    return prav;
}