Though the particulars of this model are unique, fundamentally it is closest to the work of Hu and Liu (2004a) and Carenini et al.
this optimization is well studied and can be solved with a wide variety of techniques.