from cvxopt.glpk import ilp
import numpy as np
from cvxopt import matrix

print help(ilp)

c=matrix(np.ones(6,dtype=float))

print c

coeff=np.array([[1,1,0,0,0,0],
                [1,1,0,0,0,1],
                [0,0,1,1,0,0],
                [0,0,1,1,1,0],
                [0,0,0,1,1,1],
                [0,1,0,0,1,1]
                ],dtype=float)

G=matrix(-coeff)

print G

h=matrix(-1*np.ones(6))

I=set(range(6))

B=set()

print I,B

(status,x)=ilp(c,G,h,matrix(1., (0,6)),matrix(1., (0,1)),I,B)

status

print x

(status,x)=ilp(c,G,h,matrix(1., (0,6)),matrix(1., (0,1)),B,I)

print status
print x