#include <stdio.h>
#include <stdlib.h>
#include <assert.h>
#include <math.h>
typedef double R;
#define d 50 
#define k 1000000
R scale = 0x1.p-31;
R rf(){return 0x1.p-31*random();}
int main(){int his[d]; for(int j=0; j<d; ++j) his[j]=0;
  printf("  c     PDF(c)    his    den(c)  \n");
  for(int j = 0; j<k; ++j) {int z = rf()*rf()*d; assert(0<= z && z < d); ++his[z];}
  for(int j = 0; j<d; ++j) {R c = (0.5 + (R)j)/d; printf("%5.3f %10.7f %5d %10.7f\n", 
           c,  -log(c), his[j], d*((R)his[j]/k));}
  return 0;}
/*
  c     PDF(c)    his    den(c)  
0.010  4.6051702 98107  4.9053500
0.030  3.5065579 70807  3.5403500
0.050  2.9957323 59763  2.9881500
0.070  2.6592600 53184  2.6592000
0.090  2.4079456 47994  2.3997000
0.110  2.2072749 43816  2.1908000
0.130  2.0402208 41027  2.0513500
0.150  1.8971200 38077  1.9038500
0.170  1.7719568 35866  1.7933000
0.190  1.6607312 33214  1.6607000
0.210  1.5606477 31093  1.5546500
0.230  1.4696760 29330  1.4665000
0.250  1.3862944 27836  1.3918000
0.270  1.3093333 26385  1.3192500
0.290  1.2378744 24459  1.2229500
0.310  1.1711830 23341  1.1670500
0.330  1.1086626 22282  1.1141000
0.350  1.0498221 20853  1.0426500
0.370  0.9942523 20108  1.0054000
0.390  0.9416085 18849  0.9424500
0.410  0.8915981 18029  0.9014500
0.430  0.8439701 16869  0.8434500
0.450  0.7985077 15893  0.7946500
0.470  0.7550226 15135  0.7567500
0.490  0.7133499 14380  0.7190000
0.510  0.6733446 13533  0.6766500
0.530  0.6348783 12773  0.6386500
0.550  0.5978370 12023  0.6011500
0.570  0.5621189 11217  0.5608500
0.590  0.5276327 10489  0.5244500
0.610  0.4942963  9893  0.4946500
0.630  0.4620355  9382  0.4691000
0.650  0.4307829  8512  0.4256000
0.670  0.4004776  7879  0.3939500
0.690  0.3710637  7352  0.3676000
0.710  0.3424903  6734  0.3367000
0.730  0.3147107  6339  0.3169500
0.750  0.2876821  5719  0.2859500
0.770  0.2613648  5072  0.2536000
0.790  0.2357223  4856  0.2428000
0.810  0.2107210  4196  0.2098000
0.830  0.1863296  3640  0.1820000
0.850  0.1625189  3277  0.1638500
0.870  0.1392621  2840  0.1420000
0.890  0.1165338  2348  0.1174000
0.910  0.0943107  1925  0.0962500
0.930  0.0725707  1482  0.0741000
0.950  0.0512933  1013  0.0506500
0.970  0.0304592   614  0.0307000
0.990  0.0100503   195  0.0097500
*/