MILP adjustments

- Add parameters for valid bound on phi (max_congestion)
  and for the approximation rate (r)
- Print parameters to problem.nfo
- Fix bug in relaxed constraint (must be valid for *all* s,t pairs)
- Add relaxed constraint using u[i,j] instead of all possible
  (s,t) pairs