#define N 16
#define C 15
#define TOT 50
 
int vars[TOT],parity[TOT],begin[C+1];
continuous X[N];

2lp_main()
{
	data();
	setup();
	
	and(int j=0;j<N;j++)			// Loop through variables
		or(int b=0;b<2;b++) 			// b for bit
			X[j] == b;					// Truth value of X[j] is b

	output();
}

data()
{
vars = {
1,6,0,
0,15,3,2,
4,15,14,2,7,
7,6,12,4,
7,3,
13,6,7,8,
1,9,5,
2,6,
11,0,3,
7,9,2,
11,4,13,14,
0,2,3,7,
11,12,13,
11,14,15,
1,0,13
};

parity = {
1,0,0,
0,1,0,1,
0,0,1,1,1,
0,1,0,1,
1,1,
1,0,1,1,
0,1,0,
1,0,
0,1,1,
0,0,1,
0,1,0,0,
1,1,1,1,
1,1,1,
0,0,0,
0,1,0
};

begin = {
0,
3,
7,
12,
16,
18,
22,
25,
27,
30,
33,
37,
41,
44,
47,
TOT
 };
}

setup()
{
	and(int j=0;j<N;j++)
		X[j] <= 1;

	and(int i=0;i<C;i++)
		1 <= sigma(int k=begin[i];k<begin[i+1];k++) 
			( 
			  parity[k]*X[vars[k]] + 
			      (1 - parity[k])*(1-X[vars[k]])
			); 
}

output()
{
		and(int i=0;i<N;i++)
			printf("Variable %d was set to %.0f \n",i,X[i]);
}
 