|
%I #14 Sep 08 2022 08:46:18
%S 1,2,5,7,17,13,37,43,19,73,101,61,29,157,197,211,257,29,307,181,401,
%T 421,463,53,577,601,677,73,757,421,67,37,331,151,1123,613,1297,137,67,
%U 1483,1601,547,1723,139,631,1013,109,103,461,1201,61,2551,541,919
%N The largest prime factor of (1+n^2)*(1+n^3).
%H Vincenzo Librandi, <a href="/A281794/b281794.txt">Table of n, a(n) for n = 0..1000</a>
%F a(n) = max( A014442(n), A081256(n)).
%F a(n) = A006530(A281661(n)).
%p A281794 := proc(n)
%p A006530((1+n^2)*(1+n^3)) ;
%p end proc:
%p seq(A281794(n),n=0..60) ;
%t Table[Max[Transpose[FactorInteger[(1 + n^2) (1 + n^3)]][[1]]], {n, 0, 60}] (* _Vincenzo Librandi_, Jun 03 2017 *)
%o (Magma) [#f eq 0 select 1 else f[ #f][1] where f is Factorization((1+n^2)*(1+n^3)): n in [0..60]]; // _Vincenzo Librandi_, Jun 03 2017
%o (PARI) a(n) = if (n==0, 1, my(f=factor((1+n^2)*(1+n^3))); vecmax(f[, 1])); \\ _Michel Marcus_, Jun 03 2017; corrected Jun 13 2022
%Y Cf. A006530, A014442, A081256, A281661.
%K nonn,easy
%O 0,2
%A _R. J. Mathar_, Jan 30 2017
|
