OFFSET
1,1
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
MATHEMATICA
odd[n_] := n / 2^IntegerExponent[n, 2]; a[n_] := FactorInteger[odd[(Prime[n]+1)*(Prime[n+1]+1)]][[1, 1]]; Array[a, 100] (* Amiram Eldar, Apr 06 2025 *)
PROG
(Haskell)
a079083 = a078701 . a079079 -- Reinhard Zumkeller, Oct 08 2012
(PARI) odd(n) = n >> valuation(n, 2);
a(n) = factor(odd((prime(n)+1)*(prime(n+1)+1)))[1, 1]; \\ Amiram Eldar, Apr 06 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Dec 22 2002
STATUS
approved
