Poj 3734（矩阵快速幂）

题意：给定N个方块排成一列，现在要用红，蓝，绿，黄，四种颜色的油漆给这些方块染色。求染成红色的方块数和染成绿色的方块数同时为偶数的染色方案的个数。输出对10007取膜；

思路：

我们试着从左往右染色，设染到第i块为止，红绿色都为偶数的方案数位ai，红绿色恰有一个为偶数的方案数为bi，红绿色全为奇数的方案数位ci；

这样染到第i+1块，红绿色都为偶数的有如下两种方案：

1.到第i为止红绿都是偶数，第i+1个方块染黄色或者蓝色；

2.到第i块红绿色恰有一个是奇数，第i+1块将其染成偶数

这样就能得到递推式：ai+1 = 2*ai + bi;

同样的：bi+1 = 2*ai + 2*bi + 2*ci;

              ci+1 = bi + 2*ci;

这样就能得到系数矩阵 2 1 0

                                   2 2 2

                                   0 1 2

用矩阵快速幂求解很容易就能得出答案

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <vector>

using namespace std;
typedef long long ll;
const int mod = 10007;
typedef vector<int>vec;
typedef vector<vec>mat;

mat mul(mat &A,mat &B)
{
    mat C(A.size(),vec(B[0].size()));
    for(int i = 0; i < A.size(); i++)
    {
        for(int k = 0; k < B.size(); k++)
        {
            for(int j = 0; j < B[0].size(); j++)
            {
                C[i][j] = (C[i][j]+A[i][k]*B[k][j])%mod;
            }
        }
    }
    return C;
}

mat pow(mat A,ll n)
{
    mat B(A.size(),vec(A.size()));
    for(int i = 0; i < A.size(); i++)
    {
        B[i][i] = 1;
    }
    while(n)
    {
        if(n&1)
            B = mul(B,A);
        A = mul(A,A);
        n >>= 1;
    }
    return B;
}

ll n;

void solve()
{
   mat A(3,vec(3));
   A[0][0] = 2; A[0][1] = 1; A[0][2] = 0;
   A[1][0] = 2; A[1][1] = 2; A[1][2] = 2;
   A[2][0] = 0; A[2][1] = 1; A[2][2] = 2;
   A = pow(A,n);
   cout<<A[0][0]<<endl;
}

int main()
{
    int T;
    cin >> T;
    while(T--)
    {
        cin >> n;
        solve();
    }
    return 0;
}