{
	a = vector(10 000);
	u = 1;
	f = [0];

	for (m=1, oo,
		r = m;
		w = 0;	\\ order of sequence
		s = 0;	\\ sum of previous w terms
		for (k=1, oo,
			if (k>#f,
				f = concat(f, vector(#f));
			);

			if (r,
				w++;
				s += f[k] = valuation(r+r%2, 2);
				r \= 2^f[k],

				f[k] = s;
				s = s + f[k]-f[k-w];
			);

			if (f[k]>#a,
					break,
				a[f[k]]==0,
					a[f[k]] = m
			);

			if (r==0 && w==1,
				break;	\\ order 1 = constant
			);
		);

		while (a[u],
			print (u " " a[u]);
			if (u++>#a,
				break (2);
			);
		);
	);
}

quit