#include "molecules.h"
#include <bits/stdc++.h>
using namespace std;
typedef pair<int, int> ii;
std::vector<int> find_subset(int l, int u, std::vector<int> w) {
vector<ii> wsorted;
for (unsigned int i = 0; i < w.size(); i++) wsorted.push_back({w[i], i});
sort(wsorted.begin(), wsorted.end(), greater<ii>());
int totWeight = 0;
vector<int> res;
for (auto el : wsorted) {
int weight, i;
tie(weight, i) = el;
if (totWeight + weight <= u) {
totWeight += weight;
res.push_back(i);
}
if (totWeight >= l) return res;
}
return std::vector<int>(0);
}
