生活随笔
收集整理的這篇文章主要介紹了
hdu 1521 排列组合 多重集排列 + 指数生成函数
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
傳送門
文章目錄
題意：
思路：
顯然是多重集排列數，我們考慮構造指數生成函數，讓后模擬一下多項式乘法即可啦。
由于存在分數，所以直接用doubledoubledouble即可。
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#include<random>
#include<cassert>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid ((tr[u].l+tr[u].r)>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std
;
typedef long long LL
;
typedef unsigned long long ULL
;
typedef pair
<int,int> PII
;const int N
=1000010,mod
=1e9+7,INF
=0x3f3f3f3f;
const double eps
=1e-6;int n
,m
;
int a
[N
];
double x
[N
],y
[N
];int fun(int x
) {int ans
=1;for(int i
=1;i
<=x
;i
++) ans
*=i
;return ans
;
}int main()
{
while(cin
>>n
>>m
) {for(int i
=1;i
<=n
;i
++) cin
>>a
[i
];for(int j
=0;j
<=m
;j
++) x
[j
]=y
[j
]=0;for(int i
=0;i
<=a
[1];i
++) x
[i
]=1.0/fun(i
);for(int i
=2;i
<=n
;i
++) {for(int j
=0;j
<=m
;j
++) {for(int k
=0;k
+j
<=m
&&k
<=a
[i
];k
++) {(y
[j
+k
]+=x
[j
]/fun(k
));}}for(int j
=0;j
<=m
;j
++) x
[j
]=y
[j
],y
[j
]=0;}printf("%.0f\n",(x
[m
]*fun(m
)));}return 0;
}
總結
以上是生活随笔為你收集整理的hdu 1521 排列组合 多重集排列 + 指数生成函数的全部內容，希望文章能夠幫你解決所遇到的問題。
如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。
