登录
多重数为n的素数指标的多集的多集划分数。
56

%I#10 2018年12月12日07:15:42

%S 1,1,2,2,3,4,5,5,9,7,7,11,11,12,16,15,15,26,22,21,29,19,30,36,31,30,

%电话66,38,42,52,56,52,47,45,57,92,77,67,77,74101,98135,64137,97176,

%电话:135109109118105231249,97141181139297198385195269

%N多重数为N的素数指标的多集的多集划分数。

%H Andrew Howroyd,n的表,n的a(n)=1..700</a>

%F a(n)=A001055(A181821(n))。

%F a(素数(n)^k)=A219727(n,k)-_安德鲁·霍罗伊,2018年12月10日

%e{1,1,2,3}的a(12)=11个多集分区:

%电子{{1,1,2,3}}

%e{{1},{1,2,3}}

%e{{2},{1,1,3}}

%e{{3},{1,1,2}}

%电子{{1,1},{2,3}}

%电子{{1,2},{1,3}}

%e{{1},{1},{2,3}}

%e{{1},{2},}1,3}}

%e{{1},{3},}

%e{{2},{3},{1,1}}

%e{{1},{1}、{2}、{3}}

%t nrmptn[n_]:=联接@@MapIndexed[表[#2[[1]],{#1}]&,如果[n==1,{},展平[Cases[FactorInteger[n]//反转,{p_,k_}:>表[PrimePi[p],{k}]]];

%t facs[n_]:=如果[n<=1,{{}},连接@@表[Map[Prepend[#,d]&,选择[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}];

%t表[Length[facs[Times@@Prime/@nrmptn[n]]],{n,60}]

%o(PARI)

%o置换计数(v)={my(m=1,s=0,k=0,t);对于(i=1,#v,t=v[i];k=如果(i>1&&t=v[i-1],k+1,1);m*=t*k;s+=t);s!/m}

%o sig(n)={my(f=因子(n));concat(向量(f~,i,向量(f[i,2],j,素数))}

%o计数(sig)={my(n=vecsum(sig

%o a(n)={if(n==1,1,my(s=sig(n));if(#s==1

%Y参见A001055、A001970、A007716、A181821、A219727、A255906、A269134、A305936。

%Y参见A317532、A318283、A318285、A318286、A318360、A318361。

%K nonn公司

%氧1,3

%A _Gus Wiseman_,2018年8月23日