OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
Sum_{n>=1} 1/a(n) = Sum_{q in A019434} (1/q) * Product_{primes p <= q} p/(p-1) = 2.11507281179419933506... . If there are only 5 Fermat primes, then this sum is a rational number: 2.6569019...*10^23656 / 1.2561751...*10^23656.
MATHEMATICA
fermatQ[p_] := OddQ[p] && p == 2^IntegerExponent[p-1, 2] + 1;
Select[Range[300], fermatQ[FactorInteger[#][[-1, 1]]] &]
PROG
(PARI) isfermat(p) = p % 2 && p == 1 << valuation(p-1, 2) + 1;
isok(k) = k > 1 && isfermat(vecmax(factor(k)[, 1]));
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Jan 09 2026
STATUS
approved
