|
|
A332131
|
|
a(n) = (10^(2n+1)-1)/3 - 2*10^n.
|
|
4
|
|
|
1, 313, 33133, 3331333, 333313333, 33333133333, 3333331333333, 333333313333333, 33333333133333333, 3333333331333333333, 333333333313333333333, 33333333333133333333333, 3333333333331333333333333, 333333333333313333333333333, 33333333333333133333333333333, 3333333333333331333333333333333
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
COMMENTS
|
See A183174 = {1, 3, 7, 61, 90, 92, 269, ...} for the indices of primes.
|
|
LINKS
|
|
|
FORMULA
|
G.f.: (1 + 202*x - 500*x^2)/((1 - x)(1 - 10*x)(1 - 100*x)).
a(n) = 111*a(n-1) - 1110*a(n-2) + 1000*a(n-3) for n > 2.
|
|
MAPLE
|
A332131 := n -> (10^(2*n+1)-1)/3-2*10^n;
|
|
MATHEMATICA
|
Array[3 (10^(2 # + 1)-1)/9 - 2*10^# &, 15, 0]
|
|
PROG
|
(PARI) apply( {A332131(n)=10^(n*2+1)\3-2*10^n}, [0..15])
(Python) def A332131(n): return 10**(n*2+1)//3-2*10**n
|
|
CROSSREFS
|
Cf. A138148 (cyclops numbers with binary digits), A002113 (palindromes).
Cf. A332121 .. A332191 (variants with different repeated digit 2, ..., 9).
Cf. A332130 .. A332139 (variants with different middle digit 0, ..., 9).
|
|
KEYWORD
|
nonn,base,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|