fusc(n)=local(a=1, b=0); while(n>0, if(bitand(n, 1), b+=a, a+=b); n>>=1); b \\ Charles R Greathouse IV, Oct 05 2008

{
	a = vector(2^13, n, -1);
	u = 1;
	for (k=0, oo,
		v = fusc(k) + fusc(k+1);
		if (v<=#a && a[v]<0,
			a[v] = k;
			while (a[u]>=0,
				print (u " " a[u]);
				if (u++ > #a,
					break (2);
				);
			);
		);
	);
}

quit