submodular.h

#ifndef SUBMODULAR_H
#define SUBMODULAR_H

#include<scil/scil.h>
#include<queue>

namespace SCIL {

class submodular : public sym_constraint {

   public:
      var v;
      var_map<int>& VM;
      int n;
      solution TS;
      double oldval;

   public:
      submodular(var& v_, var_map<int>& VM_, int n_)
        : v(v_), VM(VM_), n(n_) {
        feps = 1.0e-4;
      }

      virtual void init(subproblem& S) {
        return;
      }

      status standard_separation(subproblem& S);

      status feasible(solution& Sol) {
        return (Sol.value(v) > value(Sol) - .0001) ? feasible_solution : infeasible_solution;
      }

      virtual double value(solution& Sol) = 0;
   private:
      virtual double value(int p) {
        TS.set_value(VM[p], 1.);
        return value(TS);
      }
      virtual double value() {
        TS.set_zero();
        oldval = 0.;
        return 0.;
      }
};

};
#endif

---