HDU4532(组合DP)
生活随笔
收集整理的這篇文章主要介紹了
HDU4532(组合DP)
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
題目:安排座位
?
解析:http://www.douban.com/note/269136472/
?
#include <iostream> #include <string.h> #include <stdio.h>using namespace std; typedef long long LL;const LL MOD=1000000007;LL a[550]; LL A[550]; LL C[550][550]; LL dp[55][550];void Init() {for(int i=0; i<550; i++){C[i][0]=C[i][i]=1;for(int j=1; j<i; j++)C[i][j]=(C[i-1][j]+C[i-1][j-1])%MOD;}A[0]=1;for(int i=1; i<550; i++)A[i]=(A[i-1]*i)%MOD; }int main() {Init();int T,n,tt=1;cin>>T;while(T--){cin>>n;for(int i=1; i<=n; i++)cin>>a[i];memset(dp,0,sizeof(dp));dp[1][a[1]-1]=1;LL sum=a[1];for(int i=2; i<=n; i++){for(int j=0; j<sum; j++) //對每一種空位for(int k=1; k<=a[i]; k++) //將a[i]個元素分成k組for(int u=0; u<=j && u<=k; u++) //將u組放到前j個空位中dp[i][j-u+a[i]-1-(k-1)]=(dp[i][j-u+a[i]-k]+(((dp[i-1][j]*C[j][u])%MOD*C[sum+1-j][k-u])%MOD*C[a[i]-1][k-1])%MOD)%MOD;sum+=a[i];}printf("Case %d: ",tt++);LL ans=dp[n][0];for(int i=1; i<=n; i++) //對每一組,進行全排列ans=(ans*A[a[i]])%MOD;cout<<ans<<endl;}return 0; }
?
總結
以上是生活随笔為你收集整理的HDU4532(组合DP)的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: HDU4514(非连通图的环判断与图中最
- 下一篇: set与map容器