%我#45 2024年5月12日11:35:12
%S 1,2,2,4,2,5,2,6,4,5,2,10,2,5,10,5,9,2,10,10,10,5,5,2,16,45,6,10,2,14,2,
%电话:12,5,5,50,20,2,5,5,16,2,14,2,10,10,5,2,24,4,10,10,2,16,5,2,2,
%U 28,2,5,10,10,16,5,14,2,10,5,14,2,32,2,5,10,10,5,4,2,24,9,5,2,28,5,5,16,28,5
%考虑N的所有(不一定是不同的)正除数对的乘积。
%C一个长方体中具有两种不同边长的n个全等长方体的3D网格数,模数旋转(对于立方体,请参见A034836,而不是长方体,对于具有三种不同边长度的长方体请参见A007425;对于2D情况,请参阅A000005)_Manfred Boergens,2021年2月25日
%C通过将n个单位立方体排列成长方体可以获得的不同面数_克里斯·米尔森(Chris W.Milson),2021年3月14日
%H Antti Karttunen,n的表,n的a(n)=1..10000</a>
%H Chris W.Milson,<a href=“https://github.com/chrismilson/cuboid-areas(网址:https://github.com/chrismilson/cuboid-areas)“>构建长方体</a>
%H Chris W.Milson,<a href=“https://github.com/chrismilson/cuboid-areas/blob/e860c66ec6dfce910bafffe5ead15fa5548169cf/main.py#L49“>A(n)的更快算法</a>
%F a(n)=Sum_{m|n}A038548(m)=Sum _{m|n}天花板(d(m)/2),其中d(m)=m的除数(A000005)_Manfred Boergens,2021年2月25日
%F a(n)=总和{d|n}A135539(d,n/d)_Ridouane Oudra,2021年7月10日
%F a(n)=(A007425(n)+A046951(n))/2.-_Ridouane Oudra,2024年4月10日
%e 20的除数是1,2,4,5,10,20。有10对除数的乘积除20:1*1=1,1*2=2,1x4=4,1*5=5,1*10=10,1*20=20,2*2=4,2*5=10,2*10=20,4*5=20。同样,有10个乘积可以被20整除:4*5=20,2*10=20,4*10=40,10*10=100,1*20=20,2%20=40,4*20=80,5*20=100,10*20=200,20*20=400。所以a(20)=10。
%t(*first-do*)需要[“Combinatorica`”](*then*)f[n_]:=计数[n/Times@@@Union[Sort/@Tuples[Divisors@n,2]],_Integer];数组[f,91](*_Robert G.Wilson v_,2008年5月31日*)
%t d=除数[n];r=长度[d];总和[上限[长度[除数[d[[j]]]/2],{j,r}](*_Manfred Boergens_,2021年2月25日*)
%o(PARI)
%在序列作者给出两种不同的解释之后,有两种实现:
%o A140773v1(n)={my(ds=除数(n),s=0);对于(i=1,#ds,对于(j=i,#ds),如果(!(n%(ds[i]*ds[j]),s=s+1));s;}
%o A140773v2(n)={my(ds=除数(n),s=0);对于(i=1,#ds,对于(j=i,#ds),如果(!((ds[i]*ds[j])%n),s=s+1));s;}
%2017年5月19日
%o(Python)#参见C.W.Milson链接。
%Y参考A140774。
%Y参考A000005、A034836、A038548、A007425、A046951。
%Y参考A369255(奇偶校验)、A369256(奇数项的位置)。
%K nonn公司
%O 1,2号机组
%A _罗伊查询,2008年5月29日
%E由_Robert G.Wilson v_修订和扩展,2008年5月31日