【NOIP2013模拟】粉刷匠 题解代码
原題
Description
赫克托是一個魁梧的粉刷匠,而且非常喜歡思考= =
現在,神廟里有N根排列成一直線的石柱,從1到N標號,長老要求用油漆將這些石柱重新粉刷一遍。赫克托有K桶顏色各不相同的油漆,第i桶油漆恰好可以粉刷Ci根石柱,并且,C1+C2+C3…CK=N(即粉刷N根石柱正好用完所有的油漆)。長老為了刁難赫克托,要求相鄰的石柱顏色不能相同。
喜歡思考的赫克托不僅沒有立刻開始粉刷,反而開始琢磨一些奇怪的問題,比如,一共有多少種粉刷的方案?
為了讓赫克托盡快開始粉刷,請你盡快告訴他答案。
Input
第一行一個正整數T,表示測試數據組數
對于每一組測試數據數據:
第1行:一個正整數K
第2行:K個正整數,表示第i桶油漆可以粉刷的石柱個數,Ci。
Output
對于每組輸入數據,輸出一行一個整數,表示粉刷的方案數mod 1000000007。
Sample Input
3
3
1 2 3
5
2 2 2 2 2
10
1 1 2 2 3 3 4 4 5 5
Sample Output
10
39480
85937576
Data Constraint
30% N≤10,T≤5
50% N≤15,T≤5
80% K≤15,Ci≤5,T≤500
100% K≤15,Ci≤6,T≤2000
題解
每一步處理一種顏色。
設 F[i][j] 表示涂到第i種顏色,有j對相鄰的柱子顏色相同。
記涂到第i種顏色之后一共有 S 塊。
對于 i+1 這種顏色,能粉刷 A 次。
設把這種顏色分成 k 塊,插進序列里面。
插進的位置中,有 l 個位置剛好插進在j對顏色相同的柱子中間( l≤j 且 l≤k )。
新狀態(tài)的 j:原有的 j,新增的 A?k,減少的 l。
狀態(tài)轉移:
把顏色分成 k 塊的方案數:C(A?1,k?1)
把塊插進 l 個位置的方案數:C(j,l)
剩余塊的處理方式:C(S+1?j,k?l)
轉移方程如下:
F[i+1][j+A?k?l]+=F[i][j]?C(j,l)?C(A?1,k?1)?C(S+1?j,k?l)
時間復雜度O(K?N?C?C?T)。(理論復雜度與上一解法一致)
優(yōu)化
我們看到 T=2000 這一條件極不和諧。
慶幸的是,將每個詢問內的數字打亂順序不影響答案。
我們可以將每個詢問內部進行排序,再按字典序排序之后再處理,順序處理。
處理詢問 i+1 時,保留處理i時的f值。
若詢問i與 i+1 的詢問前x個數都相等,則對于x之前的f值都是相等的。
由于按字典序排序了,順序處理會減少大多數的處理。
Code
#include<cstdio> #include<cstring> using namespace std; const int mo=1e9+7,N=91; long long f[N][N],c[N][N]; inline int read() {int data=0; char ch=0;while(ch<'0' || ch>'9') ch=getchar();while(ch>='0' && ch<='9') data=data*10+ch-'0',ch=getchar();return data; } int main() {c[0][0]=1;for(int i=1;i<N;i++){c[i][i]=c[i][0]=1;for(int j=1;j<i;j++) c[i][j]=c[i-1][j]+c[i-1][j-1];}int T=read();while(T--){int n=read(),m=0;memset(f,0,sizeof(f));f[0][0]=1;for(int i=1;i<=n;i++){int x=read();for(int j=0;j<=m;j++)if(f[i-1][j])for(int k=1;k<=x;k++){int y=(x<j)?x:j;for(int l=0;l<=y;l++)f[i][j-l+x-k]=(f[i][j-l+x-k]+f[i-1][j]*c[x-1][k-1]*c[j][l]%mo*c[m-j+1][k-l])%mo;}m+=x;}printf("%lld\n",f[n][0]);}return 0; }總結
以上是生活随笔為你收集整理的【NOIP2013模拟】粉刷匠 题解代码的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 【GDOI2014模拟】旅行 题解代码
- 下一篇: 【NOIP2013模拟】Vani和Cl2