a240542[n_]:=总和[(-1)^(k+1)*天花板[(n+1)/k-(k/1)/2],{k,1,楼层[(Sqrt[8n+1]-1)/2]}]
a282198[b]:=模块[{centers={{1,1}},acc={1},k=2,cPrev=1,cCur,len},While[k<=b,cCur=a240542[k];如果[Last[acc]==cCur,AppendTo[acc,cCur],len=Length[acc];如果[First[Last[centres]]<len,AppendTo[centres,{len,k-1}]];acc={cCur};cPrev=cCur];k++];最后一个[Transpose[centers]]]
a282198[5000000](*数据*)