#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;
    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+OBoOOBleOBhDogIik7CiAgc2NhbmYoIiVsZiIsICZwKTsKICAKICBmb3Ioaj0wO2o8bjtqKyspewogICAgdT0oZG91YmxlKXJhbmQoKS8oUkFORF9NQVgrMS4wKTsKICAgIGMgPSBwLygxLXApOwogICAgaT0wOyAKICAgIHByPSgxLXApOwogICAgZm9yKHo9MTt6PG47eisrKXsKICAgICAgcHIgPSBwciooMS1wKTsKICAgIH0KICAgIGY9cHI7CiAgICB3aGlsZSh1Pj1mKXsKICAgICAgcHI9KGMqKG4taSkvKGkrMSkpKnByOyAKICAgICAgZj1mK3ByOyAKICAgICAgaT1pKzE7CiAgICB9CiAgICB4PWk7CiAgICBmcHJpbnRmKGZwLCIlZFxuIix4KTsKICB9CiAgZmNsb3NlKGZwKTsKfQo=

MTAgMC4wMQ==

10 0.01