#include "molecules.h"
#include <bits/stdc++.h>
using namespace std;
#define F first
#define S second
vector<int> find_subset(int l, int u, vector<int> w) {
vector <pair <int, int>> p;
for (int x : w) {
p.push_back(make_pair (x, (int)p.size()));
}
sort (p.begin(), p.end());
vector <long long> ps = {0};
for (pair <int, int> x : p) {
ps.push_back(ps.back() + x.F);
}
int s = 0, t = -1;
for (int i = 0; i < (int)ps.size(); i++) {
int in = lower_bound (ps.begin(), ps.end(), ps[i] - u) - ps.begin();
long long x = ps[i] - ps[in];
if (x >= l)
s = in, t = i - 1;
}
vector <int> ans;
for (int i = s; i <= t; i++) {
ans.push_back(p[i].S);
}
return ans;
}
