hdu2842（矩阵）

题目链接：hdu2842

题意是让求把n个环拆下来需要的最少步数

思路：如果要拆第n个环，那么第n-1个环就必须在竿上，前n-2个环都必须已经被拆下；假设f(n)表示拆第n个环需要的最少步数，那么拆第n个环的时候，第n-1个环在竿上，前n-2个环已经被拆下，那么f(n) = f(n-2)+1，加1是因为拆环的时候需要一步，接下来只剩下第n-1个环了，拆第n-1个环时，第n-2个环必须在竿上，那么就需要在将第n-2个环装在竿上，需要f(n-2)步，再拆第n-1个环，又需要f(n-1)步，则得出递推公式：f(n) = f(n-2)+1+ f(n-2)+f(n-1),即f(n) = 2*f(n-2)+f(n-1)+1；

接下来就是用矩阵二分幂。。。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <cstdlib>
#include <cmath>
using namespace std;
#define mod 200907
int n;
struct node
{
    __int64 map[3][3];
}unit,s;
node Mul(node a,node b)
{
    node c;
    int i,j,k;
    for(i = 0; i < 3; i ++)
    for(j = 0; j < 3; j ++)
    {
        c.map[i][j] = 0;
        for(k = 0; k < 3; k ++)
        c.map[i][j] += (a.map[i][k]*b.map[k][j])%mod;
        c.map[i][j] %= mod;
    }
    return c;
}
void Matrix()
{
    while(n)
    {
        if(n&1) unit = Mul(unit,s);
        n >>= 1;
        s = Mul(s,s);
    }
}
int main()
{
    int i,f[3] = {1,1,2};
    while(scanf("%d",&n),n)
    {
        if(n <= 2)
        {
            printf("%d\n",f[n]);
            continue;
        }
        memset(s.map,0,sizeof(s.map));
        s.map[0][0] = 1; s.map[0][1] = 2; s.map[0][2] = 1;
        s.map[1][0] = 1; s.map[2][2] = 1;
        memset(unit.map,0,sizeof(unit.map));
        for(i = 0; i < 3; i ++)
        unit.map[i][i] = 1;
        n -= 2;
        Matrix();
        __int64 ans = 0;
        for(i = 0; i < 3; i ++)
        ans += unit.map[0][i]*f[2-i], ans %= mod;
        printf("%I64d\n",ans);
    }
    return 0;
}