%I#23 2024年2月19日16:54:15
%S 2,216765109462178309149303521023341557014087331134903170,
%电话480752697620365011074329512800992258514337171548008755920,
%电话:106102098577237234602481117669030460994498455401187926421114850779780503464622906707
%N Fibonacci数等于promic数或promic数词的变位。
%C有许多术语与A162394重叠。前20个术语中的13个在A162394中_Chai Wah Wu_,2016年12月26日
%H Michael S.Branicky,n的表,n=1..28的a(n)</a>
%e a(1)=A000045(3)=2是促销号A002378(1)。
%e a(2)=A000045(8)=21是12=A002378(3)的变位符。
%e a(3)=A000045(20)=6765是7656=A002378(87)的变位。
%e a(4)=A000045(21)=10946是19460=A002378(139)和96410=A00237(310)的变位。
%e a(5)=A000045(32)=2178309是7308912=A002378(2703)、8017392=A00237.8(2831)和9731280=A0023.78(3119)的变位符。
%p位数:=200;digFreq:=处理(n)dgs:=转换(n,基数,10);c:=[序列(0,j=0..9)];对于dgs中的d,do c:=底土(d+1=op(d+1,c)+1,c);od:c;结束时间:
%p是Ana:=进程(a,b)返回(digFreq(a)=digFreg(b));结束时间:
%p为A002378:=进程(n)k:=楼层(sqrt(n));返回(k*(k+1)=n);结束时间:
%p A055642:=过程(n)最大值(1,ilog10(n)+1);结束时间:
%p代表n从1dof:=组合[fibonacci](n);fdgs:=A055642(f);对于来自楼层的k(sqrt(10^(fdgs-1)),do p:=k*(k+1);如果A055642(p)>fdgs,则断裂;fi;如果是Ana(p,f),则打印f(“%d,\n”,f);断裂;fi;od;od:#R.J.Mathar_,2009年8月14日
%o(Python)
%o从数学导入isqrt
%o从itertools导入计数,islice
%o def promics(d):#具有d位的promic数字的生成器
%o ilb,lb,ub=isqrt(10**(d-1)-1),10**(d-1),10**d
%o如果d==1:产量为0
%o对于计数中的i(ilb):
%o p=i*(i+1)
%o如果p<lb:继续
%o elif p>=ub:中断
%o产量p
%o def agen():序列项的#生成器
%o#yield 0#如果插入0,则取消注释
%o f,g,n,adict=1,2,1,dict()
%计数(1)中的d为o:
%o P=promics(d)中P的集合(“”.join(sorted(str(P)))
%o当f<10**d时:
%o如果P:yield f中的“”.join(sorted(str(f)))
%o f,g=g,f+g
%o打印(列表(islice(agen),16))#_Michael S.Branicky_,2024年2月19日
%Y参考A000045、A002378、A162394。
%K nonn,基础
%O 1,1号机组
%A _克劳迪奥·梅勒,2009年7月2日
%E 2009年8月14日,R.J.Mathar_又有六个学期
%E a(13)-a(18),来自Donovan Johnson_,2009年10月11日
%E a(19)-a(20),来自Chai Wah Wu,2017年1月2日
|