Add node in flow degree 1 constrain and solver for multiple hubs

[fhk]

This will mean that the formulation can partition the area based on cost or capacity and create non overlapping regions.

Capacity can be enforced by constraining the outflow at potential hubs.