OFFSET
1,2
COMMENTS
a(n) is the least k such that 2*s(n) < s(n+k), where s(n) is the n-th partial sum of primes.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = A231153(n) + 1.
EXAMPLE
5 <= 2 + 3 < 5 + 7, so a(2) = 2.
19 + 23 <= 2 + 3 + 5 + 7 + 11 + 13 + 17 < 19 + 23 + 27, so a(7) = 3.
MAPLE
P:= [seq(ithprime(i), i=1..10^5)]:
S:= ListTools:-PartialSums(P):
f:= proc(n) local j;
ListTools:-BinaryPlace(S, 2*S[n]+1/2)+1-n;
end proc:
map(f, [$1..200]); # Robert Israel, Mar 15 2026
MATHEMATICA
test = 0; u = Table[test = test + Prime[n]; tmp = 0; NestWhile[# + 1 &, n, test >= (tmp = tmp + Prime[#+1]) &] - n, {n, 1, 80}]; u + 1 (* Peter J. C. Moses, Nov 02 2013 *)
PROG
(PARI) a(n) = my(s=vecsum(primes(n)), k=1); while (vecsum(primes(n+k)) - 2*s <= 0, k++); k; \\ Michel Marcus, Mar 18 2026
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Clark Kimberling, Nov 09 2013
EXTENSIONS
Corrected by Robert Israel, Mar 15 2026
STATUS
approved
