OFFSET
1,4
LINKS
Raphael Ranna, Table of n, a(n) for n = 1..1000
FORMULA
a(n) = n + O(n/log n). - Charles R Greathouse IV, Jul 29 2016
EXAMPLE
9 = (2^0)*(3^2), resulting in 1 zero. So, a(4) = 4 - 1 = 3.
MATHEMATICA
Table[f = FactorInteger[2 n + 1]; n - (PrimePi[f[[1, 1]]] - 1), {n, 100}]
PROG
(PARI) lpf(n)=factor(n)[1, 1]
a(n)=n - primepi(lpf(2*n+1)) + 1 \\ Charles R Greathouse IV, Jul 29 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Raphael Ranna, Jul 28 2016
EXTENSIONS
Duplicate a(44)-a(45) removed by Andrew Howroyd, Feb 23 2018
STATUS
approved