转载请注明出处,谢谢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;
}