#include "molecules.h"
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
vector<int>find_subset(int l, int u, vector<int>w){
vector<int>p(w.size());
iota(p.begin(), p.end(), 0);
sort(p.begin(), p.end(), [&] (int i, int j) -> bool{
return w[i] < w[j];
});
int index = int(w.size()) - 1;
ll sum = 0;
while(index > -1 && sum < l){
sum += w[p[index--]];
}
int length = int(w.size()) - index - 1;
while(index > -1 && sum > u){
sum -= w[p[index + length]] - w[p[index]];
index--;
}
if(sum > u || sum < l){
return vector<int>{};
}
vector<int>ans(length);
for(int i = 1; i <= length; i++){
ans[i - 1] = p[index + i];
}
return ans;
}
