login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A320308
Numerators of the fractions a(0)/(a(1) - a(0)), a(1)/(a(2) - a(1)), a(2)/(a(3) - a(2)), ... such that the sum 1/a(0) + Sum_{n>=1} a(n-1)/(a(n) - a(n-1)) has the concatenation of these numerators as decimal part. Case a(0) = 13.
15
13, 247, 890970451, 540818819108243685091145, 77126435926826212355124730859827553231849402446824773070904238
OFFSET
0,1
COMMENTS
It appears that fractions of this kind exist only for a(0) equal to 3 (A320306), 10 (A320307), 13 (this sequence) and 38 (A320309).
Next term has 163 digits. - Giovanni Resta, Oct 11 2018
EXAMPLE
1/13 = 0.07692...
At the beginning instead of 13 we have 07 as first decimal digits. Adding the second term this is fixed.
1/13 + 13/(247 - 13) = 0.13247863...
1/13 + 13/(247 - 13) + 247/(890970451 - 247) = 0.13247890970451376...
The sum is 0.13 247 890970451 540818819108243685091145 ...
MAPLE
P:=proc(q, h) local a, b, d, n, t, x; x:=h+1; a:=1/h; b:=ilog10(h)+1;
d:=h; print(d); t:=1/a; for n from x to q do if
trunc(evalf(a+t/(n-t), 100)*10^(b+ilog10(n)+1))=d*10^(ilog10(n)+1)+n
then b:=b+ilog10(n)+1; d:=d*10^(ilog10(n)+1)+n; a:=a+t/(n-t); t:=n;
x:=n+1; print(n); fi; od; end: P(10^20, 13);
KEYWORD
nonn,base
AUTHOR
Paolo P. Lava, Oct 11 2018
EXTENSIONS
a(3)-a(5) from Giovanni Resta, Oct 11 2018
STATUS
approved