hdoj 5015 233 Matrix 【矩阵快速幂】

233 Matrix Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Submission(s): 1459    Accepted Submission(s): 855 Problem Description In our daily life we often use 233 to express our feelings. Actually, we may say 2333, 23333, or 233333 ... in the same meaning. And here is the question: Suppose we have a matrix called 233 matrix. In the first line, it would be 233, 2333, 23333... (it means a 0,1 = 233,a 0,2 = 2333,a 0,3 = 23333...) Besides, in 233 matrix, we got a i,j = a i-1,j +a i,j-1( i,j ≠ 0). Now you have known a 1,0,a 2,0,...,a n,0, could you tell me a n,m in the 233 matrix?   Input There are multiple test cases. Please process till EOF. For each case, the first line contains two postive integers n,m(n ≤ 10,m ≤ 10 9). The second line contains n integers, a 1,0,a 2,0,...,a n,0(0 ≤ a i,0 < 2 31).   Output For each case, output a n,m mod 10000007.   Sample Input 1 1 1 2 2 0 0 3 7 23 47 16   Sample Output 234 2799 72937 Hint

憋了1个半小时，终于1A了。O(∩_∩)O~~    

题意：给你一个矩阵，该矩阵a[0][0] = 0,，a[0][1] = 233，a[0][2] = 2333，a[0][3] = 2333...。

已知矩阵中元素满足a[i][j] = a[i-1][j] + a[i][j-1】。

现在给你a[1][0]，a[2][0]，...a[N][0]N个数和M，问你a[N][M] % 10000007。

纯自己的思路，大牛勿喷~~~

分析：这道题肯定是要构造矩阵来求的，关键是找到入手点。

1，首先N<=10，M<=10^9，我们必须要想到构造矩阵的行列是和N有关的，因为M过大是不能构造矩阵来求的，就算是成功构造出矩阵，求解也会TLE。

2，对于a[0][1] = 233，a[0][2] = 2333，可以得到公式a[0][i] = a[0][i-1] * 10 + 3。

3，题目已经给出递推公式a[i][j] = a[i-1][j] + a[i][j-1]，

这样可以知道矩阵求解方案这样的—— ... a[i-1][j] ... a[i][j-1] ...  * |某个矩阵| =  ... a[i][j] ...。

我们会发现，递推公式都是由前列 到 后列。这样就启发我们利用矩阵，每次由前列推导后列。

其实我们也可以这样想，用矩阵来实现这道题的公式推导肯定只有两种——行推行、列推列。对于这道题目，行推行的话必然要求出a[i][0]~a[i][M]，而M最最复杂可以为10^9，这样看来，行推行肯定不行的。而对于列推列，我们只需要求出a[0][i] ~ a[N][i]就可以了，而且N最大也只不过10，所以我们只有利用列推列的方法来实现。

因为公式a[i][j]的推导还需要知道a[i][j-1]，所以在由第j-2列求出第j-1列时，我们要保留第j-1列的数据。

有了公式和递推方向，构造矩阵的思路就已经很明了了，剩下就是用程序实现了。

说到这不得不吐槽一下，这个题目用程序构造出矩阵还真是有点麻烦。。。

这里以N = 3举个构造矩阵的例子。看图           

                                                                               初始矩阵                                F矩阵

图可能有点小，没办法啊，这样更直观些。 想仔细看的话，把图拉到一个新页面会更好一些。

下面就是求出初始矩阵的M-1次幂，然后直接与构造的F矩阵相乘，新矩阵的最后一位就是答案。

因为怕错，直接用long long了。

AC代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
#define LL long long
#define MOD 10000007
using namespace std;
struct Matrix
{
    LL a[30][30];
    int r, c;
};
Matrix ori, res;
void init(int n)
{
    ori.r = ori.c = res.r = res.c = n;
    memset(ori.a, 0, sizeof(ori.a));
    memset(res.a, 0, sizeof(res.a));
    for(int i = 0; i < n; i++)
    {
        res.a[i][i] = 1;//构造单位矩阵
        if(i <= n/2)//构造初始矩阵
            ori.a[i][i] = 1;
        else
        {
            ori.a[0][i] = 3, ori.a[n/2+1][i] = 10;
            int k = i-n/2-1;
            for(int j = n/2+2; k ; j++, k--)
                ori.a[j][i] = 1;
        }
    }
}
Matrix muitl(Matrix x, Matrix y)
{
    Matrix z;
    memset(z.a, 0, sizeof(z.a));
    z.r = x.r, z.c = y.c;
    for(int i = 0; i < x.r; i++)
    {
        for(int k = 0; k < x.c; k++)
        {
            if(x.a[i][k] == 0) continue;
            for(int j = 0; j < y.c; j++)
                z.a[i][j] = (z.a[i][j] + (x.a[i][k] * y.a[k][j]) % MOD) %MOD;
        }
    }
    return z;
}
LL F[30];
void solve(int n)
{
    while(n)
    {
        if(n & 1)
            res = muitl(ori, res);
        ori = muitl(ori, ori);
        n >>= 1;
    }
    LL ans = 0;
    for(int i  = 0; i < res.r; i++)
        ans = (ans + (F[i] * res.a[i][res.c-1]) % MOD) % MOD;
    printf("%lld\n", ans);
}
int main()
{
    int N, M;
    while(scanf("%d%d", &N, &M) != EOF)
    {
        F[0] = 1, F[1] = 0;//构造F矩阵
        for(int i = 2; i <= N+1; i++)
            scanf("%lld", &F[i]), F[i] %= MOD;
        F[N+2] = 233;
        for(int i = N+3; i < 2*(N+1)+1; i++)
            F[i] = (F[i-1] + F[i-N-1]) % MOD;
        init(2*(N+1)+1);//构造初始矩阵 和单位矩阵
        solve(M-1);//初始矩阵的M-1次幂
    }
    return 0;
}