%I#29 2023年9月15日17:09:22
%S 1,2,4,6,8,12,16,30,36,24,32,60,64,48,7221012818025612044,96,
%电话5124201619290024010243602048231028838443222604096768,
%电话:57684081927201638448018001536327684620129610811529065536
%N a(N)=分解为乘积p(i)^e(i)的最小整数,使得乘积p。
%C A025487的排列。a(n)是A025487的成员m,因此A181819(m)=n.a(n。
%C如果n=乘积_i素数(e(i)),且e(i)弱递减,则a(n)=乘积_i素数(i)^e(i)。例如,90=prime(3)*prime(2)*price(2)*prime(1),那么a(90)=prime_Gus Wiseman_,2019年1月2日
%H Alois P.Heinz,n的表格,n=1..10000的a(n)</a>
%H Eric Weisstein的数学世界,<a href=“http://mathworld.wolfram.com/ConjugatePartition.html“>共轭分区</a>
%F如果A108951(n)=乘积p(i)^e(i),则a(n)=乘积A002110(e(i))。即,a(n)=A108951(A181819(A10895l(n)))。
%F a(A181819(n))=A046523(n).-[另见A124859]_Antti Karttunen,2018年12月10日
%F a(n)=A025487(A361808(n))。-_Pontus von Brömssen,2023年3月25日
%F a(n)=A108951(A122111(n))_Antti Karttunen,2023年9月15日
%e 24的标准因式分解是2^3*3^1。因此,p(e(i))=素数(3)*素数(1)(即A000040(3)*A000040。因为24是p(e(i))=10,a(10)=24的最小整数。
%pa:=n->(l->mul(ithprime(i)^l[i],i=1..nops(l)))(排序(映射(i->
%p数理论[pi](i[1])$i[2],ifactors(n)[2]),`>`)):
%p序列(a(n),n=1..70);#_Alois P.Heinz,2018年9月5日
%t使用[{s=Array[If[#==1,1,Times@@Map[Prime@Last@#&,FactorInteger@#]]&,2^16]},Array[First@FirstPosition[s,#]&,LengthWhile[Differences@Union@s,#==1&]]](*Michael De Vlieger_,2018年12月17日*)
%t表[Times@@MapIndexed[Prime[#2[[1]]^#1&,Reverse[Flatten[Cases[FactorInteger[n],{p_,k_}:>表[PrimePi[p],{k}]]],{n,30}](*_Gus Wiseman_,2019年1月2日*)
%o(PARI)A181821(n)={my(f=因子(n),p=0,m=1_Antti Karttunen,2018年12月10日
%o(Python)
%o来自math导入prod
%o来自sympy import prime,primepi,factorint
%o def A181821(n):返回prod(prime(i)**e for i,e in enumerate(sorted(map(primepi,factorint(n,multiple=True)),reverse=True,1))#_Chai Wah Wu_,2023年9月15日
%Y A025487的其他重排包括A036035、A059901、A063008、A077569、A085988、A086141、A087443、A108951、A181822。
%Y参考A046523、A181819、A181820。
%Y参见A000040、A001221、A001225、A002110、A056239、A071625、A112798、A118914、A122111、A124859、A182850、A305936、A361808。
%K nonn公司
%O 1,2号机组
%A _马修·范德马斯特,2010年12月7日
%E由Gus Wiseman_修正的定义,2019年1月2日
|