Freda的道路

1305 Freda的道路

时间限制: 1 s
空间限制: 128000 KB
题目等级 : 大师 Master
题解
题目描述 Description
Freda要到Rainbow的城堡去玩了。我们可以认为两座城堡位于同一条数轴上，Freda的城堡坐标是0，Rainbow的城堡坐标是N。正常情况 下，Freda会朝着同一个方向（即Rainbow的城堡相对于Freda的城堡的方向）走若干步之后来到Rainbow的城堡，而且步长都为1或2。可 是，今天Freda在途中遇见了来自上海的小猫Resodo，惊奇之下，居然有一步走反了方向！不过，Freda并没有神智不清，它只有一步走反了方向， 而且这一步的步长也是1或2. 同时，Freda并不会路过Rainbow的城堡而不停下来。当然，Freda是在途中遇到Resodo的，所以它不会在 自己家门口就走错方向。
举个例子，如果Rainbow的城堡坐标是3，那么下面两条路径是合法的：
0->1->2->1->3
0->1->-1->1->3
当然，还有其它的合法路径。下面这些路径则是不合法的：
0->-1->1->3 （Freda不可能第一步就走错方向）
0->1->3（Freda一定是有一步走错方向的）
0->2->1->0->2->3（Freda只有一步是走错方向的）
0->-1->0->3（Freda每步的长度一定是1或2）
0->1->2->4->3（Freda不会越过Rainbow的城堡再回来）
0 -> 1 -> 2 -> 3 -> 2 -> 3（Freda一旦到达了Rainbow的城堡，就会停下来）
你现在需要帮助Freda求出，它一共有多少种方法能够到达Rainbow的城堡呢？

输入描述 Input Description
一行一个整数N，表示Rainbow城堡的坐标

输出描述 Output Description
一行一个整数，表示Freda到Rainbow城堡的不同路径数。由于这个数字可能很大，你只需要输出它mod 1000000007的结果。

样例输入 Sample Input
2

样例输出 Sample Output
5
数据范围及提示 Data Size & Hint
对于第一组样例，如下5条路径是合法的：
0->1->0->2
0->1->-1->0->1->2
0->1->-1->0->2
0->1->0->1->2
0->1->-1->1->2

数据范围与约定
对于10%的数据，N<=20.
对于70%的数据，N<=1000.
对于90%的数据，N<=1000000.
对于100%的数据，N<=10^15.

---

一道有点难度的矩阵乘法的题。
首先我们可以很快的想出动态规划的方法：
f[i][0]表示走到i位置的时候还没有出错时的种类数
f[i][1]表示走到i位置的时候已经出过错了的种类数
这样转移方程就是：
f[i][0]=f[i-2][0]+f[i-1][0]
f[i][1]=f[i-2][1]+f[i-1][1]+f[i+2][0]+f[i+1][0]
这样以后我们就可以过90%的数据了。
对于那个极限数据，我们就只能用矩阵乘法优化了。根据上面的dp方程我们很容易可以看出矩阵

这样我们只要求出第一个矩阵的（n-2）次方以后再乘上

既可以求出答案了。
裸dp代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define F(i,x,y) for(i=x;i<=y;++i)
#define m 1000000007
using namespace std;
int f[1000010][2]={0},n;
int main()
{
    int i,j;
    scanf("%d",&n);
    f[0][1]=1;f[1][1]=1;f[1][0]=1;
    for(i=2;i<=n+3;++i)
    {
        if(i<=n+1) f[i][0]=(f[i-2][0]+f[i-1][0])%m;
        if(i>=4) f[i-2][1]=((f[i-4][1]+f[i-3][1])%m+(f[i-1][0]+f[i][0])%m)%m;
    }
    cout<<f[n+1][1]<<endl;
}

AC代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define F(i,x,y) for(i=x;i<=y;++i)
#define m 1000000007
#define LL long long 
using namespace std;
LL a[4][4]={0},ans[4][4],c[4][4],n;
bool f=true;
long long quickcheng(long long x,long long y)
{
    long long out=0;
    x=x%m;y=y%m;
    if(x>y) swap(x,y);
    while(y)
    {
        if(y&1) out=(out+x)%m;
        y>>=1;
        x=(x+x)%m;
    }
    return out;
}
int main()
{
    LL i,j,out,y,k;
    cin>>n;
    y=n-2;
    if(n==2) cout<<5<<endl;
    else{
        a[0][0]=a[0][1]=a[0][2]=a[0][3]=a[1][0]=a[2][2]=a[2][3]=a[3][2]=1;
        while(y){
            if(y&1){
                if(f){
                    F(i,0,3)
                      F(j,0,3)
                        ans[i][j]=a[i][j];
                    f=false;
                }
                else{
                    F(i,0,3)
                      F(j,0,3){
                        c[i][j]=0;
                        F(k,0,3)
                          c[i][j]=(c[i][j]+quickcheng(a[i][k],ans[k][j]))%m;
                      }
                    F(i,0,3)
                      F(j,0,3)
                        ans[i][j]=c[i][j];
                }
            }
            y>>=1;
            F(i,0,3)
              F(j,0,3){
                c[i][j]=0;
                F(k,0,3)
                  c[i][j]=(c[i][j]+quickcheng(a[i][k],a[k][j]))%m;
              }
            F(i,0,3)
              F(j,0,3)
                a[i][j]=c[i][j];
        }
        out=((quickcheng(ans[0][0],5)+quickcheng(ans[0][1],0))%m+(quickcheng(ans[0][2],8)+quickcheng(ans[0][3],5))%m)%m;
        cout<<out<<endl;
    }
}