%I#8 2014年10月19日15:41:05
%S 1,2,3,4,5,5,6,7,8,8,9,10,11,12,13,14,15,16,17,18,19,19,20,21,
%电话22,23,23,24,25,26,27,27,28,29,30,31,32,33,34,35,36,37,38,39,
%U 40,41,41,42,43,44,45,45,46,47,48,49,49,50,51,52,52,53
%N最小k,使4*k/v(2*k)^2-Pi<1/N,其中序列v在注释中定义。
%C序列v的定义如下:v(1)=0,v(2)=1,v(n)=v(n-1)/(n-2)+v(n-2。当n>=2时,a(n+1)-a(n)似乎位于{0,1}。
%D Steven R.Finch,《数学常数》,剑桥大学出版社,2003年,第19页。
%H Clark Kimberling,<a href=“/A247971/b247971.txt”>n,a(n)表,n=1..1000</a>
%e前几个项w(n)=4*n/v(2*n)^2-Pi和1/n的近似值:
%电子。。。4*n/v(2*n)^2-Pi。。。1个/个
%e 1。。。0.858407 .......... 1
%e 2。。。0.413963 ......... 0.5
%e 3。。。0.271741 .......... 0.333333
%e 4。。。0.202081 .......... 0.25
%e 5。。。0.160801 .......... 0.2
%e 6。。。0.133508 .......... 0.166666
%e a(2)=2,因为w(2)<1/2<w(1)。
%t$RecursionLimit=无限;z=400;v[1]=0;v[2]=1;
%tv[n]:=v[n]=v[n-1]/(n-2)+v[n-2];
%t表格形式[表格[{n,n[4 n/(v[2n]^2)-Pi],n[1/n]},{n,1,10}]]
%t f[n_]:=f[n]=选择[Range[z],4#/(v[2#]^2)-Pi<1/n&,1];
%t u=压扁[表[f[n],{n,1,z}]](*A247971*)
%t d=差异[u]
%t v=压扁[位置[d,0]](*A247972*)
%Y参见A247972、A247973和A247974。
%K nonn,简单
%O 1,2号机组
%百灵鸟金伯利2014年9月28日
|