using JuMP
using CPLEX


edges = [  ...]

model = Model(solver = CplexSolver())

n, m = size(edges)

@variable(model, 0 <= x[1:n] <= 1)
@variable(model, 0 <= y[1:m] <= 1)

@constraint(model, result_constraint,  sum(x[i] for i = 1:n) >= sum(y[i] for i = 1:m) + 1)


for j in 1:m
    for i in 1:n
        if (edges[i,j] == 1)    # constraint for edges
            @constraint(model, x[i] <= y[j])
        end
    end
end

@objective(model, Min, sum(x[i] for i = 1:n))


status = solve(model)
println(status)

r= getobjectivevalue(model)
println(r)
selected   =getvalue(x)
neighbors  =getvalue(y)
println(selected)
println(neighbors)