#include <stdlib.h>

#include <stdio.h>

#include <math.h>

#include <time.h> 

​

char fileName[200];

FILE *myFile;

int K, cols;

long double * binomialProbability;//pow(2,K)

long double p;

int * temp, * zeta;//K * pow(2,K)

​

void allocatePatternsWProbability(void);

void generateAllBinaryStrings(int col, int row);

void putInZeta(int i);

void printStrings(void);

​

/***************************************************************************************************/

int main( int argc , char * argv[] ){

    int i, j, offset;

​

    if ( argc != 3 ){

        printf("Usage: K p\n");

        exit(-9);

    }

    K = (int) atoi(argv[1]);

    p = (long double) atof(argv[2]);

    cols = pow(2,K);

​

    allocatePatternsWProbability();

​

    sprintf(fileName, "binaryStrings_K%d_p%Lg.dat", K, p);

    myFile=fopen(fileName,"w");

    if(myFile == NULL){

        printf("myFile fopen failed ABORT\n");

        exit(-9);

    }

​

    fprintf(myFile,"# binary_strings.c generates binary strings of lenght K\n");

    fprintf(myFile,"# along with their binomial probability in the case there is a probability p of having a +1 bit.\n");

​

    for(i=0; i<cols; i++){

        fprintf(myFile, "%Lg ", binomialProbability[i]);

    }

    fprintf(myFile, "\n");

    for(i=0; i<cols; i++){

        for(j=0; j<K; j++){

            offset = i + cols * j;

            fprintf(myFile, "%d ", zeta[offset]);

        }   

        fprintf(myFile, "\n");

    }

​

    fclose(myFile);

​

    return 0;

}

​

/***************************************************************************************************/

​

void allocatePatternsWProbability(void){

    int i, j, sum, offset, n; 

​

    temp = (int *)malloc(K * sizeof(int));

    if(temp == NULL){

        printf("temp malloc failed ABORT\n");

        fprintf(myFile, "temp malloc failed ABORT\n");

        exit(-9);

    }

    zeta = (int *)malloc(K * cols * sizeof(int));

    if(zeta == NULL){

        printf("zeta malloc failed ABORT\n");

        fprintf(myFile, "zeta malloc failed ABORT\n");

        exit(-9);

    }

​

    binomialProbability = (long double *)malloc(cols * sizeof(long double));

    if(binomialProbability == NULL){

        printf("binomialProbability malloc failed ABORT\n");

        fprintf(myFile, "binomialProbability malloc failed ABORT\n");