HDU1521——排列組合(指數型母函數)

題目鏈接:

       最近一直在做組合數學的東西,也轉載了部分大牛的博客,這道題若有不是很理解的地方,也可看我最近轉載的有關母函數的博客。這道題是一道標準的指數型母函數裸題,直接套用母函數的公式即可,詳見我轉載的這篇博客,代碼比較好懂,就是在基礎母函數代碼的基礎上做了一點改動,在計算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;
}


發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章