fork(1) download
  1. #include<stdio.h>
  2. #include<stdlib.h>
  3. #include<math.h>
  4.  
  5. #define lambda 1.0
  6.  
  7. main(){
  8. FILE *fp;
  9. double u,p,f,c,pr;
  10. int n,i,j,x,z;
  11.  
  12. fp = stdout; /*fopen("niko.csv","wt");*/
  13.  
  14. if(fp==NULL){
  15. fprintf(stderr,"cannot open file !\n");
  16. exit(1);
  17. }
  18. printf("試行回数を入力してください: ");
  19. scanf("%d", &n);
  20. printf("成功確率を入力してください: ");
  21. scanf("%lf", &p);
  22.  
  23. for(j=0;j<n;j++){
  24. u=(double)rand()/(RAND_MAX+1.0);
  25. c = p/(1-p);
  26. i=0;
  27. pr=(1-p);
  28. for(z=1;z<n;z++){
  29. pr = pr*(1-p);
  30. }
  31. f=pr;
  32. while(u>=f){
  33. pr=(c*(n-i)/(i+1))*pr;
  34. f=f+pr;
  35. i=i+1;
  36. }
  37. x=i;
  38. fprintf(fp,"%d\n",x);
  39. }
  40. fclose(fp);
  41. }
  42.  
Success #stdin #stdout 0.01s 2008KB
stdin
1000 0.95
stdout
試行回数を入力してください: 成功確率を入力してください: 957
948
955
956
959
944
947
955
946
951
950
952
948
950
961
959
953
954
943
952
935
945
942
956
943
948
942
941
970
945
950
957
952
946
953
951
950
963
946
955
951
955
948
958
946
948
956
959
940
961
951
940
944
953
958
947
939
935
949
939
945
963
959
957
946
951
948
955
950
953
951
938
949
960
960
954
946
954
953
948
953
943
949
958
957
947
945
958
947
953
961
952
953
957
949
960
948
956
953
959
950
945
961
960
943
958
953
949
952
946
956
947
949
945
944
946
951
949
943
959
941
942
950
955
964
960
953
948
955
948
946
945
952
945
943
954
942
956
943
955
940
961
939
951
944
945
956
954
953
962
953
955
941
942
950
940
940
944
949
956
951
955
939
943
977
944
958
942
968
939
958
940
931
960
952
944
943
948
959
956
948
951
952
949
953
941
951
955
947
966
951
958
955
952
937
955
957
960
958
957
964
955
959
964
953
950
943
957
958
940
953
945
952
945
954
947
947
945
940
952
945
953
950
963
946
951
954
942
950
952
961
949
947
957
949
930
947
952
957
945
953
950
950
947
954
944
952
938
949
954
953
953
947
944
952
952
954
947
955
944
960
953
953
946
951
941
946
958
953
941
942
948
951
952
953
946
955
946
947
944
964
931
957
947
944
949
962
959
955
954
942
953
948
955
961
959
957
944
956
951
946
959
959
958
950
951
943
946
958
950
950
957
950
946
944
953
954
943
952
950
957
954
944
945
950
942
942
948
947
960
959
952
957
956
950
947
948
953
952
946
943
940
941
953
948
946
947
941
949
960
952
946
953
955
950
943
958
952
950
944
951
949
957
948
945
947
954
953
964
947
959
942
949
932
955
955
946
942
958
954
938
949
965
954
945
950
958
941
941
948
947
953
956
942
939
939
949
955
954
949
942
952
952
951
952
948
947
958
950
943
952
951
952
952
945
957
940
941
960
943
950
945
957
946
948
958
947
956
953
937
946
955
952
957
947
945
944
952
941
953
955
954
944
948
947
949
945
948
951
943
951
934
948
948
947
954
946
953
951
960
953
956
954
947
953
932
951
957
952
953
950
948
948
954
956
953
943
937
939
953
944
952
954
951
928
932
947
946
953
957
944
947
953
958
953
947
958
944
943
950
957
953
959
944
948
954
958
951
954
960
945
960
951
960
950
951
960
956
956
952
953
967
960
947
958
952
953
955
955
956
946
952
950
943
956
958
958
953
949
952
952
953
951
943
952
937
954
950
957
950
942
950
950
938
956
948
953
949
943
949
945
949
935
954
951
956
952
949
950
967
939
940
953
952
945
945
952
960
954
950
949
957
961
961
959
955
947
951
945
950
961
950
958
962
944
949
955
955
959
946
955
962
947
948
951
951
952
944
950
948
953
959
945
952
958
941
948
944
953
952
953
952
939
951
951
945
935
947
933
960
952
955
958
960
944
949
950
956
953
962
943
946
958
948
959
955
950
946
961
942
958
952
954
960
944
946
943
944
952
944
942
945
961
934
943
942
949
953
959
941
952
940
948
950
949
946
945
959
951
944
937
949
956
955
948
969
937
951
944
953
954
947
958
953
947
938
955
956
954
953
959
947
955
946
956
944
952
939
941
943
945
942
952
939
958
961
939
960
950
946
952
943
952
950
956
960
951
952
954
945
946
952
951
935
960
947
944
951
948
947
953
953
949
946
954
947
945
955
946
953
935
957
957
952
947
953
951
957
945
946
940
950
958
952
951
956
962
954
947
948
938
936
931
950
946
954
956
950
950
940
944
950
960
938
940
945
954
952
941
962
958
943
950
955
955
932
952
954
955
960
941
956
961
941
946
945
956
941
955
946
944
960
955
941
964
957
947
954
949
949
953
947
952
942
940
948
942
953
942
958
952
945
953
951
947
961
955
942
937
950
948
945
949
943
947
949
932
953
942
949
940
955
956
953
959
958
939
938
951
944
960
941
948
952
953
954
951
949
956
951
960
945
956
948
948
941
956
948
955
960
957
957
954
953
950
952
950
952
953
937
955
952
943
943
944
955
957
954
944
953
947
941
958
941
950
946
944
946
953
961
945
950
955
960
943
946
951
953
958
945
953
952
957
956
955
937
952
952
955
955
946
940
958
944
944
948
949
948
953
941
947
958
952
941
956
955
948
948
948
946
952
940
959
949
958
953
949
950
946
945
945
951
947
942
955
950
950
944
957
942
946
944
930
959
946
956
953
953
944
938
962
955
942
957
944
967
951
953
950
957
958
954
948
944
957
943
953
947
947
951