|
|
A304287
|
|
Sequence gives the denominators, in increasing values from the second term on, of Egyptian fractions whose sum has the concatenation of these denominators as decimal part. Case a(1) = 5.
|
|
16
|
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
|
|
LINKS
|
|
|
EXAMPLE
|
We start from 5: 1/5 = 0.2. At the beginning we have 2 instead of 5 as first decimal digit but the second term will fix it.
In fact the next integer is 3 and 1/5 + 1/3 = 0.5333... and so on.
The sum is 0.5 3 524 414362 964900433080 ...
|
|
MAPLE
|
P:=proc(q) local a, b, d, n; a:=1/5; b:=ilog10(1/a)+1; d:=1/a; print(d);
for n from 1 to q do if trunc(evalf(a+1/n, 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+1/n; print(n); fi; od; end: P(10^20);
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,base
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|