pdf = new double[B + 2];
cdf = new double[B + 2];
memset(pdf, 0, sizeof(double) * (B + 2));
memset(cdf, 0, sizeof(double) * (B + 2));
for (int i = 1; i <= B; i++) {
pdf = new double[B + 2];
cdf = new double[B + 2];
memset(pdf, 0, sizeof(double) * (B + 2));
memset(cdf, 0, sizeof(double) * (B + 2));
for (int i = 1; i <= B; i++) {