OFFSET
0,1
COMMENTS
Also x^2 = A091661.
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..9999
EXAMPLE
3^5 - 3 == 0 mod 10,
93^5 - 93 == 0 mod 10^2,
193^5 - 193 == 0 mod 10^3,
4193^5 - 4193 == 0 mod 10^4.
From Seiichi Manyama, Aug 01 2019: (Start)
8^(5^0) + 5^(2^0) == 3 mod 10,
8^(5^1) + 5^(2^1) == 93 mod 10^2,
8^(5^2) + 5^(2^2) == 193 mod 10^3,
8^(5^3) + 5^(2^3) == 4193 mod 10^4. (End)
PROG
(Ruby)
def P(n)
s1, s2 = 2, 8
n.times{|i|
m = 10 ** (i + 1)
(0..9).each{|j|
k1, k2 = j * m + s1, (9 - j) * m + s2
if (k1 ** 5 - k1) % (m * 10) == 0 && (k2 ** 5 - k2) % (m * 10) == 0
s1, s2 = k1, k2
break
end
}
}
s2
end
def Q(s, n)
n.times{|i|
m = 10 ** (i + 1)
(0..9).each{|j|
k = j * m + s
if (k ** 2 - k) % (m * 10) == 0
s = k
break
end
}
}
s
end
def A290375(n)
str = (P(n) + Q(5, n)).to_s.reverse
(0..n).map{|i| str[i].to_i}
end
p A290375(100)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Seiichi Manyama, Jul 28 2017
STATUS
approved