#include <stdio.h>
#include <stdlib.h>
//#include <math.h>
#define KONECPOLE 500

long double faktorial (int f);
double ztraty (int N,float y);
long double mocnina(float to,int nato);
float zatizeni (int N,double z);

int main(int PP, char *par[]){
int N,i,j;
double pole[KONECPOLE+1];
float y;
float ztraty;
float z;

N=atoi(par[1]);
if(PP>KONECPOLE){			//kontrola
    printf("je zadano moc parametru!!!!!\n");
    return 1;
}

printf("\n%Lf\n",faktorial(N));

printf("ztraty:		");

for(i=2;i<PP;i++){			//print vrchni listy
    //printf("i je:%d\n",i);
    pole[i]=atof(par[i]);
    printf("%2.3lf	",pole[i]);
}

printf("\npocet\nvedeni\n");		//makeup



for(i=1;i<=N;i++){			//for podle poctu radek
    printf("	%d",i);
    for(j=2;j<PP;j++){			//for podle poctu sloupcu
	y=zatizeni(i,pole[j]);
	printf("	%2.2f",y);
    }
    printf("\n");
}

//y=zatizeni(N,z);

//printf("zatizeni je:\n%f\n",y);

}
float zatizeni (int N,double z){		//vypocet zatizeni pomoci puleni intervalu
    int i;
    float ly=N,py=0,vz,y;
    for(i=0;i<1000;i++){
	y=(py+ly)/2;
	vz=ztraty(N,y);
	if (vz==z) break;
	if(vz<z)py=y;
	else ly=y;
	//printf("%f\n",y);
    }
    
    return py;
}


double ztraty (int N,float y){
    double suma;
    long double z,hor,dol;
    int i;
    dol=0;
    for(i=0;i<=N;i++){
	suma=mocnina(y,i)/faktorial(i);
	dol=dol+suma;
    }
    hor=mocnina(y,N)/faktorial(N);
    z=hor/dol;
    return z;

}
long double mocnina(float to,int nato) {
    long double a=to;
    int i;
    if (nato==0) return 1;
    if (nato==1) return to;
    
    for (i=1;i<nato;i++) a=a*to;
    return a;
}


long double faktorial (int f){
if(f<0)return 0;
if(f == 0)return 1;
if(f == 1)return 1;
if(f>1){
    int i;
    long double x=1;
    for(i=1;i<=f;i++){
	x=x*i;
    }
    return x;
}
}