implement DFJ constraints in algorithms.lpsolve.tsp()