原題地址首先，看這么小的范圍就知道，數(shù)學(xué)方法肯定搞不了……又想不到其它模型……只能用狀壓硬搞了囧……
問(wèn)題是，5^15穩(wěn)T，如果能倒過(guò)來(lái)，15^5，就不會(huì)T了。
可以發(fā)現(xiàn)，C值相同的顏色本質(zhì)上是一樣的……因此，只需要保存目前C值為1、2、3、4、5的顏色各有多少種就行了囧……（當(dāng)然在過(guò)程中還會(huì)出現(xiàn)C值為0的，即用完的顏色，不過(guò)0、1、2、3、4、5的和是顏色總數(shù)，而且從下面可以看出，C值為0的確實(shí)“木有用”）；
設(shè)F[s1][s2][s3][s4][s5][v]為涂完前若干個(gè)木塊（這個(gè)個(gè)數(shù)可以通過(guò)s1~s5算出，不過(guò)我們并不需要它囧……）后，C值為1~5的顏色各有s1～s5種，且這若干個(gè)中的
最后一個(gè)涂的顏色還剩的C值為v（顯然0<=v<=4）。
邊界：F[S[1]][S[2]][S[3]][S[4]][S[5]][0]=1，其余為0（S[i]為一開(kāi)始C值為i的顏色種數(shù)）。
計(jì)算F時(shí)，前推后（注意順序，是按照S[5]逆序最先，再依次是S[4]～S[1]，都是逆序，v可任意定序），枚舉下一個(gè)木塊的顏色是現(xiàn)在還剩多少的，如果它與目前的這個(gè)（最后一個(gè)）剩的相同，則要減1，否則不減。具體的方程見(jiàn)代碼。
注意細(xì)節(jié)：枚舉F的s1~s5下標(biāo)時(shí)，都要N（顏色總數(shù)）開(kāi)始枚舉，因?yàn)檫^(guò)程中某些s值會(huì)增加；
本題的啟示就是，在設(shè)計(jì)狀壓DP的時(shí)候，如果正著來(lái)不行，可以反著來(lái)，或許就能設(shè)計(jì)出符合要求的解法。
代碼：
#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
using namespace std;
#define re(i, n) for (int i=0; i<n; i++)
#define re1(i, n) for (int i=1; i<=n; i++)
#define re2(i, l, r) for (int i=l; i<r; i++)
#define re3(i, l, r) for (int i=l; i<=r; i++)
#define rre(i, n) for (int i=n-1; i>=0; i--)
#define rre1(i, n) for (int i=n; i>0; i--)
#define rre2(i, r, l) for (int i=r-1; i>=l; i--)
#define rre3(i, r, l) for (int i=r; i>=l; i--)
#define ll long long
const int n = 5, MAXM = 16, MOD = 1000000007;
int m, S[n + 1];
ll F[MAXM][MAXM][MAXM][MAXM][MAXM][n], res;
void init()
{
scanf("%d", &m); int x;
re(i, m) {scanf("%d", &x); S[x]++;}
}
void solve()
{
F[S[1]][S[2]][S[3]][S[4]][S[5]][0] = 1; int tmp;
rre3(i5, m, 0) rre3(i4, m, 0) rre3(i3, m, 0) rre3(i2, m, 0) rre3(i1, m, 0) re(v, n)
if (F[i1][i2][i3][i4][i5][v]) {
if (i1) {
if (v == 1) tmp = i1 - 1; else tmp = i1;
if (tmp) {
F[i1 - 1][i2][i3][i4][i5][0] += tmp * F[i1][i2][i3][i4][i5][v];
F[i1 - 1][i2][i3][i4][i5][0] %= MOD;
}
}
if (i2) {
if (v == 2) tmp = i2 - 1; else tmp = i2;
if (tmp) {
F[i1 + 1][i2 - 1][i3][i4][i5][1] += tmp * F[i1][i2][i3][i4][i5][v];
F[i1 + 1][i2 - 1][i3][i4][i5][1] %= MOD;
}
}
if (i3) {
if (v == 3) tmp = i3 - 1; else tmp = i3;
if (tmp) {
F[i1][i2 + 1][i3 - 1][i4][i5][2] += tmp * F[i1][i2][i3][i4][i5][v];
F[i1][i2 + 1][i3 - 1][i4][i5][2] %= MOD;
}
}
if (i4) {
if (v == 4) tmp = i4 - 1; else tmp = i4;
if (tmp) {
F[i1][i2][i3 + 1][i4 - 1][i5][3] += tmp * F[i1][i2][i3][i4][i5][v];
F[i1][i2][i3 + 1][i4 - 1][i5][3] %= MOD;
}
}
if (i5) {
F[i1][i2][i3][i4 + 1][i5 - 1][4] += i5 * F[i1][i2][i3][i4][i5][v];
F[i1][i2][i3][i4 + 1][i5 - 1][4] %= MOD;
}
}
res = F[0][0][0][0][0][0];
}
void pri()
{
cout << res << endl;
}
int main()
{
init();
solve();
pri();
return 0;
}