%I#8 2017年2月23日22:33:55
%S 971121221351441792022072142172273543577507569612,
%电话632639732832212235923622440246625172970309732473342,
%电话:33673374341934253518354535623644367239998938793939854050
%N个数x,对于某些y<sqrt(2x),x^2+x+y^2是奇数本原富足数A136476(N)。
%C A136478中列出了相应的y值。(与此处列出的x值不同,y不会随着A136476(n)而增加。)
%e 97^2+97+7^2=9555=A136476(1)是一个奇本原富足数,因此a(1)=97。
%o(PARI)是(x,n=x^2+x+1,f)={对于步骤(y=1,平方(2*x),2,σ(n+=y*4-4,-1)>2||next;对于(i=1,#f=factor(n)[,1],sigma(n\f[i],-1)>2&&next(2);return(1)))}\\_M.f.Hasler_,2017年2月22日
%Y参考A006038、A136476、A136478、A136479。
%K nonn公司
%O 1,1号机组
%A_Pierre CAMI_,2007年12月31日
%E编辑:M.F.Hasler,2017年2月22日
|