for m=1000:9999
m1=fix(m/1000);
m2=rem(fix(m/100),10);
m3=rem(fix(m/10),10);
m4=rem(m,10);
i=1;
if m==m1^4+m2^4+m3^4+m4^4
disp(m)
n(i)=m;
i=i+1;
end
end
结果
clc
clear
a=1000:9999;
% a=0:1000;
i=fix(a/1000).^3+fix(mod(a,1000)/100).^3+fix(mod(a,100)/10).^3+mod(a,10).^3;
aa=a(find(i==a));
sum1=sum(i==a);
一个也没有，程序没有问题，在0:1000测试通过的。
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ：24596024