找回密码
 欢迎注册
查看: 12414|回复: 7

[讨论] 对对碰自动连续消除次数的期望

[复制链接]
发表于 2009-11-26 14:43:57 | 显示全部楼层 |阅读模式

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

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

×
博客上有人问了一个类似的问题,觉得没有什么好的思路。

题目:
对于n x n大小的游戏地图,游戏基本元素有m个,问正常消除一次导致后续自动连续消除的次数C期望是多大?
解释:
游戏基本元素可以简单理解为游戏中有多少种不同的头像,为了简便分析可以先考虑在坐标i,j处横向消除3个头像导致自动连消的次数。
我的一点思路:假设横向消除3个的话,必然会导致消除行上方整体3列落下,然后消除发生的情况肯定是由于这个下落区域跟不动区域边界移动造成,消除点肯定发生在边界处,假设新补充进来的头像都是等概率均布,依次分析各个接触点的概率。
不过感觉好难。
我做过模拟,对于40x40的地图,如果只有6个元素的话,自动连消基本停不下来,然而对于20x20的地图,一般连消也就是几步,大家能否找到C和m,n的大致关系?
附:对对碰消除规则
二.概述

游戏在 8 × 8 格子的游戏池中进行。每个格子中有一个动物头像。鼠标连续选中两个相邻的头像,它们的位置会互换,互换后如果横排或竖排有 3 个以上相同的头像,则可以消去该头像,并得分。

三• 基本规则

交换
玩家选中相邻(横、竖)的两个头像,则这两个头像的位置发生互换,如果互换成功则消去头像,否则取消位置交换。

消去
玩家选择两个头像进行位置互换,互换后如果横排或竖排有 3 个以上相同的头像,则消去这几个相同的头像,如果互换后没有可以消去的头像,则选中的两个头像换回原来的位置。消去后的空格由上面的头像掉下来补齐。每次消去头像玩家都能得到一定的分数。

连锁
玩家消去头像后,上面的头像掉下来补充空格。如果这时游戏池中有连续摆放(横、竖)的 3 个或 3 个以上相同的头像,则可以消去这些头像,这就是一次连锁。空格被新的头像填充,又可以进行下一次连锁。每次连锁会有加分。

评分

参与人数 1威望 +1 鲜花 +1 收起 理由
KeyTo9_Fans + 1 + 1 很有趣很值得研究的问题

查看全部评分

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2009-11-26 15:11:38 | 显示全部楼层
我经常玩这个游戏,也曾经思考过自动连消与避免无路可走的问题。

但一直没有把这些现象量化,作为一道题目提出来。

感谢楼主给出了清晰的问题描述和合理的参数设置。

我有空仔细研究研究。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2009-11-26 17:49:56 | 显示全部楼层
对于不同的游戏地图,得出的期望值C不一定相同。
如果地图不确定,首先要确定一共用多少种地图?(而计算n*n一共有多少种地图,又是一个不容易求出的问题),每种地图出现的概率又各是多少?(每幅地图出现的概率跟游戏程序的设计有关)
设每幅地图出现的概率为P(x),该幅地图的C值为C(x)
那么总的期望值C=∑ P(x)*C(x)
-----------------------------------------------------------------------------------
为了便于分析,可以假设所有的P(x)都相同。
但是如果让我来设计该游戏,为了简化,我会让每个格子随机设为一个元素,然后利用规则进行自动消除,消除后稳定的状态作为一个地图。所以如果采取这种设计方法,每张地图的发生概率不会都是相同的。
-----------------------------------------------------------------
所以为了研究该问题,楼主还要确定地图产生的方法,以便确定每种地图发生的概率。只有统一地图产生的方法,才可能有相同的结论。
如果楼主只研究某种地图,那么还要确定是哪种地图。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2009-11-27 07:59:24 | 显示全部楼层
3# 056254628
我忘了说明了,
这里的地图说的是n x n的正方形的地图,
初始情况地图里面的元素都是随机的,
然后要分析两种情况,
①就是所有元素都随机的情况下,自动连消期望多大?
②当已经没有自动连消时,交换两个元素之后产生的自动连消期望多大?
个人认为情况①的期望肯定要比期望②大好多。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2010-1-3 23:30:11 | 显示全部楼层
显然地,以下两种情况不需要讨论:

①当n≤2时,C=0;

②当n≥3且m=1时,C=+∞;

对于其他情况,C为有限大的正数。

其中最简单的情况是n=3,m=2,有512种状态。

