SGU 441 Set Division(矩阵，斯特灵数)

转载请注明出处，谢谢 http://blog.csdn.net/ACM_cxlove?viewmode=contents           by---cxlove

题目：有N个物品，分到K个不可区分集合。

http://acm.sgu.ru/problem.php?contest=0&problem=441

典型的第二斯特灵数，但是时限很短，N很大，不过K很小，利用递推式

S(N,K)=K*S(N-1,K)+S(N-1.K-1)可以用矩阵解决

行向量为{S(1,0),S{1,1},……S(1,k)}，构造矩阵

1 1 0 0 0……0

0 1 1 0 0……0

0 0 2 1 0……0

0 0 0 3 1……0

………………

总共是k+1行k+1列。行向量只有S(1,1)为1，其余为0，最终的结果便是 第二行，第K列

顺便说句SGU太坑，还有个ID，每次都要重新注册，把ID放在这：053665，省得下次忘了

#include<iostream>
#include<cstring>
#include<queue>
#include<cstdio>
#include<cmath>
#include<algorithm>
#define N 15
#define inf 1<<29
#define MOD 2007
#define LL long long
using namespace std;
struct Matrix{
    int m[N][N];
}init;
Matrix Mult(Matrix m1,Matrix m2,int n){
    Matrix ans;
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++){
            ans.m[i][j]=0;
            for(int k=0;k<n;k++)
                ans.m[i][j]=(ans.m[i][j]+m1.m[i][k]*m2.m[k][j])%MOD;
        }
    return ans;
}
Matrix Pow(Matrix m1,int b,int n){
    Matrix ans;
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
            ans.m[i][j]=(i==j);
    while(b){
        if(b&1)
            ans=Mult(ans,m1,n);
        m1=Mult(m1,m1,n);
        b>>=1;
    }
    return ans;
}
void Init(int k){
    memset(init.m,0,sizeof(init.m));
    init.m[0][0]=1;
    for(int i=1;i<=k;i++){
        init.m[i-1][i]=1;
        init.m[i][i]=i;
    }
}
int main(){
    int n,k;
    while(cin>>n>>k){
        Init(k);
        init=Pow(init,n-1,k+1);
        cout<<init.m[1][k]<<endl;
    }
    return 0;
}