[HDU2294]Pendant

题目：Pendant

链接：http://acm.hdu.edu.cn/showproblem.php?pid=2294

分析：

1）f[i][j]表示长度为i,有j种珍珠的吊坠的数目。
$f[i][j] = (k - j + 1) * f[i - 1][j - 1] + j * f[i - 1][j] $

2）用矩阵来转移.

转移矩阵：

$\left[ \begin{array}{cccccc} 1 & 0 & 0 & ... & 0 & 0 \\ 0 & 1 & k-(2-1) & ... & 0 & 0 \\ 0 & 0 & 2 & ... & 0 & 0 \\ ... & ... & ... & ... & ... & ... \\ 0 & 0 & 0 & ... & k-1 & k-(k-1) \\ 0 & 0 & 0 & ... & 0 & k \end{array} \right] $

状态矩阵：

$\left[ \begin{array}{ccccc} sum[i-1] & f[i][1] & f[i][2] & ... & f[i][k] \end{array} \right] $

初始矩阵：

$\left[ \begin{array}{ccccc} 0 & k(f[i][1]的值) & 0 & ... & 0 \end{array} \right] $

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
using namespace std;
typedef long long LL;
const int MOD=;
struct Matrix{
    LL n,a[][];
    void init(int _n,int f){
        n=_n;
        memset(a,,sizeof a);
        if(f==-)return;
        for(int i=;i<=n;++i)a[i][i]=;
    }
};
Matrix operator*(Matrix& A,Matrix& B){
    Matrix C;C.init(A.n,-);
    for(int i=;i<=C.n;++i)
    for(int j=;j<=C.n;++j)
        for(int k=;k<=C.n;++k){
            C.a[i][j]+=A.a[i][k]*B.a[k][j];
            C.a[i][j]%=MOD;
        }
    return C;
}
Matrix operator^(Matrix A,int n){
    Matrix Rt;Rt.init(A.n,);
    for(;n;n>>=){
        if(n&)Rt=Rt*A;
        A=A*A;
    }
    return Rt;
}
int main(){
    int Case;scanf("%d",&Case);
    Matrix A,T;
    for(int n,k;Case--;){
        scanf("%d%d",&n,&k);
        memset(T.a,,sizeof T.a);
        T.n=k;
        T.a[][]=;T.a[k][]=;
        T.a[][]=;
        for(int i=;i<=k;++i){
            T.a[i][i]=i;
            T.a[i-][i]=k-i+;
        }
        A=T^n;
        printf("%lld\n",A.a[][]*k%MOD);
    }
    return ;
}