login
A309847
Digits of the multiplicative inverse of A309752.
1
1, 0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 0, 1, 0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 1, 1, 1, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0
OFFSET
-1
FORMULA
By definition, (Sum_{i=-1..n} a(i)*2^(i+1)) * (Sum_{i=0..floor((n+1)/2)} (-1)^i*2^(2*i)/(2*i+1)) == 1 (mod 2^(n+2)).
EXAMPLE
arctan(2) = ...0101100110111010110101111110001011001010, so 1/arctan(2) = ...100010111110000010011001110110000110110.1.
PROG
(PARI) a(n) = lift(Mod(sum(i=0, (n+1)/2, (-1)^i*2^(2*i)/(2*i+1)), 2^(n+2))^(-1))\2^(n+1)
CROSSREFS
Sequence in context: A267605 A373550 A319843 * A266786 A272532 A166946
KEYWORD
nonn,base
AUTHOR
Jianing Song, Aug 20 2019
STATUS
approved