{
	a1 = 3;
	a = vector(1 000, k, Mod(0,1));
	for (n=1, #a,
		my (v = if (n==1, a1, lift(a[n])));
		print (n " " v);
		for (k=n+1, #a,
			a[k] = chinese(a[k], Mod(v, k-n))
		)
	)
}