login
A275464
a(n) = n - A038802(n).
1
0, 0, 0, 3, 1, 1, 6, 2, 2, 9, 3, 10, 12, 5, 5, 15, 15, 7, 18, 8, 8, 21, 9, 21, 24, 11, 25, 27, 13, 13, 30, 30, 15, 33, 16, 16, 36, 35, 18, 39, 19, 40, 42, 21, 42, 45, 45, 24, 48, 25, 25, 51, 26, 26, 54, 27, 55, 57, 56, 56, 60, 60, 33, 63, 34, 63, 66, 36, 36, 69, 67, 70, 72, 40, 40, 75, 75, 42, 78, 77, 44, 81
OFFSET
1,4
LINKS
FORMULA
a(n) = n - (A049084(A020639(2n+1)) - 1).
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
Cf. A038802.
Sequence in context: A109446 A088441 A061857 * A067433 A256697 A133567
KEYWORD
nonn
AUTHOR
Raphael Ranna, Jul 28 2016
EXTENSIONS
Duplicate a(44)-a(45) removed by Andrew Howroyd, Feb 23 2018
STATUS
approved