OFFSET
0,4
COMMENTS
arctanh(x) = x + x^3/3 + x^5/5 + x^7/7 + ...
LINKS
Wikipedia, p-adic number
FORMULA
a(n) = (Sum_{i=0..floor((n-3)/4)} 4^(2*i+1)/(2*i+1)) mod 2^n.
EXAMPLE
a(3) = 4^1 mod 2^3 = 4;
a(6) = 4^1 mod 2^6 = 4
a(7) = (4^1 + 4^3/3) mod 2^7 = 68;
a(10) = (4^1 + 4^3/3) mod 2^10 = 708;
a(11) = (4^1 + 4^3/3 + 4^5/5) mod 2^11 = 1732;
a(14) = (4^1 + 4^3/3 + 4^5/5) mod 2^14 = 7876;
a(15) = (4^1 + 4^3/3 + 4^5/5 + 4^7/7) mod 2^15 = 7876.
a(18) = (4^1 + 4^3/3 + 4^5/5 + 4^7/7) mod 2^18 = 237252.
PROG
(PARI) a(n) = lift(sum(i=0, (n-3)/4, Mod(4^(2*i+1)/(2*i+1), 2^n)))
CROSSREFS
KEYWORD
nonn
AUTHOR
Jianing Song, Aug 16 2019
STATUS
approved