最近一直在做組合數學的東西,也轉載了部分大牛的博客,這道題若有不是很理解的地方,也可看我最近轉載的有關母函數的博客。這道題是一道標準的指數型母函數裸題,直接套用母函數的公式即可,詳見我轉載的這篇博客,代碼比較好懂,就是在基礎母函數代碼的基礎上做了一點改動,在計算C2的時候,C1要除以對應的階乘,最後結果乘以對應階乘即可,代碼如下:
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
double a[15]={0};
double c1[105]={0};
double c2[105]={0};
double F(double n)
{
double sum=1.0;
for(double i=1.0;i<=n;i++)
sum*=i;
return sum;
}
int main()
{
//freopen("input.in","r",stdin);
int n, m;
while(scanf("%d%d",&n,&m)!=EOF)
{
int t;
for(int i=1;i<=n;i++)
scanf("%lf",&a[i]);
for(t=1;t<=n;t++)
{
if(a[t]!=0)
{
for(int i=0;i<=a[t];i++)
c1[i]=1.0/F(i);
break;
}
}
for(int i=t+1;i<=n;i++)
{
for(int j=0;j<=100;j++)
for(int k=0;k<=a[i];k++)
c2[k+j]+=(c1[j]/F(k));
for(int j=0;j<=100;j++)
{
c1[j]=c2[j];
c2[j]=0;
}
}
printf("%.0lf\n",c1[m]*F(m));
memset(a,0,sizeof(a));
memset(c1,0,sizeof(c1));
memset(c2,0,sizeof(c2));
}
return 0;
}