OFFSET
0,3
COMMENTS
a(n) is the unique solution to x^5 == 3 (mod 2^n) in the range [0, 2^n - 1].
LINKS
Wikipedia, p-adic number
FORMULA
For n > 0, a(n) = a(n-1) if a(n-1)^5 - 3 is divisible by 2^n, otherwise a(n-1) + 2^(n-1).
EXAMPLE
For n = 2, the unique solution to x^5 == 3 (mod 4) in the range [0, 3] is x = 3, so a(2) = 3.
a(2)^5 - 3 = 240 which is divisible by 8, so a(3) = a(2) = 3;
a(3)^5 - 3 = 240 which is divisible by 16, so a(4) = a(3) = 3;
a(4)^5 - 3 = 240 which is not divisible by 32, so a(5) = a(4) + 16 = 19;
a(5)^5 - 3 = 2476096 which is divisible by 64, so a(6) = a(5) = 19.
PROG
(PARI) a(n) = lift(sqrtn(3+O(2^n), 5))
CROSSREFS
For the digits of 3^(1/5), see A325896.
Approximations of p-adic fifth-power roots:
this sequence (2-adic, 3^(1/5));
A325893 (2-adic, 5^(1/5));
A325894 (2-adic, 7^(1/5));
A325895 (2-adic, 9^(1/5));
A322157 (5-adic, 7^(1/5));
A309450 (7-adic, 2^(1/5));
A309451 (7-adic, 3^(1/5));
A309452 (7-adic, 4^(1/5));
A309453 (7-adic, 5^(1/5));
A309454 (7-adic, 6^(1/5)).
KEYWORD
nonn
AUTHOR
Jianing Song, Sep 07 2019
STATUS
approved