很简单的dp了,只是没有看清选多个时候的选法是什么,比如选6个 第一堆选 2 第二堆选 1,第三堆选3则组合数为 C(6,2)*C(4,1)*C(3,3)
知道这一点,而且实现很大,直接dp就可以了。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <queue>
#include <stack>
#include <list>
#include <cmath>
#include <sstream>
#define mem(a, val) memset((a), (val), sizeof (a))
#define all(a) (a).begin(), (a).end()
using namespace std;
typedef long long LL;
const int M = 10101;
const int N = 105;
const int MOD = 1000000007;
int c[M][N];
LL d[M][N];
void init(){
c[0][0] = 1;
for(int i=0;i<M;i++){
c[i][0]=1;
if(i<N) c[i][i] = 1;
for(int j=1;j<i&&j<N;j++){
c[i][j]=(c[i-1][j]+c[i-1][j-1])%MOD;
}
}
}
int n,a[N];
int main()
{
init();
int kase=1;
while(scanf("%d",&n)==1){
int sum = 0;
for(int i=1;i<=n;i++) {
scanf("%d",&a[i]);
sum+=a[i];
}
for(int i=n;i>=0;i--)
for(int j=0;j<=sum;j++){
if(i == n){
if(!j) d[j][i]=1;
else d[j][i]=0;
continue;
}
d[j][i] = 0;
for(int k=0;k<=a[i+1];k++){
if(j>=k)
d[j][i]=(d[j][i]+d[j-k][i+1]*c[j][k])%MOD;
else break;
}
}
printf("Case %d: ",kase++);
int res=0;
for(int i=1;i<=sum;i++){
res=(res+d[i][0])%MOD;
}
printf("%d\n",res);
}
return 0;
}