a = vector(10 000)	\\ A229037(n)-1
xx = vector(#a)		\\ forbidden values (as a bit mask)

{
	for (n=1, #a,
		a[n] = valuation(xx[n]+1, 2);
		
		for (k=1, n-1,
			if (n+k>#xx,
				break;
			);
			w = 2*a[n]-a[n-k];
			if (w>=0,
				xx[n+k] = bitor(xx[n+k], 2^w);
			);
		);
		print (n " " hammingweight(xx[n]));
	);
}

quit