有道难题练习赛 - Sibonacci

昨晚参加了有道难题的练习赛，感觉题目也都不难，参赛人数为3747人，提交次数为20000人次。但是比赛应该说很不成功，服务器非常的很卡，网络很有问题，导致大家很多题目长期刷新不出来，或者出现error页面。我的第二题就是waiting了两个多小时，最后懒得管他了，一直到今天早上发现waiting变成compile error了，真无语。不过第三道题还是做出来了，不过我感觉我的思路有些复杂。。

/* *******************************************************************************************************
有道难题练习赛 - Sibonacci

描述：
菲波那切数列可以用下列的式子表示:
f(1)=1
f(2)=1
f(n)=f(n-1)+f(n-2) (n>=3)

现在我们根据这个规则定义另一种数列 命名为"辛波那切数列", 它是这样定义的：
s(x)=0 (x<0)
s(x)=1 (0<=x<1)
s(x)=s(x-1)+s(x-3.14) (x>=1)

现在需要计算出s(x) MOD 1000000007的值。
输入
第一行有一个正整数T表示有T组测试数据。
接下来T行,每行包含一个数x。
其中 T<=10000, -1000.0<=x<=1000.0
输出
有T行，依次输出每组数据的结果。
样例输入

3
-1
0.667
3.15

样例输出

0
1
2
******************************************************************************************************* */

/* *******************************************************************************************************
解题思路：
首先通过尝试分解发现一个规律，就是：
例如：
S(8.28) 
= S(8.28 - 1) + S(8.28 - 3.14) 
= S(8.28 - 2*1) + 2*S(8.28 - 1 - 3.14) + S(8.28 - 2*3.14)
= S(8.28 - 3*1) + 3*S(8.28 - 2*1 - 3.14) + 3*S(8.28 - 1 - 2*3.14) + S(8.28 - 3*3.14)
= S(8.28 - 4*1) + 4*S(8.28 - 3*1 - 3.14) + 6*S(8.28 - 2*1 - 2*3.14) + 4*S(8.28 - 1 - 3*3.14) + S(8.28 - 4*3.14)
=....

从系数看联想到了杨辉三角形：
　　1
　　1 1
　　1 2 1
　　1 3 3 1
　　1 4 6 4 1
　　1 5 10 10 5 1
　　1 6 15 20 15 6 1
　　1 7 21 35 35 21 7 1 

我们就可以这样定义个杨辉三角形二维数组num[1001][1001]：
    0   1   2   3   4   5   6 .... 1000
0   1   1   1   1   1   1   1 ....
1   1   2   3   4   5   6
2   1   3   6   10  15
3   1   4   10  20 
4   1   5   15
5   1   6
6   1 
...
1000

就有如下规律：
S(x) = ......... = num[a1][b1]*S(x - a1*1 - b1*3.14) + num[a2][b2]*S(x - a2*1 - b2*3.14) + ... + num[an][bn]*S(x - an*1 - bn*3.14)
在根据题目中所给的规律，我们就找出所有合适的ai,bi使得x - ai*1 - bi*3.14在[0, 1)之间，然后将所有的
num[ai][bi]相加并且模除1000000007就可以了。
******************************************************************************************************* */

#include  < iostream >
#include  < cmath >
#include  < cctype >
#include  < string >
#include  < map >
#include  < set >
#include  < vector >
#include  < algorithm >
#include  < list >
#include  < stack >
// #include <stdlib.h>
// #include <iomanip>

using   namespace  std;

// 二维数组保存杨辉三角形
#define  NUM 1001
int  num[NUM][NUM];

int  compute( double  x)
{
     int  res  =   0 , a, b, i, j;
     double  temp;
    a  =  ( int )x;
    b  =  x / 3.14   +   1 ;

     // i表示1的个数，j表示3.14的个数
     for  (j  =  b; j  >=   0 ; j -- )
    {
        temp  =  x -  j * 3.14 ;
         if (temp  >=   0 )
        {
            i  =  ( int )temp;
            res  +=  num[i][j];
             if (res  >=   1000000007 )
                res  %=   1000000007 ;
        }
    }

     return  res;
}

int  main()
{
     int  t, res;
     double  x;
    cin  >>  t;

     // 杨辉三角形
     for ( int  i  =   0 ; i  <  NUM; i ++ )
    {
         for ( int  j  =   0 ; j  <  NUM  -  i; j ++ )
        {
             if (i  ==   0   ||  j  ==   0 )
                num[i][j]  =   1 ;
             else
            {
                num[i][j]  =  num[i - 1 ][j]  +  num[i][j - 1 ];
                 if (num[i][j]  >=   1000000007 )
                    num[i][j]  %=   1000000007 ;
                num[j][i]  =  num[i][j];
            }
        }
    }

     for  ( int  i  =   0 ; i  <  t; i ++ )
    {
        cin  >>  x;
        x  +=   0.000000001 ;   // 防止精度损失
         if (x  <   0 )
            res  =   0 ;
         else   if (x  >=   0   &&  x  <   1.0 )
            res  =   1 ;
         else
            res  =  compute(x);
        cout  <<  res  <<  endl;
    }

     return   0 ;
}

最后结果为：

内存: 4416kB  时间:230ms

还有一种方法感觉很不错：

模拟Fibonacci数列枚举的计算方法，把Sibonacci看成是：s(x) = s(x - 100) + s(x - 314)。
因为s(x) = s(x + 0.01)，所以小数点两位以后的可以忽略。这样就一共有1000*100种情况，全部枚举出来。