for L in s1:
  i=L[0][0]
  for (ii,j) in P:
   if ii!=i: continue
   for w in (0,1):
    t=L[2]+sm(i)+tau(i,j)-w
    e=L[3]+ff(i,j)+2*w+1
    mu=vm(L[5],j)
    r=L[6]+bb(i,j)+(2*j-1)*w-pi[j-1]
    att((i,L[2],L[3],L[4],L[5],L[6]),j,t,e,6,mu,r,(i,j))