Make a list of years instead of 100 year variables:
years = [lp.LpVariable(str(2011+i), 0, None, lp.LpInteger) for i in xrange(n)]
Note that lists are 0-indexed, so what used to be year_1
is now years[0]
.
You can loop over it for the "declare constraints" part of the script:
for year, next_year in zip(years[:-1], years[1:]):
prob += next_year - year >= 0