prev = 1	\\ product of previous terms

{
	for (n=1, 25,
		f = factor((1+prev) * if (type(prev)==type(42), I, 1));
		v = 0;
		for (i=1, #f~,
			p = f[i,1];
			if (norm(p) > 1,
				if (v==0 || norm(v) > norm(p) || (norm(v) == norm(p) && imag(v) > imag(p)),
					v = p;
				);
			);
		);
		prev *= v;
		print (n " " norm(v));
	);
}

quit