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[N];
long long dp[N][N];
long long c[N][N];
long long d[N];

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;
}