a(n, base=2) = {
	my (d=if (n, digits(n, base), [0]), nb=vector(#d));
	for (i=1, #d,
		for (j=0, i-1,
			my (m=fromdigits(Vecrev(d[j+1..i]), base),
			    x=base^(i-j));
			nb[i]=if (j==0, [m], setunion(nb[i], Set(apply(p -> p*x+m, nb[j]))));
		);
	);
	nb[#d][1];
}

for (n=0, 2^14, print (n " " a(n)))

quit