%I#25 2022年4月25日08:04:48
%S 1,1,2,1,3,1,2,3,5,1,1,7,5,2,1,3,5,7,11,1,5,13,3,7,1,1,1,2,11,
%T 17,7,1,1,19,13,1,7,11,11,1,23,1,7,5,17,13,13,11,19,29,1,1,1,31,
%U 1,2,13,11,1,17,23,1,1,1,37,5,19,11,13,1,1,3,3,41,1,17,43,29,11,1,1,13
%N a(N)=rad(N!)/rad(A001142(N)),其中rad(N)是N的无平方核,A007947(N)。
%C与卢克·卢梭观察到的数据不匹配的前一个名称是:A002944(n)和A001405的无平方核的商。
%C a(n)是A001142(n)的素因式分解中缺失的不大于n的唯一素数p,如果存在这样的素数;否则,a(n)为1_卢克·卢梭,2019年1月1日
%H Luc Rousseau,n表,n=1..1000的a(n)(来自Labos Elemer的前90个术语)
%F a(n)=A034386(n)/A056606(n).-_Sean A.Irvine_,2022年4月24日
%e来自卢克·卢梭,2019年1月2日:(开始)
%e在帕斯卡三角形中,
%e行n=3(1 3 3 1)不包含素数因子为2的数,因此a(3)=2;
%e-行n=4(1 4 6 4 1)包含,对于所有p素数<=4,p的倍数,因此a(4)=1;
%e-行n=5(1 5 10 10 5 1)不包含素数因子为3的数字,因此a(5)=3;
%e等。
%e(结束)
%t L[n_]:=表[二项式[n,k],{k,1,Floor[n/2]}]
%t c[n_]:=补码[Prime/@Range[PrimePi[n]],First/@FactorInteger[Times@@L[n]]]
%t a[n_]:=模[{x=c[n]},如果[x=={},1,第一个[x]]]
%t表[a[n],{n,1,100}]
%t(*卢梭,2019年1月1日*)
%o(PARI)rad(n)=factorback(factorint(n)[,1]);\\A007947号
%o b(n)=prod(m=1,n,二项式(n,m));\\A001142号
%o a(n)=拉德(n!)/拉德(b(n));\\_米歇尔·马库斯,2019年1月2日
%Y参见A002944、A001405、A003418、A002110、A001142、A034386、A056606。
%K nonn公司
%氧1,3
%A _Labos Elemer,2000年8月7日
%E定义和示例由_Luc Rousseau_更改,2019年1月2日
|