#include<stdio.h>
#include<stdlib.h>
#include<math.h>
 
#define lambda 1.0
 
main(){
  FILE *fp;
  double u,p,f,c,pr;
  int n,i,j,x,z;
 
  fp = stdout; /*fopen("niko.csv","wt");*/
 
  if(fp==NULL){
    fprintf(stderr,"cannot open file !\n");
    exit(1);
  }
  printf("試行回数を入力してください: ");
  scanf("%d", &n);
  printf("成功確率を入力してください: ");
  scanf("%lf", &p);
 
  for(j=0;j<n;j++){
    u=(double)rand()/(RAND_MAX+1.0);
    c = p/(1-p);
    i=0; 
    pr=(1-p);
    for(z=1;z<n;z++){
      pr = pr*(1-p);
    }
    f=pr;
    printf("%g\n",pr);
    fflush(stdout);
    while(u>=f){
      pr=(c*(n-i)/(i+1))*pr; 
      f=f+pr; 
      i=i+1;
    }
    x=i;
    fprintf(fp,"%d\n",x);
  }
  fclose(fp);
}
 

I2luY2x1ZGU8c3RkaW8uaD4KI2luY2x1ZGU8c3RkbGliLmg+CiNpbmNsdWRlPG1hdGguaD4KCiNkZWZpbmUgbGFtYmRhIDEuMAoKbWFpbigpewogIEZJTEUgKmZwOwogIGRvdWJsZSB1LHAsZixjLHByOwogIGludCBuLGksaix4LHo7CiAgCiAgZnAgPSBzdGRvdXQ7IC8qZm9wZW4oIm5pa28uY3N2Iiwid3QiKTsqLwogIAogIGlmKGZwPT1OVUxMKXsKICAgIGZwcmludGYoc3RkZXJyLCJjYW5ub3Qgb3BlbiBmaWxlICFcbiIpOwogICAgZXhpdCgxKTsKICB9CiAgcHJpbnRmKCLoqabooYzlm57mlbDjgpLlhaXlipvjgZfjgabjgY/jgaDjgZXjgYQ6ICIpOwogIHNjYW5mKCIlZCIsICZuKTsKICBwcmludGYoIuaIkOWKn+eiuueOh+OCkuWFpeWKm+OBl+OBpuOBj+OBoOOBleOBhDogIik7CiAgc2NhbmYoIiVsZiIsICZwKTsKICAKICBmb3Ioaj0wO2o8bjtqKyspewogICAgdT0oZG91YmxlKXJhbmQoKS8oUkFORF9NQVgrMS4wKTsKICAgIGMgPSBwLygxLXApOwogICAgaT0wOyAKICAgIHByPSgxLXApOwogICAgZm9yKHo9MTt6PG47eisrKXsKICAgICAgcHIgPSBwciooMS1wKTsKICAgIH0KICAgIGY9cHI7CiAgICBwcmludGYoIiVnXG4iLHByKTsKICAgIGZmbHVzaChzdGRvdXQpOwogICAgd2hpbGUodT49Zil7CiAgICAgIHByPShjKihuLWkpLyhpKzEpKSpwcjsgCiAgICAgIGY9ZitwcjsgCiAgICAgIGk9aSsxOwogICAgfQogICAgeD1pOwogICAgZnByaW50ZihmcCwiJWRcbiIseCk7CiAgfQogIGZjbG9zZShmcCk7Cn0K

MTAwMCAwLjk1

1000 0.95