标签:des blog http os io strong ar div 问题
Time Limit: 1000MS | Memory Limit: 10000K | |
Total Submissions: 29247 | Accepted: 9575 |
Description
Input
Output
Sample Input
2 1 8 4 4 7
Sample Output
0 1 0
转自http://yjq24.blogbus.com/logs/42826226.html
大致上是这样的:有两堆石子,不妨先认为一堆有10,另一堆有15个,双方轮流取走一些石子,合法的取法有如下两种:
1)在一堆石子中取走任意多颗;
2)在两堆石子中取走相同多的任意颗;
约定取走最后一颗石子的人为赢家,求必败态(必胜策略)。
这个可以说是MR.Wythoff(Wythoff于1907年提出此游戏)一生全部的贡献吧,我在一篇日志里就说完有点残酷。这个问题好像被用作编程竞赛的题目,网上有很多把它Label为POJ1067,不过如果学编程的人不知道Beatty定理和Beatty序列 ,他们所做的只能是找规律而已。不熟悉的人可以先在这里 玩几局~
简单分析一下,容易知道两堆石头地位是一样的,我们用余下的石子数(a,b)来表示状态,并画在平面直角坐标系上。
用之前的定理: 有限个结点的无回路有向图有唯一的核 中所述的方法寻找必败态。先标出(0,0),然后划去所有(0,k),(k,0),(k,k)的格点;然后找y=x上方未被划去的格点,标出(1,2),然后划去(1,k),(k,2),(1+k,2+k),同时标出对称点(2,1),划去(2,k),(1,k),(2+k,1+k);然后在未被划去的点中在y=x上方再找出(3,5)。。。按照这样的方法做下去,如果只列出a<=b的必败态的话,前面的一些是(0,0),(1,2),(3,5),(4,7),(6,10),…
接下来就是找规律的过程了,忽略(0,0),记第n组必败态为(a[n],b[n])
命题一:a[n+1]=前n组必败态中未出现过的最小正整数
[分析]:如果a[n+1]不是未出现的数中最小的,那么可以从a[n+1]的状态走到一个使a[n+1]更小的状态,和我们的寻找方法矛盾。
命题二:b[n]=a[n]+n
[分析]:归纳法:若前k个必败态分别为 ,下证:第k+1个必败态为
从该第k+1个必败态出发,一共可能走向三类状态,从左边堆拿走一些,从右边堆拿走一些,或者从两堆中拿走一些.下面证明这三类都是胜态.
情况一:由命题一,任意一个比a[k+1]小的数都在之前的必败态中出现过,一旦把左边堆拿少了,我们只要再拿成那个数相应的必败态即可。
情况二(从右边堆拿走不太多):这使得两堆之间的差变小了,比如拿成了 ,则可再拿成 ;
情况二(从右边堆拿走很多):使得右边一堆比左边一堆更少,这时类似于情况一,比如拿成了 (其中a[m] ;
情况三:比如拿成 ,则可再拿成 .
综上所述,任何从 出发走向的状态都可以走回核中.故原命题成立.
以上两个命题对于确定(a[n],b[n])是完备的了,给定(0,0)然后按照这两个命题,就可以写出(1,2),(3,5),(4,7),…
这样我们得到了这个数列的递推式,以下我们把这两个命题当成是(a[n],b[n])的定义。
先证明两个性质:
性质一:核中的a[n],b[n]遍历所有正整数。
[分析]:由命题一,二可得a[n],b[n]是递增的,且由a[n]的定义显然。
性质二:A={a[n]:n=1,2,3,…},B={b[n]:n=1,2,3,…},则集合A,B不交。
[分析]:由核是内固集,显然。
看到这里大家有没有想到Beatty序列呢,实际上a[n]和b[n]就是一个Beatty序列。
,有 ,解方程
得 ,到此,我们找到了该必败态的通项公式。
实际上这组Beatty序列还有一些别的性质,比如当一个数是Fibonacci数的时候,另一个数也是Fibonacci数;而且两者的比值也越来越接近黄金比,这些性质在得到通项公式之后不难证明。
总的来说,这个问题给我们了哪些启示呢?首先用定理所说的方法找核,然后给出核的规律(递推,或是通项)并且证明。最后附上一张对应的必败态图.
转自http://yjq24.blogbus.com/logs/42653430.html
上次说了胜态和必败态,还记得最后的练习么?桌子上有15个石子,每人每次可以拿去1个或3个石子,拿走最后一个石子的人赢,列出所有的必败态:0,2,4,6,8,10,12,14。说过了状态作为结点可以画一张有向图,下面这张图就是这个游戏所对应的:
定理:有限个结点的无回路有向图有唯一的核。
证明:核可以用如下的方式找出:首先找出没有后继结点的点集P[1](最基本的必败态,比如上图中的结点0),然后找到那些指向P[1]的结点集合为N[1](最基本的胜态,比如上图的结点1和结点3);然后,除去P[1]和N[1]中的点并除去和这些点关联的边,继续寻找没有后继结点的点集P[2](更高级的必败态,比如上图中的结点2),依次类推,则最后的核为P=P[1]并P[2]并…并P[n]。
很容易说明如此找到的核是内固集,也是外固集,满足核的定义,下面说明一下核为什么不是空集:实际上P[1]就不是空集,对一个没有回路的有向图来说,从图上的某一点出发,就无法回到原来到过的点。而图中的点又是有限的,所以最后必将在某个结点终止,故P不是空集。
针对不同的游戏,找核是一个麻烦事。首先生成图,有向边取决于游戏规则,然后当我们要找某个必败态的时候,是要先找到之前所有的必败态的,而这正是一个数学问题和一个编程问题的关键差别。在立方和分解问题[unsolved]中,我的问题的提法都是针对某一个特定输入的n来看是否存在(x,y)满足立方和或者平方和等于n.实际上,如果提法换成,输出对所有不大于n的数中可以被分解的数,那么这种提法更适合计算机去解决,因为本质上来说,两个问题是不一样的。对于前者我只需要知道有关n的情况就可以了,而对后者,却调动了资源去计算所有不大于n的数的情况。虽然他们看起来很相近,但是从道理上来说应该后者的劳动量要大得多,可悲的事情就在于,有时候你要算出n的情况,就不得不算一些比n小的数的情况,而这个计算的数目通常是随着n增大而增大的;另一个可悲的事情是,程序员往往已经习惯了第二种提问方式。数学家希望找到某些必要条件或者充分条件来确定n能否被分解,同样的道理,我们也希望能直接找到必败态的规律,而不真正依赖于象上述定理那样递归的思想从P[1]开始找起,这样来解决问题。
但是,必败态的规律是严格依赖于规则的,这一点对找出必败态的规律来说造成了很大的局限性。这个图的模型在以后还会遇到,到时有更好的方法来寻找必败态。
转自http://yjq24.blogbus.com/logs/42304551.html
高斯取整函数又叫向下取整函数,常见的记法如下: ,既然是向下取整,也就是说[-3.5]=-4,这个取整对负数来说就不是简单地扔掉小数部分,这是要注意的。可以说,高斯取整是联系连续和离散的重要桥梁。
小知识:高斯函数性质
1) x-1<[x]<=x<[x]+1
2) [x+n]=[x]+n,(n为整数)
3) [x]+[y]<=[x+y]<=[x]+[y]+1 //左边由性质2易证,右边利用[x+y]<=x+y<[x]+[y]+2
4) [nx]>=n[x],(n为正整数) //反复利用性质3左边
5) [x/n]=[[x]/n],(n为正整数)
// 换元后等价于证[ny]/n-1<[y]<=[ny]/n,右边由性质4易证,左边有 [ny]/n<=[y]+{y}<[y]+1
欧拉给出过一个很经典的多项式: ,该多项式在n=0,1,2,…,39时产生40个素数。利用高斯取整函数,可以做 一件差不多的事:
,这个函数跳过所有的平方数,而且值域覆盖所有非完全平方数构成的集合,有了上面的这些性质作武器,证明并不难,这里就略去了。
今天的主题还是 Beatty定理 :
正无理数 满足 , 则数列 ; 严格递增, 并且这两个数列构成Z+上的一个分划(也就是它们无交地遍历全体正整数)。
[题解]:
其实作为习题是不难的,显然 ,于是 ,故
Step1.先证明两数列不交:[反证]若 ,有 ,即有 ,
两式相加:得k<m+n<k+1,这和m,n,k都是自然数矛盾;
Step2.再证两数列能取遍所有的正整数:[反证]若k不在 中,则有
于是
两式分别除以 和 后相加:得 ,这和m,n,k都为自然数矛盾.
证毕.
由Beatty定理得到的两个数列称为互质数列,不过别被名称所欺骗,a[n]和b[n]并不能保证对应互质。
//再继续看看必败点
//看完上面的有了点了解吧,代码如下:
#include<iostream> #include<cmath> using namespace std; int main () { int a,b,dif; double p=(sqrt((double)5)+1)/double(2); while(cin>>a>>b) { dif=abs(a-b);//取差值 a=a<b?a:b;//取较小的值 if(a==(int)(p*dif))//判断是不是奇异局势 printf("0\n"); else printf("1\n"); } return 0; }
标签:des blog http os io strong ar div 问题
原文地址:http://www.cnblogs.com/chenyouyou/p/3948009.html