紫罗兰算式（数学代数方面的）

mathematica

我把num01做最大值进行求解，结果发现不存在解，当然这是计算机程序或精度的问题吧。但是我确实提供了一种求解的办法！

wayne

9# mathematica 1、本题是整数规划 2、lingo好像不能求出所有解吧？

mathematica

1. model:
2. sets:
3. lie/1..10/:shuzi;!行标题是10个数字,按行排放在表的上侧;
4. hang/C,L,O,V,E,R,U,S,I,T/:zimu;!列标题10个字母,按列排放在表的左侧;
5. links(hang,lie):bianliang;!形成100个变量;
6. endsets
7. data:
8. shuzi=1,2,3,4,5,6,7,8,9,0;
9. enddata
10. !表的每行变量的求和都是1,用来使得每个字母只被分到一个数字去;
11. @for(hang(i):@sum(lie(j):bianliang(i,j))=1);
12. !表的每列变量的求和都是1,用来保证使得每个数字只接收一个字母;
13. @for(lie(j):@sum(hang(i):bianliang(i,j))=1);
14. !每行字母所得到的数字;
15. @for(hang(i):@sum(lie(j):bianliang(i,j)*shuzi(j))=zimu(i));
16. !每个变量都是0/1变量;
17. @for(links(i,j):@bin(bianliang(i,j)));
18. !形成第一个数字;
19. num01=10^5*zimu(1)+10^4*zimu(2)+10^3*zimu(3)+100*zimu(4)+10*zimu(5)+zimu(6);
20. !形成第二个数字;
21. num02=10^5*zimu(1)+10^4*zimu(6)+10^3*zimu(3)+100*zimu(1)+10*zimu(7)+zimu(8);
22. !形成第三个数字;
23. num03=10^5*zimu(4)+10^4*zimu(9)+10^3*zimu(3)+100*zimu(2)+10*zimu(5)+zimu(10);
24. num01+num02-num03=0;
25. !min=num01;!不需要目标函数其实也能求解的;
26. end

复制代码

其实不需要目标函数反而计算得更快了，真是无法理解，在我的电脑上22秒这样就计算出来了

yongle

你好，昨天晚上睡前在读者上看到，刚才算出来了，其实方法不难， 先把几个数字确定，E=1,U=9，O=0。然后再把几个逻辑搞清楚，首先，R+S肯定是进位的，E+U肯定是进位的，V+C肯定是不进位的，O+O是不进位的。不能确定的是L+R进不进位，这要做两种情况分别演算。假设L+R不进位，则V=2C,L=3C+1,则C=2,L=7,因为L+R不进位，所以R只能是0、1或者2，而这三个数字都已经使用，所以L+R不进位是不成立的，也即L+R必须进位。那么则有：V=2C+1,L=3C+2，C只能是0、1、2之中数字，0、1已经使用，所以C=2,V=5,L=8，这样没使用的数字只有3、4、6、7对应I、R、S、T这其中的数字，很容易推出最终结果。

wayne

14# yongle E为啥是1呢

yongle

呵呵,我试了一下，只有1加别的数字才有可能等于自己，这个没什么道理可言吧。

wayne

16# yongle 没明白： E+U=E ，只能说明 U=9，R+S=T+10 或者 U=0，R+S=T 而E可以是任意的。

yongle

有道理，比较巧合吧，我开始就认定E=1，看来还得再研究。

没——问题

2. /**************************************************************************
3. * 题目:紫罗兰算式
4. * 设第一个数是“CLOVER”；
5. * 第二个数是“CROCUS”；
6. * 第三个数是“VIOLET”；
7. * 其中第一个数与第二个数的和是第三个数，
8. * 注意：相同的字母表示同一个数字，不同的字母
9. * 表示不同的数字。每个字母的取值是从0到9的整数。
10. * http://bbs.emath.ac.cn/viewthread.php?tid=2606&extra=&page=1
11. * 分析:
12. * 设所求的三个数为X,Y,Z,通过写做10的幂的和的形式可得到方程:
13. * -99900*V+10*U-T+S+10001*R+1000*O+9900*L-10000*I+200100*C==0
14. * 根据系数的正负分为两组:
15. * 10*U+S+10001*R+1000*O+9900*L+200100*C==99900*V+T+10000*I
16. * 对等式两边分别将所有可能的值写入数组left和right,并分别排序
17. * 数组的每个元素除存储等式一侧的值外还需存储生成此值的变元的值
18. * 寻找两组中相等的数,时间复杂度:10^6+10^3
19. *************************************************************************/
20. #include <stdlib.h>
21. #include <stdio.h>
23. typedef struct {int value,id;} group;//value为等式一侧的值,id为一侧变量的值依次排列在一起组成的十进制数字的值
24. int cmp(const void *a, const void *b){return ((group *)a)->value - ((group *)b)->value;}//从小到大排序
26. int main()
27. {
28. int u,s,r,o,l,c,v,t,i;
29. group left[1000000],*pl=left,right[1000],*pr=right;
30. for(v=0;v<10;v++)for(t=0;t<10;t++)for(i=0;i<10;i++)
31. {
32. pr->value=99900*v+t+10000*i;
33. pr->id=v*100+10*t+i;
34. pr++;
35. for(o=0;o<10;o++)for(l=0;l<10;l++)for(c=0;c<10;c++)
36. {
37. pl->value=10*v+t+10001*i+1000*o+9900*l+200100*c;
38. pl->id=v*100000+t*10000+i*1000+o*100+l*10+c;
39. pl++;
40. }
41. }
43. qsort(left, 1000000, sizeof(group), cmp);
44. qsort(right, 1000, sizeof(group), cmp);
46. printf("ulrolc,vti\n----------\n");
47. for(pl--,pr--;pl!=left&&pr!=right;)
48. {
49. if(pl->value==pr->value)
50. //输出变元的值,此处没有排除使得X,Y,Z中出现前导零的答案.
51. printf("%06d,%03d\n",pl->id,pr->id),pl--,pr--;
52. if(pl->value > pr->value)pl--;else pr--;
53. }
55. return 0;
56. }

复制代码

我这里是1秒以内 real 0m0.446s user 0m0.340s sys 0m0.012s

没——问题

上面的代码没有输出e的值,因为e是任意的 这从注释中的方程中就可以看出