数学研发论坛

 找回密码
 欢迎注册
楼主: mathematica

[提问] 紫罗兰算式(数学代数方面的)

[复制链接]
 楼主| 发表于 2010-8-16 20:18:36 | 显示全部楼层
我把num01做最大值进行求解,结果发现不存在解,当然这是计算机程序或精度的问题吧。但是我确实提供了一种求解的办法!
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2010-8-17 08:54:17 | 显示全部楼层
9# mathematica
1、本题是整数规划
2、lingo好像不能求出所有解吧?
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2010-8-17 17:33:33 | 显示全部楼层
  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秒这样就计算出来了
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2010-8-19 12:11:26 | 显示全部楼层
你好,昨天晚上睡前在读者上看到,刚才算出来了,其实方法不难, 先把几个数字确定,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这其中的数字,很容易推出最终结果。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2010-8-19 12:13:51 | 显示全部楼层
14# yongle

E为啥是1呢
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2010-8-19 12:17:44 | 显示全部楼层
呵呵,我试了一下,只有1加别的数字才有可能等于自己,这个没什么道理可言吧。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2010-8-19 12:38:20 | 显示全部楼层
16# yongle
没明白:
E+U=E  ,只能说明
U=9,R+S=T+10  或者  U=0,R+S=T
而E可以是任意的。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2010-8-19 12:56:37 | 显示全部楼层
有道理,比较巧合吧,我开始就认定E=1,看来还得再研究。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2010-8-19 14:45:59 | 显示全部楼层
本帖最后由 没——问题 于 2010-8-19 15:44 编辑

  1. /**************************************************************************
  2. * 题目:紫罗兰算式
  3. *      设第一个数是“CLOVER”;
  4. *      第二个数是“CROCUS”;
  5. *      第三个数是“VIOLET”;
  6. *      其中第一个数与第二个数的和是第三个数,
  7. *      注意:相同的字母表示同一个数字,不同的字母
  8. *      表示不同的数字。每个字母的取值是从0到9的整数。
  9. *      /viewthread.php?tid=2606&extra=&page=1
  10. * 分析:
  11. *      设所求的三个数为X,Y,Z,通过写做10的幂的和的形式可得到方程:
  12. *      -99900*V+10*U-T+S+10001*R+1000*O+9900*L-10000*I+200100*C==0
  13. *      根据系数的正负分为两组:
  14. *      10*U+S+10001*R+1000*O+9900*L+200100*C==99900*V+T+10000*I
  15. *      对等式两边分别将所有可能的值写入数组left和right,并分别排序
  16. *      数组的每个元素除存储等式一侧的值外还需存储生成此值的变元的值
  17. *      寻找两组中相等的数,时间复杂度:10^6+10^3
  18. *************************************************************************/
  19. #include <stdlib.h>
  20. #include <stdio.h>

  21. typedef struct {int value,id;} group;//value为等式一侧的值,id为一侧变量的值依次排列在一起组成的十进制数字的值
  22. int cmp(const void *a, const void *b){return ((group *)a)->value - ((group *)b)->value;}//从小到大排序

  23. int main()
  24. {
  25.     int u,s,r,o,l,c,v,t,i;
  26.     group left[1000000],*pl=left,right[1000],*pr=right;
  27.     for(v=0;v<10;v++)for(t=0;t<10;t++)for(i=0;i<10;i++)
  28.     {
  29.         pr->value=99900*v+t+10000*i;
  30.         pr->id=v*100+10*t+i;
  31.         pr++;
  32.         for(o=0;o<10;o++)for(l=0;l<10;l++)for(c=0;c<10;c++)
  33.         {
  34.             pl->value=10*v+t+10001*i+1000*o+9900*l+200100*c;
  35.             pl->id=v*100000+t*10000+i*1000+o*100+l*10+c;
  36.             pl++;
  37.         }
  38.     }

  39.     qsort(left, 1000000, sizeof(group), cmp);
  40.     qsort(right, 1000, sizeof(group), cmp);

  41.     printf("ulrolc,vti\n----------\n");
  42.     for(pl--,pr--;pl!=left&&pr!=right;)
  43.     {
  44.         if(pl->value==pr->value)
  45.             //输出变元的值,此处没有排除使得X,Y,Z中出现前导零的答案.
  46.             printf("%06d,%03d\n",pl->id,pr->id),pl--,pr--;
  47.         if(pl->value > pr->value)pl--;else pr--;
  48.     }

  49.     return 0;
  50. }
复制代码
我这里是1秒以内
real        0m0.446s
user        0m0.340s
sys        0m0.012s
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2010-8-19 14:51:42 | 显示全部楼层
上面的代码没有输出e的值,因为e是任意的
这从注释中的方程中就可以看出
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

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

GMT+8, 2019-4-25 05:59 , Processed in 0.052632 second(s), 15 queries .

Powered by Discuz! X3.4

© 2001-2017 Comsenz Inc.

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