#include<stdio.h> 

//原文ママ
#include<math.h>
 
	double AttackerSuccessProbability(double q, int z) 
        { 
                double p = 1.0 - q; 
                double lambda = z * (q / p); 
                double sum = 1.0; 
                int i, k; 
                for (k = 0; k <= z; k++) 
                { 
                        double poisson = exp(-lambda); 
                        
			for (i = 1; i <= k; i++)
				poisson *= lambda / i; 
                        
		        sum -= poisson * (1 - pow(q / p, z - k)); 
			
                } 
		
                return sum; 
        } 
//ここまで原文ママ
	

	//ここから先が追加したmain関数
	void main(){
		
		double sumP;//表示されるP値
		double qValue=0.3;//表示されるq値
		int z=0;//表示されるz値
		int num=50;//この値がz値のMAX

		for(z=0; z<=num; z++){
		   sumP= AttackerSuccessProbability(qValue,z);//関数に値を渡す

		   if(z%5==0)//５の倍数で表示
			printf("q=%2.1f   z=%d  Pro= %10.8f\n",qValue,z,sumP);
		

		}

	}
	//ここまで