|
|
A375525
|
|
a(n) = numerator of Sum_{i=1..n} 1/A031216(i).
|
|
2
|
|
|
1, 5, 31, 247, 529, 559, 1747, 181, 4651, 147331, 1276726, 10379813, 10527373, 501426731, 507546131, 14874882499, 120130384067, 8121140473289, 597694372183697, 7829048769596861, 7886595152525441, 341425446875737163, 26430832817384951, 2579018441538210647, 137426518826740878691
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
EXAMPLE
|
The first few fractions are 1/2, 5/6, 31/30, 247/210, 529/420, 559/420, 1747/1260, 181/126, 4651/3150, 147331/97650, 1276726/830025, 10379813/6640200, 10527373/6640200, ...
|
|
MAPLE
|
b:= n-> (l-> add(l[i]*11^(i-1), i=1..nops(l)))(convert(ithprime(n), base, 10)):
g:= proc(n) option remember; `if`(n<1, 0, g(n-1)+1/b(n)) end:
a:= n-> numer(g(n)):
|
|
PROG
|
(PARI) a(n) = numerator(sum(i=1, n, 1/fromdigits(digits(prime(i)), 11))); \\ Michel Marcus, Aug 31 2024
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,frac,new
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|