#include "molecules.h"
#include <bits/stdc++.h>
using namespace std;
// #define int long long
#define ff first
#define ss second
#define pb push_back
#define ppb pop_back
#define meta int tm = (tl + tr) / 2, x = i * 2 + 1, y = x + 1
const int N = 2e5 + 7;
const int TN = 4 * N;
// const int oo = 1e18;
const int mod = 1e9 + 7;
typedef pair<int, int> pii;
typedef vector<int> vi;
typedef vector<pii> vii;
std::vector<int> find_subset(int l, int u, std::vector<int> w) {
vii v;
vi ans;
for (int i = 0; i < w.size(); i++) {
v.pb({w[i], i});
}
sort (v.begin(), v.end());
int sum = v[0].ff, tl = 0, tr = 0;
while (1) {
if (l <= sum && sum <= u) {
for (int i = tl; i <= tr; i++) {
ans.pb (v[i].ss);
} return ans;
}
if (sum < l) {
tr++;
if (tr >= w.size()) return vector<int> ();
sum += v[tr].ff;
}
else {
sum -= v[tl].ff;
tl++;
if (tl >= w.size()) return vector<int> ();
}
// cout << tl << ' ' << tr << ' ' << sum << '\n';
}
}
// signed main() {
// int n, l, u;
// assert(3 == scanf("%d %d %d", &n, &l, &u));
// vector<int> w(n);
// for (int i = 0; i < n; i++)
// assert(1 == scanf("%d", &w[i]));
// vector<int> result = find_subset(l, u, w);
// cout << result.size() << '\n';
// for (int i = 0; i < result.size(); i++) {
// cout << result[i] << ' ';
// }
// // printf("%d\n", (int)result.size());
// // for (int i = 0; i < (int)result.size(); i++)
// // printf("%d%c", result[i], " \n"[i == (int)result.size() - 1]);
// }
