%I#25 2023年5月2日11:11:46
%S 3258513768656368187073131231337819683655386579366048,
%电话7209877235378658131073131328137633196608390627390882391137,
%电话:397187397442403747456162456417462722521697781251781506787811846786
%N是3个非零八次幂之和的数字。
%H R.J.Mathar,n表,n=1..8000的a(n)(替换缺少术语的早期b文件)。
%p A003381:=进程(nmax::整数)
%p局部xyzmax,ins,x,x8,y,y8,z,z8;
%p xyzmax:=细胞(根[8](nmax/3));
%pa:={};
%x的p从1到xyzmax do
%px8:=x^8;
%p如果3*x8>nmax,则
%p断裂;
%p end if;
%x do中y的p
%py8:=y^8;
%p如果x8+2*y8>nmax,则
%p断裂;
%p end if;
%y do中z的p
%pz8:=z^8;
%p如果x8+y8+z8>nmax,则
%p断裂;
%p end if;
%p如果x8+y8+z8<=nmax,则
%p a:=并集{x8+y8+z8};
%p结束if;
%p端do:
%p端do:
%p端do:
%p排序(convert(a,list));
%p端程序:
%p最大值:=6755626171875;
%p L:=A003381(nmax);
%p列表文件(L,“b003381.txt”,1);#_R.J.Mathar,2020年8月1日
%t kmax=4*10^12;
%t m=kmax^(1/8)//天花板;
%t表[k=x^8+y^8+z^8;如果[k<=kmax,k,Nothing],{x,1,m},{y,x,m},{z,y,m}]//扁平//联合(*_Jean-François Alcover_,2023年5月2日*)
%Y参见A001016(8次方)。
%K nonn公司
%O 1,1号机组
%A _N.J.A.斯隆_