subject to c1a {k in K, o in O, n in N: n!=t[k,o]}:
(f[n,o]==1) ==>
sum{e in E}
(a[n,e]*x[e,k,o]) -
sum{e in E}
(b[n,e]*x[e,k,o]) =
(if r[n,k]==1
then d[k,o]*(1-f[k,o])
else 0);
subject to c1b {k in K, o in O, n in N: n!=t[k,o]}:
(f[n,o]==0) ==>
sum{e in E}
(a[n,e]*x[e,k,o]) -
sum{e in E}
(b[n,e]*x[e,k,o]) = 0;
The ==> operator is used to create an indicator constraint (if expression-involving-variables is true then constraint holds else constraint need not hold).
Paul