Bzoj 1079 着色方案 题解
1079: [SCOI2008]着色方案
Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 2237 Solved: 1361
[Submit][Status][Discuss]
Description
有n个木块排成一行，从左到右依次编号为1~n。你有k种颜色的油漆，其中第i种颜色的油漆足够涂ci个木块。
所有油漆刚好足够涂满所有木块，即c1+c2+...+ck=n。相邻两个木块涂相同色显得很难看，所以你希望统计任意两
个相邻木块颜色不同的着色方案。
Input
第一行为一个正整数k，第二行包含k个整数c1, c2, ... , ck。
Output
输出一个整数，即方案总数模1,000,000,007的结果。
Sample Input
3
1 2 3
1 2 3
Sample Output
10
HINT
100%的数据满足：1 <= k <= 15, 1 <= ci <= 5
这道题挺有趣的，他将DP与记忆化搜索结合在了一起……
我们做一个假设，假设ci都为1，那么这就是一道状压题了，但是ci<=5，虽然仍然不大，但是状压15位显然扑街。
让我们回过头在看最基础的暴力，也就是我们去枚举每一位放的颜色，并将它传递给下一层dfs，如果我们分析一下的话我们会发现每一种剩下可涂数量相同的颜色都可以看作等价的，换句话说涂谁都行。那么，我们将状压的方式换一换，不对，不能叫状压了。改为ci剩余1 2 3 4 5 个的颜色有几种，上一个是谁，也就是f[15][15][15][15][15][7]。空间没问题，至于转移，我们利用dfs的思想，使用记忆化搜索，然后对于每一个f，假设当前状态为f[a][b][c][d][e][la]
则若d!=0且la!=3那么f[a][b][c][d][e][la]+=f[a][b][c][d-1][e+1][2]*d。
若la=3那么f[a][b][c][d][e][la]+=(d-1)*f[a][b][c][d-1][e+1][2]。
其余同理。
不得不说转移数组挺像一道概率DP“抵制克苏恩”的，可惜没有能够应用到这道题来啊。
1 #include <iostream> 2 #include <cstdlib> 3 #include <cstdio> 4 #include <cstring> 5 #include <queue> 6 #include <algorithm> 7 #include <cmath> 8 #include <map> 9 using namespace std; 10 int n,t,a[20],js[10]; 11 long long f[16][16][16][16][16][7],p=1000000007; 12 bool vi[16][16][16][16][16][7]; 13 long long dfs(int a,int b,int c,int d,int e,int la) 14 { 15 16 if(a==0&&b==0&&c==0&&d==0&&e==0)return 1; 17 if(vi[a][b][c][d][e][la])return f[a][b][c][d][e][la]; 18 vi[a][b][c][d][e][la]=1; 19 if(a)f[a][b][c][d][e][la]+=a*dfs(a-1,b+1,c,d,e,5),f[a][b][c][d][e][la]%=p; 20 if(b) 21 { 22 if(b!=1)f[a][b][c][d][e][la]+=(b-(la==5))*dfs(a,b-1,c+1,d,e,4); 23 else if(la!=5)f[a][b][c][d][e][la]+=dfs(a,b-1,c+1,d,e,4); 24 f[a][b][c][d][e][la]%=p; 25 } 26 if(c) 27 { 28 //cout<<a<<' '<<b<<' '<<c<<' '<<d<<' '<<e<<endl; 29 if(c!=1)f[a][b][c][d][e][la]+=(c-(la==4))*dfs(a,b,c-1,d+1,e,3); 30 else if(la!=4) f[a][b][c][d][e][la]+=dfs(a,b,c-1,d+1,e,3); 31 f[a][b][c][d][e][la]%=p; 32 } 33 if(d) 34 { 35 36 if(d!=1)f[a][b][c][d][e][la]+=(d-(la==3))*dfs(a,b,c,d-1,e+1,2); 37 else if(la!=3)f[a][b][c][d][e][la]+=dfs(a,b,c,d-1,e+1,2); 38 f[a][b][c][d][e][la]%=p; 39 } 40 if(e) 41 { 42 if(e!=1)f[a][b][c][d][e][la]+=(e-(la==2))*dfs(a,b,c,d,e-1,1); 43 else if(la!=2)f[a][b][c][d][e][la]+=dfs(a,b,c,d,e-1,1); 44 f[a][b][c][d][e][la]%=p; 45 } 46 return f[a][b][c][d][e][la]; 47 } 48 int main() 49 { 50 scanf("%d",&t); 51 for(int i=1;i<=t;i++) 52 { 53 scanf("%d",&a[i]); 54 n+=a[i]; 55 js[a[i]]++; 56 } 57 long long ans=dfs(js[5],js[4],js[3],js[2],js[1],0); 58 printf("%lld\n",ans); 59 return 0; 60 }
顺便提一句，我之所以还要多开一个vi数组，是为了防止f模p之后为0的情况。
