\\ A343404
a(n) = {
	my (v=Mod(0,1));
	forprime (p=2, oo,
		if (n==0,
				return (lift(v)),
				v=chinese(v, Mod(n, p));
				n\=p;
		);
	);
}

vv = vector(1000)
n = 0
emit(v) = {
	print (n++ " " vv[n]=v);
	if (n==#vv,
		quit;
	);
}

{
	pp = 1;
	emit(0);
	forprime (p=2, oo,
		nn = n;
		for (d=1, p-1,
			for (m=1, nn,
				v = d*pp + vv[m];
				if (a(v)==v,
					emit(v);
				);
			);
		);
		pp*=p;
	);
}

quit