%I#40 2013年10月31日12:17:17
%S 1,1,2,2,3,2,3,3,4,3,3,1,3,2,4,4,5,4,4,1,4,2,3,1,5,3,5,5,6,5,
%温度:5,5,4,4,4,3,4,5,5,4,6,5,4,4,5,5,4,1,6,6,5,
%U 5,5,5,1,7,5,4,6,4,5,4,4,4,5,5,6,4,5,7,5,4,4,6,6,6,6,5,7,5,5,6,6
%N(N基2)的最长回文子字的长度。
%Ca(A083318(n-1))=n;a(A193159(k))=3,1<=k<=26。[_Reinhard Zumkeller_,2011年7月17日]
%H Reinhard Zumkeller,n的表,n的a(n)=1.16385=2^14+1</a>
%F a(n)<=最小值(a(2*n),a(2*n+1))。[_Reinhard Zumkeller_,2011年7月31日]
%e(11基2)=1011,包含回文101,因此a(11)=3。
%p#A050430[M1..M2]中n的n的最长回文因子的长度-摘自n.J.A.Sloane_,2012年8月7日,于2012年8月月11日修订
%p isPal:=进程(L)
%p局部d;
%p代表d从1到nops(L)/2 do
%p如果op(d,L)<>op(-d,L),则
%p返回false;
%p结束if;
%p端do:
%p返回true;
%p结束过程:
%p#主程序启动
%答案:=[];
%p M1:=0;M2:=64;
%从M1到M2的n的p do
%p t1:=换算(n,基数,2);
%p记录:=0:
%pl1:=nops(t1);
%j1的p从0到l1-1 do
%p代表j2从j1+1到l1 do
%p F1:=[op(j1+1..j2,t1)];
%p如果(isPal(F1)和j2-j1>rec),则rec:=j2-j1;fi;
%日期:
%日期:
%p ans:=[操作(ans),记录]:
%日期:
%p an;
%t f[n_]:=块[{id=整数位数[n,2]},k=长度@id;而[Union[#==Reverse@#&/@Partition[id,k,1]][[-1]]!=真,k--];k];数组[f,105](*_Robert G.Wilson v_,2011年7月16日*)
%o(哈斯克尔)
%o导入数据。字符(intToDigit,digitToInt)
%o导入数字(showIntAtBase)
%o a050430 n=a050430_列表!!(n-1)
%o a050430_list=f 1,其中
%o f n=g(showIntAtBase 2 intToDigit n“”):f(n+1)
%o g zs | zs==反向zs=长度zs
%o |否则=最大值(h$init zs)(h$tail zs
%o h zs@('0':_)=g zs
%o h zs@('1':_)=a050430$foldl(\v d->digitalToInt d+2*v)0 zs
%o--_Reinhard Zumkeller,2011年7月16日
%Y参考A007088;A050431(基数3)、A050432(基数4)、A050433(基数5)。
%Y参见A215244、A215467、A215256。
%K非n,碱基
%氧1,3
%百灵鸟金伯利_
%E由雷·钱德勒于2010年3月11日延期