%I#27 2021年5月20日08:19:55
%S 0,0,0,1,0,2,0,3,0,10,5,2,0,1,3,4,9,0,1,11,0,5,1,0,10.8,
%T 15,2,1,0,17,10,3,0,4,0,9,2,21,0,10,4,14,11,0,6,5,16,27,0,2,9,4,
%U 0,8,8,0,15,20,6,0,10,35,2,17,4,10,0,3,0,39,0,4,12,41,26,9,0,2,6,21,28,45,14,0,5,8,3,0,14,0,14,0,11
%N的连续不同素因子之间距离的乘积;如果n只有1个不同的素因子,则为零。
%C对于n≤41,a(n)=A049087(n)。
%H Antti Karttunen,n的表,n=1..65537的a(n)</a>
%t f[n_]:=模[{ps},如果[n<=1,0,ps=转置[FactorInteger[n]][[1];时间@@差异[ps]]];表[f[n],{n,100}](*T.D.Noe_,2012年8月20日*)
%t数组[Apply[Times,Differences@FactorInteger[#][[All,1]]/。{}->0]&,105](*Michael De Vlieger_,2018年9月10日*)
%o(Python)
%o来自sympy import primerage
%o素数=列表(素数范围(2500))
%o表示范围(1100)内的n:
%o d=n
%o上一个=0
%o产品=1
%o表示质数中的p:
%o如果d%p==0:
%o如果之前:
%o产品*=p-prev
%o当d%p==0时:
%o d//=p
%o如果d==1:
%o中断
%o上一个=p
%o如果prev==0:
%o产品=0
%o打印(产品,结束=',')
%o(PARI)A178921(n)=如果(1>=ω(n),0,my(ps=系数(n)[,1],m=1);对于(i=2,#ps,m*=(ps[i]-ps[i-1]));(m) );\\_Antti Karttune_,2018年9月7日
%Y参考A081060、A049087。
%Y另请参阅A137795。
%K nonn,简单
%O 1,10号
%A _Alex Ratushnyak,2012年8月18日
%E 2018年9月7日,来自安蒂·卡图宁的更多条款
|