一個有趣的數學題：歐拉計劃第33題，有興趣的同學看過來

火 [復制鏈接]

pacelife

11^#

樓主| 發表于 2015-12-14 22:16:08 | 只看該作者

程序見圖片，沒有進行優化

本帖子中包含更多資源

您需要登錄才可以下載或查看，沒有賬號？注冊會員

回復問題專業，描述清楚伸手黨/灌水/看不懂

使用道具舉報

鋼鐵俠2014

12^#

發表于 2015-12-16 22:16:16 | 只看該作者

鋼鐵俠2014 發表于 2015-12-14 10:46
- H: ~5 ]# M$ P! f/ |很有意思的小題目啊，剛好最近在學python，晚上回家寫個小程序算算

我的程序丑爆了。。。matlab
clear all
a=(10:99)';
b=10:99;
k=1;
c={};
for i=1:90
for j=1:90
      if a(i)<b(j)
         if floor(a(i)/10)==floor(b(j)/10)
            if a(i)/b(j)==(a(i)-10*floor(a(i)/10))/(b(j)-floor(b(j)/10)*10)
                  c{k}=[a(i);b(j)];
                  k=k+1;
            end
         end
      end
end
end
for i=1:90
for j=1:90
      if a(i)<b(j)
         if (a(i)-floor(a(i)/10)*10)==floor(b(j)/10)
            if a(i)/b(j)==floor(a(i)/10)/(b(j)-floor(b(j)/10)*10)
                  c{k}=[a(i);b(j)];
                  k=k+1;
            end
         end
      end
end
end
for i=1:90
for j=1:90
      if a(i)<b(j)
         if floor(a(i)/10)==(b(j)-floor(b(j)/10)*10)
            if a(i)/b(j)==(a(i)-10*floor(a(i)/10))/(floor(b(j)/10))
                  c{k}=[a(i);b(j)];
                  k=k+1;
            end
         end
      end
end
end
for i=1:90
for j=1:90
      if a(i)<b(j)
         if (a(i)-floor(a(i)/10)*10)==(b(j)-floor(b(j)/10)*10)
            if a(i)/b(j)==(floor(a(i)/10))/(floor(b(j)/10))
                  c{k}=[a(i);b(j)];
                  k=k+1;
            end
         end
      end
end
end

[16;64]
[19;95]
[26;65]
[49;98]