aaaa.py

n = 1000
values = [5+(i**2 % 111) for i in range(1, n+1)]
weights = [20+(i % 11) for i in range(1, n+1)]
W=800

print("maximize")
print("-", " + ".join(str(v)+" x("+str(i+1)+")" for i, v in enumerate(values)))
print("subject to")
print("c1:", " + ".join(str(w)+" x("+str(i+1)+")" for i, w in enumerate(weights)), "<=", str(W))
print("bounds")
print("binary")
print(" ".join("x("+str(i)+")" for i in range(1, n+1)))
print("end")