%I#27 2023年3月17日05:27:40
%S 1,3,4,5,7,8,9,11,13,15,16,17,19,21,23,24,25,27,29,31,32,33,35,36,37,
%电话:39、41、43、45、47、48、49、51、53、55、57、59、60、61、63、64、65、67、69、71、73、75、77,
%U 79,80,81,83,85,87,89,91,93,95,96,97,99100101103107107109111112
%N对k进行编号,使(k的最大除数<=sqrt(k))和(k的最小除数>=sqrt(k)的和为偶数。
%所有奇数正整数和所有完美平方都包含在这个序列中。
%C A139710包含不在此序列中的所有正整数,反之亦然。
%H Harvey P.Dale,n的表,n=1..1000的a(n)</a>
%F{k:A000035(A033676(k)+A033677(k))=0}.-_R.J.Mathar,2008年5月11日
%e 24的除数是1,2,3,4,6,8,12,24。中间的两个除数是4和6。这些的总和是10,这是偶数。因此序列中包含24个。
%p A033676:=进程(n)局部d;对于从地板(sqrt(n))到1乘-1的d,如果n mod d=0,则返回(d);fi;od:结束:A033677:=进程(n)n/A033676(n);结束:isA139711:=进程(n)返回((A033676(n)+A033677(n))模块2=0);结束:对于从1到300的n,如果是A139711(n),则执行printf(“%d,”,n);fi;od:#_R.J.Mathar_,2008年5月11日
%t evdQ[n_]:=模块[{divs=Divisors[n],sr=Sqrt[n]},EvenQ[Max[Select[divs,#<=sr&]]+Min[Select[Civs,#>=sr&]]];选择[Range[120],evdQ](*哈维·P·戴尔,2012年3月5日*)
%t选择[Range[112],IntegerQ[Median[Divisors[#]]&](*_Stefano Spezia_,2023年3月14日*)
%o(PARI)A063655(n)={局部(d);d=除数(n);d[(长度(d)+1)\2]+d[长度(d;
%o表示(n=1120,if(A063655(n)%2==0,print1(n,“,”))\\_G.C.Greubel_,2019年5月31日
%Y参考A000035、A033676、A03367、A063655、A139710。
%K nonn公司
%O 1,2号机组
%2008年4月30日
%E更多来自R.J.Mathar_的条款,2008年5月11日
|