a000043 = [ 2, 3, 5, 7, 13, 17, 19, 31, 61, 89, 107, 127, 521, 607, 1279, 2203, 2281, 3217, 4253, 4423, 9689, 9941, 11213, 19937, 21701, 23209, 44497, 86243, 110503, 132049, 216091, 756839, 859433, 1257787, 1398269, 2976221, 3021377, 6972593, 13466917, 20996011, 24036583, 25964951, 30402457, 32582657, 37156667, 42643801, 43112609 ]

explore(sig, tau, i, v) = {
	if (tau && sig%tau==0 && sig/tau<=#a,
		my (n=sig/tau);
		if (n>1,
			a[n]=concat(a[n], [v]);
		);
	);
	for (j=i, oo,
		if (j>#a000043,
			quit
		);

		my (xsig = sig + a000043[j], xtau = tau + 1);
		if (xsig <= #a * xtau,
			explore(xsig, xtau, j+1, v*(2^a000043[j]-1)),
			break
		);
	);
}

a = vector(500, n, [])
a[1] = [1]

explore(0, 0, 1, 1)

{
	m=0;
	for (n=1, #a,
		apply (v -> print (m++ " " v), Set(a[n]));
	);
}

quit