找回密码
 注册会员

QQ登录

只需一步,快速开始

搜索
楼主: pacelife

一个有趣的数学题:欧拉计划第33题,有兴趣的同学看过来

  [复制链接]
 楼主| 发表于 2015-12-14 22:16:08 | 显示全部楼层
程序见图片,没有进行优化
3 x2 j: y" @5 j( z* Q+ ^7 [3 l6 q- g9 U

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册会员

×
发表于 2015-12-16 22:16:16 | 显示全部楼层
钢铁侠2014 发表于 2015-12-14 10:46
1 H7 j( o* p$ t( J1 t/ d很有意思的小题目啊,刚好最近在学python,晚上回家写个小程序算算
1 X" k1 }' U" B% \1 r
我的程序丑爆了。。。matlab
# ]6 `) }  E5 B: ^& z3 P! \6 ^( Iclear all
; D4 g  ?; ]+ v. s* U" f# \a=(10:99)';) l3 Z6 T. v6 y
b=10:99;
6 R3 U" d8 I1 O$ O& ?/ Vk=1;
' ?) q; U6 _1 Tc={};
0 n/ A" ~3 p8 I: \( _for i=1:90
+ t( o  H4 F7 E    for j=1:90
) B+ r% u# D0 m* r2 `; `4 U        if a(i)<b(j), W, j% Y% s* {9 U7 X* ^" n
            if floor(a(i)/10)==floor(b(j)/10)  q, c' l7 w/ M/ @
                if a(i)/b(j)==(a(i)-10*floor(a(i)/10))/(b(j)-floor(b(j)/10)*10)
, H0 N3 s  B7 g6 K+ a                    c{k}=[a(i);b(j)];5 ~; |: _- w$ K9 r; w! V0 i# Y
                    k=k+1;: h7 _1 j5 v$ {* L
                end) z8 A% ^: W) J  e5 p
            end
# j! Z9 A: y' m2 j8 N! F  s: C        end
  {$ s- X: e4 O+ W' I5 _4 w    end; m5 j3 b5 T# Q
end: `5 w) o/ J. p* p# E2 R' Z, E; |2 j
for i=1:90
/ W! a4 a3 |( Z; S+ Y1 y    for j=1:90
- V2 D0 g0 u( O% @2 I2 ?1 f' ~        if a(i)<b(j)
0 F' n5 a. ~; I6 y! h2 |" ?  t0 k3 R            if (a(i)-floor(a(i)/10)*10)==floor(b(j)/10): _2 B* E4 b, i# _% A- T) A
                if a(i)/b(j)==floor(a(i)/10)/(b(j)-floor(b(j)/10)*10)0 I7 j: a, h9 B. G* a7 b  g
                    c{k}=[a(i);b(j)];4 Y+ a' _" ~5 v& W" d3 F
                    k=k+1;+ \. F1 i* U& G
                end/ G$ U, {* L. _3 p* p
            end
; r+ |1 F3 P" m( [' I        end
9 p% x: i7 K* V# N" j3 K: H1 H    end
2 ~- v' Z- n! R' z6 o2 nend
: }/ u' I1 {$ A8 k7 Lfor i=1:907 Q# g4 O. _8 h. I0 c1 m' q
    for j=1:90& R. T, [1 p6 Z. ?& }% [! g% D
        if a(i)<b(j): N0 ~$ _: o8 _$ o- d
            if floor(a(i)/10)==(b(j)-floor(b(j)/10)*10)$ v2 C! Z) @% X# e6 @- P
                if a(i)/b(j)==(a(i)-10*floor(a(i)/10))/(floor(b(j)/10))
) `# \0 v4 L) ?5 w9 |- O                    c{k}=[a(i);b(j)];7 y+ k# C' T, |  `& I
                    k=k+1;  ~1 E. n8 H, v
                end
  {0 n# \: Y! F' ?$ }" I# {5 c4 n/ z" \% g            end
# p: ?$ \5 n0 ~" ?# N0 @, T% i        end
" a( V4 k7 V- A    end
- y  D8 ?. y; H! Qend
- m4 ^% e% ~0 E$ ?8 {3 @for i=1:900 m; L. S3 J8 T* @4 i3 p6 d& q
    for j=1:901 t7 _# j( \7 {$ P5 U$ O9 _  ]
        if a(i)<b(j)4 w, ]$ t/ v% M
            if (a(i)-floor(a(i)/10)*10)==(b(j)-floor(b(j)/10)*10)( w) ]! @; n6 B; m8 b0 g
                if a(i)/b(j)==(floor(a(i)/10))/(floor(b(j)/10))
) }$ K  p* V/ \% v  K( R                    c{k}=[a(i);b(j)];
8 ^2 I7 r# Q  h9 Z                    k=k+1;$ Q- H1 j4 y2 P! l
                end
. z" Z* J, E! M5 ~            end% T! W* }: S$ f5 X) z- t
        end: h# X! g; F) F1 }7 E9 A9 J
    end
0 x; f) H* v/ E0 a# [& P1 q9 L% ~end
0 {# K7 Z; z& w3 |
5 W2 @# E8 ~4 Y[16;64]. f1 c# c) B1 E, u% ]8 D
[19;95]
+ K5 a0 J& \! N2 C7 X" ?6 |$ ?/ Q0 _& P( D[26;65]
' f* h/ X+ B7 @[49;98]% |* S  S0 z, ]  Y; l% V: k
* \( Z' R+ G7 D

# f; @% H9 Z& M8 a
发表于 2015-12-18 09:05:43 | 显示全部楼层
牛!!!!!!
发表于 2016-11-18 15:09:19 | 显示全部楼层
都是牛人
回复

使用道具 举报

发表于 2018-4-12 16:59:52 | 显示全部楼层
回复

使用道具 举报

发表于 2019-2-18 10:10:48 | 显示全部楼层
学习了都是大神啊& q( {8 U: |- T! r' i; I
发表于 2019-2-20 15:16:14 | 显示全部楼层
前排都是大神,我只会用笔算排除得到答案。
发表于 2019-2-22 11:23:24 | 显示全部楼层
惭愧,连题目都没看懂
发表于 2019-6-6 11:43:40 | 显示全部楼层
感谢分享
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册会员

本版积分规则

Archiver|手机版|小黑屋|机械社区 ( 京ICP备10217105号-1,京ICP证050210号,浙公网安备33038202004372号 )

GMT+8, 2025-9-14 21:03 , Processed in 0.073995 second(s), 13 queries , Gzip On.

Powered by Discuz! X3.5 Licensed

© 2001-2025 Discuz! Team.

快速回复 返回顶部 返回列表