|
%I #17 May 14 2022 11:20:43
%S 2,5,5,17,17,37,37,37,41,101,101,101,101,197,197,257,257,257,257,401,
%T 401,401,401,577,577,677,677,677,677,677,677,677,677,677,677,1297,
%U 1297,1297,1297,1601,1601,1601,1601,1601,1601,1601,1601,1601,1601,1601,1601
%N Greatest prime factor of Product_{x=1..n} (x^2 + 1).
%H C. Hooley, <a href="http://link.springer.com/content/pdf/10.1007%2FBF02395047">On the greatest prime factor of a quadratic polynomial</a>, Acta Mathematica July 1967, Volume 117, Issue 1, pp 281-299.
%p a:= proc(n) option remember; `if`(n=0, 0,
%p max(a(n-1), numtheory[factorset](n^2+1)[]))
%p end:
%p seq(a(n), n=1..55); # _Alois P. Heinz_, Jan 03 2021
%t a[n_] := a[n] = If[n == 1, 2, Max[a[n-1], FactorInteger[n^2+1][[-1, 1]]]];
%t Table[a[n], {n, 1, 55}] (* _Jean-François Alcover_, May 14 2022, after _Alois P. Heinz_ *)
%o (PARI) a(m) = {for (n=1, m, f = factor(prod(x=1, n, x^2+1)); print1(f[length(f~), 1], ", "););}
%Y Cf. A002496, A002522, A005574, A014442.
%K nonn
%O 1,1
%A _Michel Marcus_, Nov 23 2012
|
