login
A084812
a(1)=1; for n >= 2, choose a(n) so that |Sum_{k=1..n} 1/a(k)^z| < |Sum_{k=1..n-1} 1/a(k)^z|, where z = 1/2 + i*2*Pi.
7
1, 4, 5, 11, 14, 19, 29, 43, 59, 78, 99, 123, 150, 179, 212, 247, 286, 327, 372, 419, 470, 523, 580, 639, 702, 767, 836, 907, 982, 1059, 1140, 1223, 1310, 1399, 1492, 1587, 1686, 1787, 1892, 1999, 2110, 2223, 2340, 2459, 2582, 2707, 2836, 2967, 3102, 3239
OFFSET
1,2
PROG
(PARI) z=1/2 + 2*Pi*I; S=0; w=1; a=0; for(n=1, 100, b=a+1; while(abs(S+exp(-z*log(b)))>w, b++); S=S+exp(-z*log(b)); w=abs(S); a=b; print1(b, ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jun 04 2003
STATUS
approved