OFFSET
1,1
COMMENTS
10's complement of A306555.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = 9 - A306555(n) for n >= 2.
EXAMPLE
3^3 == 7 == 1/13 (mod 10).
53^3 == 77 == 1/13 (mod 100).
653^3 == 77 == 1/13 (mod 1000).
4653^3 == 3077 == 1/13 (mod 10000).
...
...047823914653^3 = ...923076923077 = 1 + (...999999999999)*(12/13) = 1 - 12/13 = 1/13.
MAPLE
op([1, 3], padic:-rootp(13*x^3-1, 10, 100)); # Robert Israel, Mar 24 2019
PROG
(PARI) seq(n)={Vecrev(digits(lift(chinese( Mod((1/13 + O(5^n))^(1/3), 5^n), Mod((1/13 + O(2^n))^(1/3), 2^n)))), n)} \\ Following Andrew Howroyd's code for A319740.
CROSSREFS
10-adic cube root of p/q:
q=13: A306555 (p=-1), this sequence (p=1).
KEYWORD
nonn,base
AUTHOR
Jianing Song, Feb 23 2019
STATUS
approved