#include <stdlib.h>
#include <stdio.h>
#include <math.h>

#define MAX_N 24

double threshold = 0.8;
int N = 20;
double cr[MAX_N];

void generateCredences() {
    
    double sum = 0;
    
    for (int i=0;i<N;i++) {
        cr[i] = rand()/(double)RAND_MAX;
        sum += cr[i];
    }
    for (int i=0;i<N;i++) {
        cr[i] /= sum;
    }
}

double P(char* p) {
    double sum = 0;
    
    for (int i=0;i<N;i++)
        if (p[i])
            sum += cr[i];
        
    return sum;
}

void generateBelief(char* out) {
    do {
        for (int i=0;i<N;i++) {
            out[i] = rand() % 2;
        }
    }
    while (P(out) < threshold);
}

int closure() {
    char p1[MAX_N];
    char p2[MAX_N];
    generateCredences();
    generateBelief(p1);
    generateBelief(p2);
    for (int i=0;i<N;i++)
        p1[i] = p1[i] && p2[i];
    return P(p1) >= threshold;
}

double calculate() {
    int count=0;
    int n;
    for (n=0;n<10000;n++) {
        count += closure();
    }
    return count / (double)n;
}

int main(int argc, char** argv) {
    if (argc>1)
        threshold = atof(argv[1]);
    printf("N,p\n");
    for (int i=1; i<=MAX_N; i++) {
        N = i;
        printf("%d,%g\n", i, calculate());
    }
    return 0;
}