登录
(以n为基数2)的最长回文子字的长度。
14

%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日延期