Implementing a constraint based on previous variable's value in GNU Mathprog/AMPL

The implication you want to implement is:

 x(i,t) = 1 ==> x(i,t+1) = 1, x(i,t+2) = 1

AMPL supports implications (with the ==> operator), so we can write this directly. MathProg does not.

A simple way to implement the implication as straightforward linear inequalities is:

x(i,t+1) >= x(i,t)    
x(i,t+2) >= x(i,t)

This can easily be expressed in AMPL, MathProg, or any modeling tool.

This is the pure, naive translation of the question. This means however that once a single x(i,t)=1 all following x(i,t+1),x(i,t+2),x(i,t+3)..=1. That could have been accomplished by just the constraint x(i,t+1) >= x(i,t).

A better interpretation would be: we don't want very short run lengths. I.e. patterns: 010 and 0110 are not allowed. This is sometimes called a minimum up-time in machine scheduling and can be modeled in different ways.

1. Forbid the patterns 010 and 0110:

    (1-x(i,t-1))+x(i,t)+(1-x(i,t+1)) <= 2
    (1-x(i,t-1))+x(i,t)+x(i,t+1)+(1-x(i,t+2)) <= 3
   

2. The pattern 01 implies 0111:

     x(i,t+1)+x(i,t+2) >= 2*(x(i,t)-x(i,t-1))
   

Both these approaches will prevent patterns 010 and 0110 to occur.