首页 > 代码库 > ACM-SG函数之Fibonacci again and again——hdu1848

ACM-SG函数之Fibonacci again and again——hdu1848

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;
}