#include "molecules.h" #include <vector> #include <algorithm> using namespace std; vector<int> W; int n; vector<int> find_subset(int l, int u, vector<int> w) { W = w; n = w.size(); int i, x; vector<int> I(n); for(i = 0; i < n; i++) I[i] = i; sort(I.begin(), I.end(), [] (int x, int y) { return W[x] < W[y]; }); int sum = w[I[0]]; vector<int> empty_res; vector<int> res; if(sum > u) return empty_res; for(x = 1; x < n && sum + w[I[x]] < l; x++) { i = I[x]; sum += w[i]; res.push_back(i); } if(x == n-1) return empty_res; res.push_back(I[n-1]); return res; }