P4609 [FJOI2016]建筑师（第一类斯特林数）

传送门

没想到连黑题都会有双倍经验的

其实这题本质上是和CF960G Bandit Blues一样的，不过那里是要用分治FFT预处理第一类斯特林数，这里直接打表预处理第一类斯特林数就可以了

//minamoto
#include<bits/stdc++.h>
#define R register
#define fp(i,a,b) for(R int i=a,I=b+1;i<I;++i)
#define fd(i,a,b) for(R int i=a,I=b-1;i>I;--i)
#define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)
using namespace std;
char buf[1<<21],*p1=buf,*p2=buf;
inline char getc(){return p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++;}
int read(){
    R int res,f=1;R char ch;
    while((ch=getc())>'9'||ch<'0')(ch=='-')&&(f=-1);
    for(res=ch-'0';(ch=getc())>='0'&&ch<='9';res=res*10+ch-'0');
    return res*f;
}
char sr[1<<21],z[20];int K=-1,Z=0;
inline void Ot(){fwrite(sr,1,K+1,stdout),K=-1;}
void print(R int x){
    if(K>1<<20)Ot();if(x<0)sr[++K]='-',x=-x;
    while(z[++Z]=x%10+48,x/=10);
    while(sr[++K]=z[Z],--Z);sr[++K]='\n';
}
const int N=50000,M=200,P=1e9+7;
inline int add(R int x,R int y){return x+y>=P?x+y-P:x+y;}
inline int dec(R int x,R int y){return x-y<0?x-y+P:x-y;}
inline int mul(R int x,R int y){return 1ll*x*y-1ll*x*y/P*P;}
int ksm(R int x,R int y){
    R int res=1;
    for(;y;y>>=1,x=mul(x,x))if(y&1)res=mul(res,x);
    return res;
}
int fac[M+5],inv[M+5],s[N+5][M+5],n,a,b;
inline int C(R int n,R int m){return 1ll*fac[n]*inv[m]%P*inv[n-m]%P;}
void init(){
	inv[0]=fac[0]=1;fp(i,1,M)fac[i]=mul(fac[i-1],i);
	inv[M]=ksm(fac[M],P-2);fd(i,M-1,1)inv[i]=mul(inv[i+1],i+1);
	s[0][0]=1;
	fp(i,1,N)fp(j,1,min(i,M))s[i][j]=add(s[i-1][j-1],1ll*s[i-1][j]*(i-1)%P);
}
int main(){
//  freopen("testdata.in","r",stdin);
    int T=read();init();
    while(T--){
    	n=read(),a=read(),b=read();
    	print(mul(s[n-1][a+b-2],C(a+b-2,b-1)));
	}return Ot(),0;
}