%I#10 2018年9月23日20:57:43
%S 0,0,0,1,0,0,0,1,0,0,0,0,
%温度2,0,0,1,1,3,0,10,3,4,1,0,2,0,2,1,0,1,1,1,0,4,0,12,0,3,1,0,1,
%U 2,3,0,1,0,1,0,3,4,1,0,2,0,1,0,5,2,1,2,3,0,1,6,3,1,4,0,2,0,0,0,0,5,0,5,1
%N余数,当N的最小除数大于或等于N的平方根(A033677(N)),除以N的最大除数小于或等于N(A03367(N)的平方根时。
%C给定n个点,将它们排序为可能的最方形矩形点阵。现在将点排序为正方形点格,其尺寸等于先前矩形的较小尺寸。a(n)是剩余的点数。质数n的a(n)通常为0(质数点上的最平方且唯一的矩形点阵是线性点阵)。a(n)!=0 iff n是A080363的成员。
%H Antti Karttunen,n表,n=1..65537的a(n)</a>
%H Eric Weistein的《数学世界》,<a href=“http://mathworld.wolfram.com/PointLattice.html“>点格</a>。
%H Eric Weistein的《数学世界》,<a href=“http://mathworld.wolfram.com/Divisor.html“>除数。
%e a(6)=1,因为6大于sqrt(6)的最小除数是3,6小于sqrt的最大除数是2,3模2=1
%t num[n_]:=If[OddQ[DivisorSigma[0,n]],Sqrt[n],Divisors[n][[Divisor Sigma[0,n]/2+1]]]den[n]:=如果[OddQ[DivisortSigma\0,n],Sqrt[n]
%o(PARI)
%o A033676(n)=如果(n<2,1,my(d=除数(n));d[(长度(d)+1)\2]);\\来自A033676
%o A033677(n)=(n/A033676(n));
%o A101057(n)=(A033677(n)%A033676(n));\\_Antti Karttune_,2018年9月23日
%Y参见A033676、A033677、A080363。
%K nonn,看
%O 1,15号
%A Joseph Biberstine(jrbibers(AT)indiana.edu),2004年12月17日
%E定义由_Antti Karttunen修订,2018年9月23日
|