Fibonacci again and again
Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 4248 Accepted Submission(s): 1768
Problem Description
任何一个大学生对菲波那契数列(Fibonacci numbers)应该都不会陌生，它是这样定义的：
F(1)=1;
F(2)=2;
F(n)=F(n-1)+F(n-2)(n>=3);
所以，1,2,3,5,8,13……就是菲波那契数列。
在HDOJ上有不少相关的题目，比如1005 Fibonacci again就是曾经的浙江省赛题。
今天，又一个关于Fibonacci的题目出现了，它是一个小游戏，定义如下：
1、 这是一个二人游戏;
2、 一共有3堆石子，数量分别是m, n, p个；
3、 两人轮流走;
4、 每走一步可以选择任意一堆石子，然后取走f个；
5、 f只能是菲波那契数列中的元素（即每次只能取1，2，3，5，8…等数量）；
6、 最先取光所有石子的人为胜者；
假设双方都使用最优策略，请判断先手的人会赢还是后手的人会赢。
Input
输入数据包含多个测试用例，每个测试用例占一行，包含3个整数m,n,p（1<=m,n,p<=1000）。
m=n=p=0则表示输入结束。
Output
如果先手的人能赢，请输出“Fibo”，否则请输出“Nacci”，每个实例的输出占一行。
Sample Input
1 1 1 1 4 1 0 0 0
Sample Output
Fibo Nacci
Author
lcy
Source
ACM Short Term Exam_2007/12/13
题目：http://acm.hdu.edu.cn/showproblem.php?pid=1848
用SG函数做的第一道题。
对于SG函数，还是有些不太懂，
但是，我看下面说的，就有些明白了：
首先定义mex(minimal excludant)运算，这是施加于一个集合的运算，表示最小的不属于这个集合的非负整数。
例如mex{0,1,2,4}=3、mex{2,3,5}=0、mex{}=0。
对于一个给定的有向无环图，定义关于图的每个顶点的Sprague-Grundy函数g如下：g(x)=mex{ g(y) | y是x的后继 },这里的g（x）即sg[x]
例如：取石子问题，有1堆n个的石子，每次只能取{1，3,4}个石子，先取完石子者胜利，那么各个数的SG值为多少？
sg[0]=0，f[]={1,3,4},
x=1时，可以取走1-f{1}个石子，剩余{0}个，mex{sg[0]}={0}，故sg[1]=1;
x=2时，可以取走2-f{1}个石子，剩余{1}个，mex{sg[1]}={1}，故sg[2]=0；
x=3时，可以取走3-f{1,3}个石子，剩余{2,0}个，mex{sg[2],sg[0]}={0,0}，故sg[3]=1;
x=4时，可以取走4-f{1,3,4}个石子，剩余{3,1,0}个，mex{sg[3],sg[1],sg[0]}={1,1,0},故sg[4]=2;
x=5时，可以取走5-f{1,3,4}个石子，剩余{4,2,1}个，mex{sg[4],sg[2],sg[1]}={2,0,1},故sg[5]=3；
以此类推.....
x 0 1 2 3 4 5 6 7 8....
sg[x] 0 1 0 1 2 3 2 0 1...
计算从1-n范围内的SG值。
f(存储可以走的步数，f[0]表示可以有多少种走法)
f[]需要从小到大排序
1.可选步数为1~m的连续整数，直接取模即可，SG(x) = x % (m+1);
2.可选步数为任意步，SG(x) = x;
3.可选步数为一系列不连续的数，用GetSG()计算
上述是自jumping_frog博文的建立SG模板时的解释，稍后我也会做个SG函数的模板。
这道题，有了上述方法，就简单了。
首先建立f数组，就是Fibonacci数列。
然后预处理求1000以内的SG数组，通过模板：
// 获得SG数组函数模板，t代表f数组的个数，n代表要求的sg数组上限
// f数组就是能取的个数（对于此题就是Fibonacci数列
// 有时，对于t已知就不需要单独传参
void get_sg(int t,int n)
{
int i,j;
memset(sg,0,sizeof(sg));
for(i=1;i<=n;i++)
{
memset(mex,0,sizeof(mex));
// 对于属于g(x)后继的数置1
for( j=1 ;j<=t && fib[j]<=i ;j++ )
mex[sg[i-fib[j]]]=1;
// 找到最小不属于该集合的数
for( j=0 ; j<=n ; j++ )
if(!mex[j])
break;
sg[i] = j;
}
}
SG的题，很多都可以看成是多个Nim博弈。
然后就可以分析奇异态，非奇异态来确定答案了。
关于，Nim博弈可见博客：http://blog.csdn.net/lttree/article/details/24874819
然后就是此题完整代码：
/************************************************
*************************************************
* Author:Tree *
*From :http://blog.csdn.net/lttree *
* Title : Fibonacci again and again *
*Source: hdu 1848 *
* Hint : SG *
*************************************************
*************************************************/
#include <stdio.h>
#include <string.h>
int fib[21]; //fib保存Fibonacci数列
int sg[1001];//sg[]来保存SG值
bool mex[1001];//mex{}
// 构建SG数组，函数各步骤意义详见上面模板
void get_sg(int n)
{
int i,j;
memset(sg,0,sizeof(sg));
for(i=1;i<=n;i++)
{
memset(mex,0,sizeof(mex));
for( j=1 ;fib[j]<=i ;j++ )
mex[sg[i-fib[j]]]=1;
for( j=0 ; j<=n ; j++ )
if(!mex[j])
break;
sg[i] = j;
}
}
int main()
{
int i,m,n,p;
// 构建Fibonacci数列
fib[0]=1,fib[1]=1;
for(i=2;i<21;++i) fib[i]=fib[i-1]+fib[i-2];
// 预处理获得sg数组
get_sg(1000);
while( scanf("%d%d%d",&m,&n,&p) && m+n+p )
{
if( (sg[m]^sg[n]^sg[p])==0 ) printf("Nacci\n");
else printf("Fibo\n");
}
return 0;
}
ACM-SG函数之Fibonacci again and again——hdu1848,布布扣,bubuko.com