用matlab的脚本文件编辑程序找出1000内的勾股数 即找a,b,c，使a^2+b^2=c^2

用matlab的脚本文件编辑程序找出1000内的勾股数 即找a,b,c，使a^2+b^2=c^2

a=[];
b=[];
c=[];
for i=1:1000
for j=1:1000
for k=1:1000
if i^2+j^2==k^2
a=[a i];
b=[b j];
c=[c k];
end
end
end
end

用for循环太慢了，看我的 >> a = 1:1000; >> repa = repmat(a,1000,1);%复制向量成矩阵1000x1000 >> repb = repa'; >> c = sqrt(repa.^2 + repb.^2); >> csign = (c == fix(c));%判别是否a^2+b^2=c^2 >> sum(csign(:))/2 ans = 1034 >> csign = tril(csign);%为避免重复，只取上半部分 >> tmp = [repa(csign),repb(csign),c(csign)];%三列，第一列为a，第二列为b，第三列为c