base = 2

a = vector(2^13)

compute(n) = {
	if (n<base,
		return (n),
		my (s=Set(), ln=#digits(n,base));
		for (l=1, ln-1,
			for (o=0, ln-l,
				s = setunion(s, Set((n\(base^o))%(base^l)));
			);
		);
		return (sum(i=if (s[1]==0, 2, 1), #s, a[s[i]]))
	)
}

for (n=1, #a, a[n] = compute(n); print (n " " a[n]))

quit