hdu4045 第二类斯特林，至少k个间隔选数

来源：互联网发布：新手淘宝客推广编辑：程序博客网时间：2024/05/21 16:54

分两步，第一步n个数选r个，r个数的间隔数大于等于k-1，这个有公式！搞不懂，只能写进模板

r个不同数分成相同m组，可以为空，箱子装球问题7，第二类斯特林数的部分和

#include<iostream>#include<stdio.h>using namespace std;typedef long long  ll;const int mod=1e9+7;ll s2[1005][1005];//第二类斯特灵void initS(){for(int i=1; i<=1000; i++){s2[i][i]=s2[i][1]=1;//边界}for(int i=2; i<=1000; i++){for(int j=1; j<i; j++){s2[i][j]=s2[i-1][j-1]+1ll*j*s2[i-1][j];//递推式s2[i][j]%=mod;}}}ll B[1005][1005];void initB(){for(int i=1; i<=1000; i++){B[i][1]=1;for(int j=2; j<=i; j++){B[i][j]=B[i][j-1]%mod+s2[i][j]%mod;B[i][j]%=mod;}}}ll c[2005][1005];void initC(){for(int i=0; i<=2000; i++){c[i][0]=c[i][i]=1;}for(int i=1; i<=2000; i++){for(int j=1; j<=i&&j<=1000; j++){c[i][j]=c[i-1][j-1]+c[i-1][j];//递推公式，1000比较小还能用 mod可以不是质数c[i][j]%=mod;}}}int main(){ll n,r,k,m;initS();initB();initC();while(~scanf("%I64d%I64d%I64d%I64d",&n,&r,&k,&m)){ll sum=0;ll rem=n-((r-1)*k+1);if(rem<0){printf("0\n");continue;}//for(int i=1;i<=min(r,m);i++)//{//sum+=s2[r][i]%mod;//sum%=mod;//}if(r<m)m=r;sum=c[rem+r][r]*B[r][m]%mod;cout<<sum%mod<<endl;}}

阅读全文

0 0