%I#10 2021年6月23日09:58:59
%编号:6,23123341131293912939572278964332988110730831197553,
%电话:717107917107957301247208327509117178057712196841371478297171,
%电话:258757143327958359798663327331366640941110307711027332114885171036993025586797
%N a(1)=6,a(N)=a(N-1)的两个最接近的除数的串联,其乘积等于a(N-1),或者如果a(N-l)是素数,则1和a(N-1.)的串联。
%H Harvey P.Dale,n表,n=1..100的a(n)</a>
%t f[n_Integer]:=(d=除数[n];l=长度[d];If[EvenQ[l],ToExpression[ToString[d[l/2]]]<>ToString[d[[l/2+1]]],ToExpression[ToString[d[l/2+.5]]]<>ToString[d[l/2+/.5]]]);嵌套列表[f,6,25]
%t tcf[n_]:=模[{d=除数[n],len},len=长度[d]/2;FromDigits[Flatten[IntegerDigits/@Take[d,{len,len+1}]]];ctc[n_]:=如果[PrimeQ[n],10^整数长度[n]+n,tcf[n]];NestList[ctc,6,30](*哈维·P·戴尔,2019年5月19日*)
%o(Python)
%o来自symby导入除数,isprime
%o定义aupton(术语):
%o alst=[6]
%o表示范围(2,术语+1)中的n:
%o if isprime(alst[-1]):alst.append(int('1'+str(alst[-1]))
%o其他:
%o divs=除数(alst[-1])
%o d1=divs[(长度(divs)-1)//2]
%o d2=alst[-1]//d1
%o alst.append(int(字符串(d1)+字符串(d2)))
%o返回alst
%o打印(aupton(27))#_Michael S.Branicky_,2021年6月23日
%Y参考A003681、A062094和A062095。
%K基,nonn
%O 1,1号机组
%A _Robert G.Wilson v_,2001年8月8日
%E定义由Harvey P.Dale_于2019年5月19日澄清
|