hdu 4532 湫秋系列故事——安排座位

http://acm.hdu.edu.cn/showproblem.php?pid=4532

细节决定成败呀 ( x * ( y * z ) % mod ) % mod

刚开始想的是 先 ( y * z ) 对 mod 取模 然后再 *x 对mod取模

后来一想不对呀 （y * z ） 要先 *x 后取模呀

脑子短路伤不起呀

dp+数学

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<set>
#include<map>
#include<vector>
#include<queue>
#include<stack>

#define LL long long

using namespace std;
const int N=500;
const long long MOD = 1000000007;
int a
;
long long dp

;
long long c

;
long long d
;

int main()
{
//freopen("data.in","r",stdin);
for(int i=0;i<N;++i)
for(int j=0;j<=i;++j)
{
if(j==0||i==j) c[i][j]=1;
else c[i][j]=(c[i-1][j-1]+c[i-1][j])%MOD;
}
d[0]=1;
for(int i=1;i<N;++i)
d[i]=(d[i-1]*i)%MOD;
int T;
cin>>T;
for(int w=1;w<=T;++w)
{
int n;
cin>>n;
for(int i=0;i<n;++i)
cin>>a[i];
memset(dp,0,sizeof(dp));
dp[0][a[0]-1]=d[a[0]];
int sum=a[0];
for(int i=0;i<n-1;++i)
{
for(int j=0;j<sum;++j)
{
if(dp[i][j]>0)
{
for(int l=0;l<=j&&l<=a[i+1];++l)
{
for(int r=0;r<=(sum+1-j)&&l+r<=a[i+1];++r)
{
if(l+r==0) continue;
int J=(j-l)+(a[i+1]-l-r);
//dp[i+1][J]=(dp[i+1][J]+(dp[i][j]*(c[j][l]*(c[sum+1-j][r]*(d[a[i+1]]*c[a[i+1]-1][l+r-1])%MOD)%MOD)%MOD)%MOD)%MOD;
dp[i+1][J]=(((((c[a[i+1]-1][l+r-1]*d[a[i+1]])%MOD*c[sum+1-j][r])%MOD*c[j][l])%MOD*dp[i][j])%MOD+dp[i+1][J])%MOD;
}
}
}
}
sum+=a[i+1];
}
printf("Case %d: ",w);
cout<<dp[n-1][0]<<endl;
}
return 0;
}