%I#57 2021年9月22日05:57:52
%S 1,2,2,4,4,3,2,6,9,5,11,4,3,11,6,13,9,11,7,3,20,14,5,2,9,23,20,12,
%电话:4,31,17,5,23,12,32,17,22,32,15,26,14,42,2,11,37,29,46,27,14,9,48,6,
%U 40,2,43,22,51,18,12,43,17,39,56,14,32,45,6,50
%N设S_N是以二进制形式写入的前N个素数的列表,左边是最低有效位,并按字典顺序排序;a(n)=第n素数在S_n中的位置。
%C a(A264647(n))=n和a(m)!=n表示m<A264647(n)。-_Reinhard Zumkeller_,2015年11月19日
%C A264662(n,a(n))=A000040(n):a(n_Reinhard Zumkeller_,2015年11月20日
%H John Bodeen,<a href=“/A263856/b263856.txt”>n,a(n)表,n=1..24771</a>
%H John Bodeen,<a href=“https://github.com/jbodeen/ava/blob/master/camlib/primes.ml“>OCaml程序生成序列</a>
%e S_1=[01],a(1)=1;
%e S_2=[01,11],a(2)=2;
%e S_3=[01,101,11],a(3)=2;
%e S_4=[01,101,11,111],a(4)=4;
%e S_5=[01,101,11,1101,111],a(5)=4;
%e S_5=[0101110111111111],a(6)=3;
%e。。。
%p s:=proc(n)s(n):=cat(“”,转换(ithprime(n),base,2)[])结束:
%pa:=n->ListTools[BinarySearch](排序([seq(s(i),i=1..n)]),s(n)):
%p序列(a(n),n=1..100);#_阿洛伊斯·海因茨(Alois P.Heinz),2015年11月19日
%t S[n_]:=S[n]=SortBy[Prime[Range[n]],StringJoin@@ToString/@Reverse[IntegerDigits[#,2]]&];
%t a[n_]:=第一位置[S[n],素数[n]][1];
%t表[a[n],{n,1100}](*Jean-François Alcover_,2021年9月22日*)
%o(哈斯克尔)
%o导入数据。列表(insertBy);导入数据。功能(打开)
%o导入数据。列表(elemIndex);导入数据。也许(来自Just)
%o a263856 n=a263856列表!!(n-1)
%o a263856_list=f[]a004676_list,其中
%o f bps(x:xs)=y:f bps'xs其中
%o y=fromJust(elemIndex x bps')+1
%o bps'=插入依据(compare `on`(reverse.show))x bps
%o--_Reinhard Zumkeller_,2015年11月19日
%o(Python)
%o来自sympy import prime
%o定义A263856(n):
%o返回1+排序(范围(1,n+1)中i的格式(prime(i),'b')[::-1]
%Y A004676是字典排序所依据的序列。
%Y参考A264596。
%Y参见A264647和A264662。
%K nonn,基础
%O 1,2号机组
%A _John Bodeen_,2015年10月28日
|