WWWDOT - GOOGLE = DOTCOM

wayne

10# chyanog 这个容易解释。 这个Mathematica应该是用通用算法解出来的，具体应该就是线性规划。 如果你还加了互不相等的条件的话，那么这个线性规划的约束条件 一下子就多了C(9,2) =36个不等式，加上原先的9个(或者10个)， 就是46个。 一个是10个方程的线性规划， 一个是46个方程的线性规划， 应该高下立判了吧

wayne

11# wayne 加了互不相等的条件，Mathematica的AAS(算法自动选择) 不一定仍然使用线性规划算法。 在找不到匹配算法的时候， Mathematica只能brute force了

wayne

基于hujunhua的代码， 我稍稍做了改动，经过漫长的等待，Mathematica终于解出答案了， 发现Mathematica好容易被骗哦， 。

2. f = FromDigits;(*函数名缩写*)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] + Sin[wayne]^2 && wayne == 0 && And @@ (0 <= # <= 9 & /@ v), Flatten[{v, wayne}], Integers];
7. sele = Select[solv, Length[Union[v /. #]] == 9 &];
8. f /@ {x, y, z} /. sele

复制代码

chyanog

11# wayne Oh,知道了

chyanog

Solve有些时候很慢，记得仿真论坛有过一道题，用Solve似乎就不适合

用1到9这九个数不重不漏组成等式( )/( )( )+( )/( )( )+( )/( )( )=1

1. Compile[{}, Module[{a, b, c, d, e, f, g, h, i},
2. Do[{a, b, c, d, e, f, g, h, i} = x;
3. If[a < d < g && a/(10 b + c) + d/(10 e + f) + g/(10 h + i) == 1,
4. Print@x], {x, Permutations@Range@9}]]][] // Timing

复制代码