根据移动规则可以得到各种状态之间的转移规律。

于是得到一个512元1次方程组:

512eqs.PNG

解之,得:

512sols.PNG

所以对于问题①,C= cell.PNG -1=1.92936482260326=17891/9273

问题②还不够清晰:如果有多种可行的交换方案应该交换哪两个?

评分

参与人数 3威望 +6 贡献 +3 鲜花 +3 收起 理由
wayne + 3 + 3 强,又是写bmp文件!!!
mathe + 3 非常强大
northwolves + 3

查看全部评分

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2010-1-3 23:39:46 | 显示全部楼层
显然地,以下两种情况不需要讨论:

①当n≤2时,C=0;

②当n≥3且m=1时,C=+∞;

对于其他情况,C为有限大的正数。

其中最简单的情况是n=3,m=2,有512种状态。

根据移动规则可以得到各种状态之间的转 ...
KeyTo9_Fans 发表于 2010-1-3 23:30

怎么这样强!
请教这是个什么软件做的效果?
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2010-1-3 23:46:13 | 显示全部楼层
是C++代码,才刚把图片熬出锅,还是热的,就放上来了……


  1. #include<cstdio>
  2. #include<memory>

  3. char hd[54]={66,77,38,106,0,0,0,0,0,0,54,0,0,0,40,0,0,0,0,8,0,0,0,6,
  4. 0,0,1,0,24,0,0,0,0,0,-16,105,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};

  5. char nm[10][7][4]={

  6. 0,1,1,0,
  7. 1,0,0,1,
  8. 1,0,0,1,
  9. 1,0,0,1,
  10. 1,0,0,1,
  11. 1,0,0,1,
  12. 0,1,1,0,

  13. 0,0,1,0,
  14. 0,1,1,0,
  15. 0,0,1,0,
  16. 0,0,1,0,
  17. 0,0,1,0,
  18. 0,0,1,0,
  19. 0,0,1,0,

  20. 0,1,1,0,
  21. 1,0,0,1,
  22. 0,0,0,1,
  23. 0,0,1,0,
  24. 0,1,0,0,
  25. 1,0,0,0,
  26. 1,1,1,1,

  27. 0,1,1,0,
  28. 1,0,0,1,
  29. 0,0,0,1,
  30. 0,1,1,0,
  31. 0,0,0,1,
  32. 1,0,0,1,
  33. 0,1,1,0,

  34. 0,0,1,0,
  35. 0,1,1,0,
  36. 1,0,1,0,
  37. 1,0,1,0,
  38. 1,1,1,1,
  39. 0,0,1,0,
  40. 0,0,1,0,

  41. 1,1,1,1,
  42. 1,0,0,0,
  43. 1,1,1,0,
  44. 0,0,0,1,
  45. 0,0,0,1,
  46. 1,0,0,1,
  47. 0,1,1,0,

  48. 0,1,1,0,
  49. 1,0,0,1,
  50. 1,0,0,0,
  51. 1,1,1,0,
  52. 1,0,0,1,
  53. 1,0,0,1,
  54. 0,1,1,0,

  55. 1,1,1,1,
  56. 0,0,0,1,
  57. 0,0,1,0,
  58. 0,0,1,0,
  59. 0,1,0,0,
  60. 0,1,0,0,
  61. 0,1,0,0,

  62. 0,1,1,0,
  63. 1,0,0,1,
  64. 1,0,0,1,
  65. 0,1,1,0,
  66. 1,0,0,1,
  67. 1,0,0,1,
  68. 0,1,1,0,

  69. 0,1,1,0,
  70. 1,0,0,1,
  71. 1,0,0,1,
  72. 0,1,1,1,
  73. 0,0,0,1,
  74. 1,0,0,1,
  75. 0,1,1,0

  76. };

  77. char R[12000][600],G[12000][600],B[12000][600];
  78. int g,h,i,j,k,l,m,ln,c[3][3];

  79. double a[512][513],b[512][513],d;

  80. void pm(int p,int x,int y)
  81. {
  82.         int i,j,q[3][3];
  83.         for(i=2;i>=0;i--)
  84.                 for(j=2;j>=0;j--)
  85.                 {
  86.                         q[i][j]=p%2;
  87.                         p/=2;
  88.                 }
  89.         for(i=0;i<10;i++)
  90.                 for(j=0;j<10;j++)
  91.                         if(i%3&&j%3)
  92.                         {
  93.                                 B[x+i][y+j]=0;
  94.                                 if(q[i/3][j/3])G[x+i][y+j]=0;
  95.                                 else R[x+i][y+j]=0;
  96.                         }
  97.                         else
  98.                         {
  99.                                 R[x+i][y+j]=0;
  100.                                 G[x+i][y+j]=0;
  101.                                 B[x+i][y+j]=0;
  102.                         }
  103. }

  104. void pn(int p,int x,int y)
  105. {
  106.         int i,j;
  107.         for(i=0;i<7;i++)
  108.                 for(j=0;j<4;j++)
  109.                         if(nm[p][i][j])
  110.                         {
  111.                                 R[x+i][y+j]=0;
  112.                                 G[x+i][y+j]=0;
  113.                                 B[x+i][y+j]=0;
  114.                         }
  115. }

  116. void pj(int x,int y)
  117. {
  118.         R[x][y+1]=0;
  119.         G[x][y+1]=0;
  120.         B[x][y+1]=0;
  121.         R[x+1][y]=0;
  122.         G[x+1][y]=0;
  123.         B[x+1][y]=0;
  124.         R[x+1][y+1]=0;
  125.         G[x+1][y+1]=0;
  126.         B[x+1][y+1]=0;
  127.         R[x+1][y+2]=0;
  128.         G[x+1][y+2]=0;
  129.         B[x+1][y+2]=0;
  130.         R[x+2][y+1]=0;
  131.         G[x+2][y+1]=0;
  132.         B[x+2][y+1]=0;
  133. }

  134. void pd(int x,int y)
  135. {
  136.         int i,j;
  137.         for(i=0;i<4;i+=2)
  138.                 for(j=0;j<3;j++)
  139.                 {
  140.                         R[x+i][y+j]=0;
  141.                         G[x+i][y+j]=0;
  142.                         B[x+i][y+j]=0;
  143.                 }
  144. }

  145. int main()
  146. {
  147.         memset(R,-1,sizeof(R));
  148.         memset(G,-1,sizeof(G));
  149.         memset(B,-1,sizeof(B));
  150.         FILE *f=fopen("512eqs.bmp","wb");
  151.         hd[18]=88;
  152.         hd[19]=2;
  153.         hd[22]=-32;
  154.         hd[23]=46;
  155.         for(i=0;i<54;i++)
  156.                 fprintf(f,"%c",hd[i]);
  157.         ln=1;
  158.         for(h=0;h<512;h++)
  159.         {
  160.                 pm(h,ln,16);
  161.                 pd(ln+4,27);
  162.                 g=0;
  163.                 if((h&7)==7||(h&7)==0)g|=7;
  164.                 if((h&56)==56||(h&56)==0)g|=56;
  165.                 if((h&448)==448||(h&448)==0)g|=448;
  166.                 if((h&73)==73||(h&73)==0)g|=73;
  167.                 if((h&146)==146||(h&146)==0)g|=146;
  168.                 if((h&292)==292||(h&292)==0)g|=292;
  169.                 if(!g)
  170.                 {
  171.                         a[h][h]=1;
  172.                         pn(0,ln+2,31);
  173.                         ln+=12;
  174.                 }
  175.                 else
  176.                 {
  177.                         k=256;
  178.                         m=1;
  179.                         for(i=0;i<3;i++)
  180.                                 for(j=0;j<3;j++)
  181.                                 {
  182.                                         c[i][j]=0;
  183.                                         if(h&k)c[i][j]=1;
  184.                                         if(g&k)
  185.                                         {
  186.                                                 c[i][j]=-1;
  187.                                                 m+=m;
  188.                                         }
  189.                                         k/=2;
  190.                                 }
  191.                         a[h][h]=m;
  192.                         a[h][512]=m;
  193.                         if(m<9)pn(8,ln+2,11);
  194.                         else
  195.                                 if(m<99)
  196.                                 {
  197.                                         pn(m/10,ln+2,6);
  198.                                         pn(m%10,ln+2,11);
  199.                                 }
  200.                                 else
  201.                                 {
  202.                                         pn(m/100,ln+2,1);
  203.                                         pn(m/10%10,ln+2,6);
  204.                                         pn(m%10,ln+2,11);
  205.                                 }
  206.                         for(j=0;j<3;j++)
  207.                         {
  208.                                 k=2;
  209.                                 for(i=2;i>=0;i--)
  210.                                         if(c[i][j]>=0)c[k--][j]=c[i][j];
  211.                                 for(;k>=0;k--)c[k][j]=-1;
  212.                         }
  213.                         k=256;
  214.                         g=0;
  215.                         l=0;
  216.                         for(i=0;i<3;i++)
  217.                                 for(j=0;j<3;j++)
  218.                                 {
  219.                                         if(c[i][j]<0)g|=k;
  220.                                         if(c[i][j]>0)l|=k;
  221.                                         k/=2;
  222.                                 }
  223.                         j=-1;
  224.                         for(i=0;i<512;i++)
  225.                                 if((i&g)==i)
  226.                                 {
  227.                                         a[h][i|l]-=1;
  228.                                         if(j>=0)pj(ln+4,j*15+27);
  229.                                         else j=0;
  230.                                         pm(i|l,ln,j*15+31);
  231.                                         j++;
  232.                                         if(j>37)
  233.                                         {
  234.                                                 j=0;
  235.                                                 ln+=12;
  236.                                         }
  237.                                 }
  238.                         pj(ln+4,j*15+27);
  239.                         if(m<9)pn(8,ln+2,j*15+31);
  240.                         else
  241.                                 if(m<99)
  242.                                 {
  243.                                         pn(m/10,ln+2,j*15+31);
  244.                                         pn(m%10,ln+2,j*15+36);
  245.                                 }
  246.                                 else
  247.                                 {
  248.                                         pn(m/100,ln+2,j*15+31);
  249.                                         pn(m/10%10,ln+2,j*15+36);
  250.                                         pn(m%10,ln+2,j*15+41);
  251.                                 }
  252.                         ln+=12;
  253.                 }
  254.         }
  255.         for(i=11999;i>=0;i--)
  256.                 for(j=0;j<600;j++)
  257.                         fprintf(f,"%c%c%c",B[i][j],G[i][j],R[i][j]);
  258.         fclose(f);
  259.         memcpy(b,a,sizeof(a));
  260.         for(i=0;i<512;i++)
  261.                 for(j=i+1;j<512;j++)
  262.                 {
  263.                         for(k=i+1;k<513;k++)
  264.                                 a[j][k]-=a[i][k]*a[j][i]/a[i][i];
  265.                         a[j][i]=0;
  266.                 }
  267.         for(i=511;i>=0;i--)
  268.         {
  269.                 for(j=i+1;j<512;j++)
  270.                 {
  271.                         a[i][512]-=a[j][512]*a[i][j];
  272.                         a[i][j]=0;
  273.                 }
  274.                 a[i][512]/=a[i][i];
  275.                 a[i][i]=1;
  276.         }
  277.         for(i=0;i<512;i++)
  278.         {
  279.                 printf("%d: %.16lf ",i,a[i][512]);
  280.                 d=0;
  281.                 for(j=0;j<512;j++)
  282.                         d+=a[j][512]*b[i][j];
  283.                 printf("%lf %lf\n",d,b[i][512]);
  284.         }
  285.         f=fopen("512sols.bmp","wb");
  286.         memset(R,-1,sizeof(R));
  287.         memset(G,-1,sizeof(G));
  288.         memset(B,-1,sizeof(B));
  289.         for(i=0;i<54;i++)
  290.                 fprintf(f,"%c",hd[i]);
  291.         ln=1;
  292.         for(h=0;h<512;h++)
  293.         {
  294.                 pm(h,ln,1);
  295.                 pd(ln+4,12);
  296.                 d=a[h][512];
  297.                 pn(int(d),ln+2,16);
  298.                 R[ln+8][21]=0;
  299.                 G[ln+8][21]=0;
  300.                 B[ln+8][21]=0;
  301.                 for(j=23;j<103;j+=5)
  302.                 {
  303.                         d=(d-int(d))*10;
  304.                         pn(int(d),ln+2,j);
  305.                 }
  306.                 ln+=12;
  307.         }
  308.         for(i=11999;i>=0;i--)
  309.                 for(j=0;j<600;j++)
  310.                         fprintf(f,"%c%c%c",B[i][j],G[i][j],R[i][j]);
  311.         fclose(f);
  312.         return 0;
  313. }

复制代码

评分

参与人数 1威望 +2 鲜花 +2 收起 理由
winxos + 2 + 2

查看全部评分

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2010-1-4 14:21:59 | 显示全部楼层
效果太好了,假期我一定好好看一下代码
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

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

GMT+8, 2024-5-5 02:10 , Processed in 0.057622 second(s), 20 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

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