#include <cstdio>
#include <cmath>
#include <cstring>
#include <vector>
#include <iostream>
#include <algorithm>
using namespace std;
int c[11][11];//组合数
int f[11],a[11];
void init()
{
int i,j;
memset(c,0,sizeof(c));
c[1][0]=c[1][1]=1;
for(i=2;i<=10;i++)
{
c[i][0]=c[i][i]=1;
for(j=1;j<i;j++)
c[i][j]=c[i-1][j]+c[i-1][j-1];
}
}
int main()
{
init();
int n,m;
while(cin>>n>>m)
{
int i,j,k;
for(i=0;i<n;i++)
cin>>a[i];
memset(f,0,sizeof(f));
f[0]=1;
for(i=0;i<n;i++)
{
for(j=m;j>=0;j--)//枚举已经选择了j件物品
{
for(k=1;k<=min(a[i],m-j);k++)//k不能为0不然就会加上本身
f[k+j]+=c[m-j][k]*f[j]; //每次都只处理用f[j]处理f[j+k],j+k>j,使得不会因为前面处理影响后面的处理
}
}
cout<<f[m]<<endl;
}
return 0;
}
/*
有n种物品,每种取ai种(∑ai==m),则方案数有
ans=m!/(a1!*a2!*...an!)=m!/(a1!*(m-a1)!)*(m-a1)!/(a2!*(m-a1-a2)!)...=c(m,a1)*c(m-a1,a2)*...
含义就是第一种物品有a1个,由于相同,所以就是在m个位置上选a1个,共c(m,a1)种放法,剩下m-a1个位置,
再第二种物品有a2个,再在m-a1个位置选a2个放,有c(m-a1,a2)。。。以此类推就是方案数了。。
然后用变形的背包或则说dp..来解,
f[i]表示已经选好了i件的排列数。。应该这么叫吧,如果没有排列这个就是方案数了。
这里的好处就是(a+b)*c=a*c+b*c;原本应该独立处理每种{ai}序列,最后加上排列数,但这种耗时太大,不能用。
所以就要压缩时间,对于队列两个序列{ai},{bi},如果a1+a2=b1+b2=p,a3=b3=q,;则c(m,a1)*c(m-a1,a2)*c(m-a1-a2,a3)+
c(m,b1)*c(m-b1,b2)*c(m-b1-b2,b3)=(c(m,a1)*c(m-a1,a2)+c(m,b1)*c(m-b1,b2))*c(m-p,q);
这就解释了,为什么上面不用求出每个方案,然后求排列了。。
*/
hdu 1521 排列组合 组合
最新推荐文章于 2021-05-17 06:25:43 发布