%I#10 2019年9月8日04:34:21
%S 5085013002210610085001450018100221002210022100202465026500,
%电话:3255042250042100481004810048010068500689008410092500103700,
%电话:110500110500110500 140500158600174100201110186100221000224500
%N平方数的σ_2的公共值,σ_2(m)=σ_2。
%H Amiram Eldar,n的表,n的a(n)=1..10000</a>
%F西格玛_2(m)=西格玛_(n),m<n。
%e sigma_2(30)=1^1+2^2+3^2+5^2+6^2+10^2+15^2+30^2=1300,sigma_2(35)=1^2+5^2+7^2+35^2=1300。
%p与(数字理论);sigmap:=proc(p,n)convert(map(proc(z)z^p end,divisors(n)),`+`)end;SA2:=[]:对于z从1到1,do对于m到1500 do m:=sigmap(2,m);对于从m+1到1500的n,n:=σap(2,n);如果N=M,则SA2:=[op(SA2),[M,N,N]]fiod od;SA2;选择(proc(z)z[1]<=1000 end,SA2)#只是想把它缩短一点
%t a[n_]:=模块[{s=DivisorSigma[2,n],ans={}},kmax=天花板[Sqrt[s]];Do[If[DivisorSigma[2,k]==s,AppendTo[ans,s]],{k,n+1,kmax}];ans];s={};Do[v=a[n];Do[AppendTo[s,v[[k]]],{k,1,Length[v]}],{n,1,400}];s(*阿米拉姆·埃尔达尔,2019年9月8日*)
%Y参见A001157、A002025、A002046、A063990。
%Y参见A110926、A110927、A110928。
%K nonn公司
%O 1,1号机组
%A _Walter Kehowski_,2005年9月23日
|