%I#14 2022年9月8日08:46:24
%S 1,1,1,1,2,1,1,1,1,1,1,36,1,1,1,1,1,1,1,2,1,1,11,1,36,1,1,4,1,10,1,1,1,
%T 1,1,36,1,1,1,1,1,2,1,1,1,1,1,1,36,1,11,1,1,2,1784,1,1180,1,1,11,1,
%U 1,2,1,1,1,1,1,1,36,1,1,1,1,1,2,1,1,1,1,1
%N a(N)是N的除数d的乘积,因此sigma(d)除以N。
%C a(A097603(n))>1。
%C有关此类除数的数量和总和,请参见A173441和A326697。
%H Antti Karttunen,n的表,n=1..20000的a(n)</a>
%e对于n=12,除数d为12:1,2,3,4,6,12;
%e对应的σ(d):1、3、4、7、12、28;
%eσ(d)将n除以4个除数d:1,2,3,6;
%e a(12)=1*2*3*6=36。
%t a[n_]:=倍@@Select[Divisors[n],Divisible[n,DivisorSigma[1,#]&]];阵列[a,100](*_Amiram Eldar_,2019年7月21日*)
%o(岩浆)[&*[除数(n)中的d:d | IsIntegral(n/SumOfDivisors(d))]:[1..100]]中的n
%o(PARI)a(n)=我的(p=1);fordiv(n,d,if(!(n%西格玛(d)),p*=d));p、 \\_米歇尔·马库斯,2019年7月19日
%Y参见A000203、A000961、A173441、A326697。
%K nonn公司
%O 1,6型
%A_Jaroslav Krizek,2019年7月19日
|