|
|
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
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
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
|
|
|
STATUS
|
approved
|
|
|
|