bzoj1079: [SCOI2008]着色方案(dp)

/*
先确定这是一道dp 
再看范围得结论这是一道多维dp，想起了王八棋。
然后定状态，很多维的话只能一种一维咯，因为要转移所以还要加一维
f[a][b][c][d][e][k] 铅丝维是种类最后是上一块放的是第几种
转移略麻烦，比如这次三块的减了一，两块的就加了一，乘的时候就要减回来。
乘法原理，自行体会，，，，，， 
*/
#include<iostream>
#include<cstdio>

#define ll long long
#define mod 1000000007

using namespace std;
ll f[16][16][16][16][16][6];
int x[6],n;
bool L[16][16][16][16][16][6];

ll dp(int a,int b,int c,int d,int e,int k)
{
    ll t=0;
    if(L[a][b][c][d][e][k])return f[a][b][c][d][e][k];
    if(a+b+c+d+e==0)return 1;
    if(a)  
        t+=(a-(k==2))*dp(a-1,b,c,d,e,1  );
    if(b)  
        t+=(b-(k==3))*dp(a+1,b-1,c,d,e,2);  
    if(c)  
        t+=(c-(k==4))*dp(a,b+1,c-1,d,e,3);  
    if(d)  
        t+=(d-(k==5))*dp(a,b,c+1,d-1,e,4);  
    if(e)  
        t+=e*dp(a,b,c,d+1,e-1,5);  
    L[a][b][c][d][e][k]=1;
    return f[a][b][c][d][e][k]=(t%mod);
}

int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        int t;
        scanf("%d",&t);
        x[t]++;
    }
    printf("%lld",dp(x[1],x[2],x[3],x[4],x[5],0));
    return 0;
}