

A213062


Minimal sum x(1) +...+ x(n) such that 1/x(1) +...+ 1/x(n) = 1, the x(i) being n distinct positive integers.


3



1, 0, 11, 24, 38, 50, 71, 87, 106, 127, 151, 185, 211, 249, 288, 325, 364, 406, 459, 508, 550, 613, 676, 728
(list;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

1,3


COMMENTS

The term a(2)=0 corresponds to the fact that 1 cannot be written as Egyptian fraction with 2 (distinct) terms.


LINKS

Table of n, a(n) for n=1..24.
Eric Weisstein's World of Mathematics, Egyptian Fraction
Wikipedia, Egyptian fraction
Index entries for sequences related to Egyptian fractions


EXAMPLE

a(3) = 11 = 2 + 3 + 6, because 1/2+1/3+1/6 is the only Egyptian fraction with 3 terms having 1 as sum.
a(4) = 24 = 2 + 4 + 6 + 12 is the smallest sum of denominators among the six 4term Egyptian fractions equal to 1.
a(5) = 38 = 3 + 4 + 5 + 6 + 20, least sum of denominators among 72 possible 5term Egyptian fractions equal to 1.
a(6) = 50 = 3 + 4 + 6 + 10 + 12 + 15, least sum of denominators among 2320 possible 6term Egyptian fractions equal to 1.
a(7) <= 71 = 3 + 5 + 20 + 6 + 10 + 12 + 15 (obtained from n=6 using 1/4 = 1/5 + 1/20).
a(8) <= 114 = 3 + 5 + 20 + 7 + 42 + 10 + 12 + 15 (obtained using 1/6 = 1/7 + 1/42).
a(9) <= 145 = 3 + 6 + 30 + 20 + 7 + 42 + 10 + 12 + 15 (obtained using 1/5 = 1/6 + 1/30).
a(10) <= 202 = 3 + 6 + 30 + 20 + 8 + 56 + 42 + 10 + 12 + 15 (obtained using 1/7 = 1/8 + 1/56).


PROG

(PARI) a(n, M=9e9, s=1, m=2)={ n==1 & return((numerator(s)==1 & 1 >= m*s  s==1)/s); sum( k=m, m+n1, 1/k ) < s & return; for(x=max(m, 1\s+1), n\s, n*(x+(n1)/2)>=M & break; (m=a(n1, Mx, s1/x, x+1)) & M=min(M, x+m)); M} /* For n>6, a good upper bound must be given as 2nd (optional) argument. Such a bound can be obtained using 1/x = 1/(x+1) + 1/x(x+1) in solutions for n1, cf. Examples. */


CROSSREFS

Cf. A030659.  Alois P. Heinz, Sep 21 2012
Sequence in context: A063307 A213991 A269100 * A157756 A061043 A098603
Adjacent sequences: A213059 A213060 A213061 * A213063 A213064 A213065


KEYWORD

nonn,more


AUTHOR

M. F. Hasler, Jun 03 2012


EXTENSIONS

a(11)a(24) from Robert Price, Aug 26 2012  Sep 21 2012


STATUS

approved



