%I#27 2020年7月24日07:12:45
%S 15、18、21、26、28、30、31、32、44、45、50、52、56、57、62、64、68、75、76、85、86、91、92,
%电话93,98,99110111116117129133146147148153164175183188207,
%电话:2122152192362442452592612682752792843143163253323383413433563636363693387388
%N个巴西数字,正好有两个巴西表示。
%C这些数字可以称为2-巴西数字。
%C该序列的最小数为15,这也是A257521中最小的奇数巴西复合物,15=11111 _2=33_4。A329383中的数字15高度巴西化。
%C根据古尔马提格猜想,只有两个素数,即31和8191,它们都是梅森数,在两个不同的基数下是巴西的(A119598)。
%H维基百科,<a href=“http://en.wikipedia.org/wiki/Goormaghtigh_consustructure(维基百科网)“>古尔马提格猜想</a>
%e 18=2*9=22_8=3*6=33_5。
%e 26=2*13=2*111_3=222_3=22_12。
%e 31=11111 _2=111_5;
%e 8191=11111111111 _2=111_90。
%对短烯基:=过程(n,b)
%p局部r、q、coupleq:
%p如果n<b则
%p返回[1,n]
%p其他
%p r:=(n模块b):
%p q:=(n-r)/b
%p耦合q:=短烯基(q,b):
%p如果r=耦合器[2]且r>0,则
%p返回[1]+1,r]
%p其他
%p返回[0,0]
%p结束if
%p结束if
%p端程序;
%p bresil:=程序(n)
%p局部b,L,k,t:
%p k:=0:
%b的p从2到(n-2)do
%p t:=短烯基(n,b):
%p如果t[1]>0,则
%p k:=k+1
%p L[k]:=[b,t[1],t[2]:
%p结束,如果:
%p端do:
%p序列(L[i],i=1..k);
%p端程序;
%p nbbresil:=n->nops([bresil(n)]);
%p#巴西数字的2倍
%p表示n从1到100 do,如果nbbresil(n)=2,则打印(n,bresil,n))else fi;日期:
%t压扁@位置[#,2]和@表[Count[Range[2,n-2],_?(和[Length@#!=1,Length@Union@#==1]&@IntegerDigits[n,#]&)],{n,400}](*_Michael De Vlieger_,2017年7月18日*)
%Y参考A085104、A125134、A220570、A220571、A257521、A284758、A288783、A290016、A290017、A290018、A329383。
%K nonn,基础
%O 1,1号机组
%A _伯纳德·肖特,2017年7月17日
|