OFFSET
1,1
COMMENTS
a(n) >= 5 except for n = 3. - Robert Israel, Sep 19 2016
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
EXAMPLE
n=100: p(100)=541, m = 3*541 + 1 = 1624 = 2*2*2*7*29 so a(100)=29.
MAPLE
seq(max(numtheory:-factorset(3*ithprime(i)+1)), i=1..100); # Robert Israel, Sep 19 2016
MATHEMATICA
ffi[x_] := Flatten[FactorInteger[x]] ma[x_] := Part[Reverse[ffi[x]], 2] Table[ma[3*Prime[w]+1], {w, 1, 100}]
Table[FactorInteger[3Prime[n]+1][[-1, 1]], {n, 70}] (* Harvey P. Dale, Apr 10 2012 *)
PROG
(PARI) a(n) = vecmax(factor(3*prime(n)+1)[, 1]); \\ Michel Marcus, Sep 19 2016
(Magma) [#f eq 0 select 1 else f[#f][1] where f is Factorization(3*NthPrime(n)+1): n in [1..90]]; // Vincenzo Librandi, Sep 20 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Sep 18 2003
STATUS
approved