|
Search: id:A134411
|
|
|
| A134411 |
|
a(n) is the smallest positive integer such that the numerator of (sum{k=1 to n} 1/a(k)) is prime (or 1), for all positive integers n. |
|
+0
5
|
|
| 1, 1, 1, 2, 3, 1, 3, 1, 1, 2, 6, 1, 1, 2, 4, 1, 1, 1, 3, 1, 1, 2, 3, 1, 2, 3, 2, 1, 4, 6, 1, 3, 1, 3, 2, 8, 3, 2, 3, 3, 1, 2, 6, 2, 1, 3, 3, 1, 5, 4, 3, 2, 1, 3, 1, 4, 2, 1, 3, 2, 1, 3, 1, 1, 3, 1, 1, 2, 6, 2, 4, 5, 3, 1, 3, 2, 1, 3, 3, 1, 3
(list; graph; listen)
|
|
|
OFFSET
|
1,4
|
|
|
LINKS
|
Leroy Quet, Home Page (listed in lieu of email address)
|
|
EXAMPLE
|
The sum of the reciprocals of the first 9 terms is 1 + 1 + 1 + 1/2 + 1/3 + 1 + 1/3 + 1 + 1 = 43/6. (And the numerator, 43, is prime.) Adding the reciprocal of 1 to this gets 49/6 (in reduced form). But 49 is composite. However, adding the reciprocal of 2 to 43/6 gets 23/3 (when written in reduced form). 23 is a prime, so therefore a(10) = 2.
|
|
MATHEMATICA
|
a = {1}; s = 1; Do[i = 1; While[ ! PrimeQ[Numerator[s + 1/i]], i++ ]; s = s + 1/i; AppendTo[a, i], {80}]; a - Stefan Steinerberger (stefan.steinerberger(AT)gmail.com), Oct 27 2007
|
|
CROSSREFS
|
Cf. A134412, A134413.
Sequence in context: A118846 A082503 A064442 this_sequence A126044 A114899 A023678
Adjacent sequences: A134408 A134409 A134410 this_sequence A134412 A134413 A134414
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
Leroy Quet Oct 24 2007
|
|
EXTENSIONS
|
More terms from Stefan Steinerberger (stefan.steinerberger(AT)gmail.com), Oct 27 2007
|
|
|
Search completed in 0.002 seconds
|
