base = 2

a = vector(2^14, n, b=digits(n,base); if (b==Vecrev(b), n, 0))

{
	print ("0 0");
	for (n=1, #a,
		b=digits(n, base);

		for (k=2, #b,
			t=b[k..#b];
			if (t==Vecrev(t),
				h=fromdigits(b[1..k-1],base);
				a[n]=max(a[n], a[h]+fromdigits(t,base));
			);
		);

		print (n " " a[n]);
	);
}

quit