OFFSET
0,4
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..10000
FORMULA
Define the sequence {b(n)} by the recurrence b(0) = 0 and b(1) = 1, b(n) = b(n-1) + 3 * b(n-1)^3 + 7 mod 10^n for n > 1, then a(n) = (b(n+1) - b(n))/10^n.
EXAMPLE
1^3 == 1 (mod 10).
11^3 == 31 (mod 10^2).
11^3 == 331 (mod 10^3).
4011^3 == 3331 (mod 10^4).
64011^3 == 33331 (mod 10^5).
764011^3 == 333331 (mod 10^6).
MAPLE
op([1, 3], padic:-rootp(x^3+7/3, 10, 100)); # Robert Israel, Aug 09 2019
PROG
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Aug 09 2019
STATUS
approved