using JuMP
using MathProgBase
using GLPKMathProgInterface

m = Model(solver = GLPKSolverMIP())
@variable(m, 0 <= x[i=1:8] <= 9, Int)
@variable(m, y[i=1:8,j=0:9], Bin)

@constraint(m, 1000*x[1] + 100x[2] + 10x[3] + 1x[4] + 1000x[5] + 100x[6] + 10x[7] + 1x[2] == 10000x[5] + 1000x[6] + 100x[7] + 10x[2] + 1x[8])

@constraint(m, [i=1:8], sum([j*y[i,j] for j=0:9]) == x[i])
@constraint(m, [i=1:8], sum([y[i,j] for j=0:9]) == 1)
@constraint(m, [j=0:9], sum([y[i,j] for i=1:8]) <= 1)

@constraint(m, x[1] >= 1)
@constraint(m, x[5] >= 1)

print(m)

status = solve(m)

println("S = ", getvalue(x[1]))
println("E = ", getvalue(x[2]))
println("N = ", getvalue(x[3]))
println("D = ", getvalue(x[4]))

println("M = ", getvalue(x[5]))
println("O = ", getvalue(x[6]))
println("R = ", getvalue(x[7]))

println("Y = ", getvalue(x[8]))