OFFSET
0,2
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) + 7 * (9 * b(n-1)^3 + 71) mod 10^n for n > 1, then a(n) = (b(n+1) - b(n))/10^n.
EXAMPLE
1^3 == 1 (mod 10).
61^3 == 81 (mod 10^2).
361^3 == 881 (mod 10^3).
1361^3 == 8881 (mod 10^4).
61361^3 == 88881 (mod 10^5).
761361^3 == 888881 (mod 10^6).
PROG
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Seiichi Manyama, Aug 11 2019
STATUS
approved