找回密码
 欢迎注册
查看: 48275|回复: 14

[提问] WWWDOT - GOOGLE = DOTCOM

[复制链接]
发表于 2013-2-25 19:23:01 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

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

×
C,D,E,G,L,M,O,T,W分别代表代表0-9之间的数字,彼此不能相等,W!=0, E、M可以互换 请问,满足WWWDOT - GOOGLE = DOTCOM的有多少种情况?
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-2-26 08:51:25 | 显示全部楼层
再限定首位的G≠0,D≠0,可得3个解: 555378-177104=378274(E=M) 777589-188103=589486 777589-188106=589483
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-2-26 09:02:37 | 显示全部楼层
允许首位的G,d取0也无增解
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2013-2-26 15:39:45 | 显示全部楼层
2# hujunhua 555378-177104=378274(E=M) 这个应该不算的,不同字母不能相等
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2013-2-27 17:25:56 | 显示全部楼层
2# hujunhua 是用C/C++穷举,还是Mathematica,或者是分析出来的?
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-2-27 23:16:19 | 显示全部楼层
Mathematica 9
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2013-2-28 17:24:13 | 显示全部楼层
  1. Module[{c, d, e, g, l, m, o, t, w, A = 10^Range[5, 0, -1]},
  2. Do[{c, d, e, g, l, m, o, t, w} = x;
  3. If[{w, w, w, d, o, t}.A - {g, o, o, g, l, e}.A == {d, o, t, c, o,
  4. m}.A, Print@x], {x, Permutations[Range[0, 9], {9}]}]]
复制代码
Permutations有了第二个参数时无法Compile
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2013-3-2 20:48:43 | 显示全部楼层
本帖最后由 chyanog 于 2013-3-2 20:50 编辑 想到了一个能Compile的办法,不知道能否更快
  1. cf = Compile[{{arr, _Integer, 2}},
  2. Module[{c, d, e, g, l, m, o, t, w, A, B, C, T = 10^Range[5, 0, -1]},
  3. Do[{c, d, e, g, l, m, o, t, w} = x;
  4. A = {w, w, w, d, o, t};
  5. B = {g, o, o, g, l, e};
  6. C = {d, o, t, c, o, m};
  7. If[A.T - B.T == C.T, Print@{A, B, C}], {x, arr}]],
  8. CompilationTarget -> "C", RuntimeOptions -> "Speed"];
  9. cf[Permutations[Range[0, 9], {9}]] // AbsoluteTiming
复制代码
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-3-2 22:33:57 | 显示全部楼层
还是直接解方程来得快。
  1. f = FromDigits;(*函数名缩写*)
  2. v = {c, d, e, g, l, m, o, t, w};
  3. x = {w, w, w, d, o, t};
  4. y = {g, o, o, g, l, e};
  5. z = {d, o, t, c, o, m};
  6. solv = Solve[f[x] == f[y] + f[z] && And @@ (0 <= # <= 9 & /@ v), v, Integers];
  7. sele = Select[solv,Length[Union[v /. #]] == 9 &];
  8. f /@ {x, y, z} /. sele
复制代码
  1. {{777589, 188103, 589486}, {777589, 188106, 589483}}
复制代码
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2013-3-3 00:23:53 | 显示全部楼层
9# hujunhua 很好。不过我尝试过这样解方程,这个速度就极慢了
  1. Solve[f[x] - f[y] == f[z] && And @@ (0 <= # <= 9 & /@ v) &&
  2. LogicalExpand[Unequal @@ v], v, Integers] // Timing
复制代码
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-11-22 02:01 , Processed in 0.031889 second(s), 16 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

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