題目

【題目描述】

#### 題目背景
一年一度的綜藝節(jié)目《中國好碼農(nóng)》又開始了。本季度，好碼農(nóng)由 Yazid、Zayid、小 R、大 R 四位夢想導(dǎo)師坐鎮(zhèn)，他們都將組建自己的夢想戰(zhàn)隊，并率領(lǐng)隊員向夢想發(fā)起沖擊。

四位導(dǎo)師的**派系**不盡相同，節(jié)目組為了營造看點，又將導(dǎo)師分成了不同的**陣營**，與此同時對不同陣營、不同派系都作出了戰(zhàn)隊總?cè)藬?shù)限制：

- 四位導(dǎo)師分成兩個**陣營**：
- Yazid、小 R 兩位導(dǎo)師組成**藍陣營**，他們兩位的戰(zhàn)隊人數(shù)**總和**不得超過 $C_0$。
- Zayid、大 R 兩位導(dǎo)師組成**紅陣營**，他們兩位的戰(zhàn)隊人數(shù)**總和**不得超過 $C_1$。
- 四位導(dǎo)師分成兩個**派系**：
- Yazid、Zayid 兩位導(dǎo)師屬于**鴨派系**，他們兩位的戰(zhàn)隊人數(shù)**總和**不得超過 $D_0$。
- 小 R、大 R 兩位導(dǎo)師屬于 **R 派系**，他們兩位的戰(zhàn)隊人數(shù)**總和**不得超過 $D_1$。

#### 題目描述
本季好碼農(nóng)邀請到了全國各路學(xué)生精英參賽。他們來自全國 $c$ 個城市的 $n$ 所不同學(xué)校（城市的編號從 $1$ 至 $c$，學(xué)校的編號從 $1$ 至 $n$）。其中，第 $i$ 所學(xué)校所屬的城市編號為 $b_i$，且共有 $s_i$ 名選手參賽。

在「題目背景」中提到的各總?cè)藬?shù)限制之外，本季度《中國好碼農(nóng)》的導(dǎo)師選擇階段有額外規(guī)則如下：

- 來自同**城市**的所有選手必須加入相同的**陣營**。
- 來自同**學(xué)校**的所有選手必須選擇相同的**導(dǎo)師**。

對于導(dǎo)師，大部分學(xué)校的學(xué)生對導(dǎo)師沒有**偏好**。但是有 $k$ 所學(xué)校，其中每所學(xué)校的學(xué)生有且僅有一位他們不喜歡的導(dǎo)師。同一所學(xué)校的學(xué)生不喜歡的導(dǎo)師相同，他們**不會加入他們不喜歡的導(dǎo)師的戰(zhàn)隊**。

面對琳瑯滿目的規(guī)則和選手的偏好，作為好碼農(nóng)忠實觀眾的你想計算出，在所有選手都進行了戰(zhàn)隊選擇后，戰(zhàn)隊組成共有多少種可能的局面？

- 兩種戰(zhàn)隊組成的局面被認為是不同的，當且僅當在存在一所學(xué)校，使得在這兩種局面中這所學(xué)校的選手加入了不同導(dǎo)師的戰(zhàn)隊。
- 由于答案可能很大，你只需輸出可能局面數(shù)對 $998,244,353$ 取模的結(jié)果即可。

【輸入格式】

從標準輸入讀入數(shù)據(jù)。

單個測試點中包含多組數(shù)據(jù)，輸入的第一行包含一個非負整數(shù) $T$ 表示數(shù)據(jù)組數(shù)。接下來依次描述每組數(shù)據(jù)，對于每組數(shù)據(jù)：

- 第 $1$ 行 $2$ 個正整數(shù) $n,c$，分別表示學(xué)校數(shù)目、城市數(shù)目。
- 第 $2$ 行 $4$ 個正整數(shù) $C_0,C_1,D_0,D_1$，分別表示題目中所描述的四個限制。
- 接下來 $n$ 行每行 $2$ 個正整數(shù)：
- 這部分中第 $i$ 行的兩個數(shù)依次為 $b_i,s_i$，分別表示第 $i$ 所學(xué)校的所屬城市以及選手數(shù)目。
- 保證 $b_i \leq c$，$s_i \leq \min\{M, 10\}$。其中 $M=\max{\left\{ C_0,C_1,D_0,D_1\right\}}$。
- 接下來 $1$ 行一個非負整數(shù) $k$，表示選手有偏好的學(xué)校數(shù)目。
- 接下來 $k$ 行，每行 $2$ 個整數(shù) $i,p$，描述編號為 $i$ 的學(xué)校選手有偏好：
- 其中，$p$ 為一個 $0$ 至 $3$ 之間的整數(shù)，描述該校選手不喜歡的導(dǎo)師：0 代表 Yazid，1 代表小 R，2 代表 Zayid，3 代表大 R。
- 保證 $1\leq i\leq n$，且各行的 $i$ 互不相同。

對于輸入的每一行，如果其包含多個數(shù)，則用單個空格將它們隔開。

【輸出格式】

輸出到標準輸出。

依次輸出每組數(shù)據(jù)的答案，對于每組數(shù)據(jù)：

- 一行一個整數(shù)，表示可能局面數(shù)對 $998,244,353$ 取模的結(jié)果。

【樣例輸入】

2
2 1
3 2 2 2
1 1
1 2
1
1 0
4 2
10 30 20 30
1 6
2 4
1 7
2 4
2
2 3
3 1

【樣例輸出】

1
22

【樣例解釋】

對于第 1 組數(shù)據(jù)：

- 唯一的城市 1 包含共 $3$ 名選手，但紅陣營的總?cè)藬?shù)限制為 $2$，無法容納這些選手，因此他們被迫只能選擇藍陣營。
- 在此基礎(chǔ)上，由于 1 號學(xué)校的選手不喜歡 Yazid 老師，因此他們就必須加入 R 派系的小 R 老師麾下。
- 由于 R 派系總?cè)藬?shù)限制為 $2$，因此小 R 老師戰(zhàn)隊無法容納 2 號學(xué)校的選手，所以他們只能被迫加入 Yazid 老師戰(zhàn)隊。
- 綜上所述，可能的局面僅有這一種。

對于第 2 組數(shù)據(jù)：

- 一個顯然的事實是，1 號城市的所有選手都無法加入藍陣營，這是因為 1 號城市的選手總?cè)藬?shù)超過了藍陣營的總?cè)藬?shù)限制，因此他們被迫全部加入紅陣營。
- 對于 2 號城市選手加入藍陣營的情況，稍加計算可得出共有 $15$ 種可能的局面。
- 對于 2 號城市選手加入紅陣營的情況，稍加計算可得出共有 $7$ 種可能的局面。
- 綜上所述，可能的局面數(shù)為 $15+7=22$ 種。

【數(shù)據(jù)范圍與提示】

|測試點|$n$|$c$|$k$|$M$|
|:-:|:-:|:-:|:-:|:-:|
|$1$|$=1$|$=n$|$\le 1$|$=1$|
|$2$|$=10$|$=n$|$\le 10$|$\le 100$|
|$3$|$=20$|$=n$|$=0$|$\le 100$|
|$4$|$=30$|$=n$|$=0$|$\le 100$|
|$5$|$=30$|$=n$|$\le 30$|$\le 500$|
|$6$|$=500$|$=n$|$=0$|$\le 1000$|
|$7$|$=500$|$=30$|$= 30$|$\le 1000$|
|$8$|$=500$|$=n$|$= 30$|$\le 1000$|
|$9$|$=1000$|$=n$|$=0$|$\le 2500$|
|$10$|$=1000$|$=n$|$=30$|$\le 2500$|

其中，$M=\max{\left\{ C_0,C_1,D_0,D_1\right\}}$。

對于所有測試點，保證 $T\leq 5$。

對于所有測試點中的每一組數(shù)據(jù)，保證 $c\leq n\leq 1000$，$k\leq 30$，$M\leq 2500$，$1\leq s_i \leq \min\{M, 10\}$。

**另外，請你注意，數(shù)據(jù)并不保證所有的 $c$ 個城市都有參賽學(xué)校。**

#### 提示
十二省聯(lián)考命題組溫馨提醒您：

**數(shù)據(jù)千萬條，清空第一條。**

**多測不清空，爆零兩行淚。**

題解

這是一個 t1 難度 $ > $ t2 的故事

首先，是一個 $ 50 \% $ 的暴力

記 $ f[i][j][k] $ 表示前 $ i $ 個人，有 $ j $ 個在 $ 1 $ 陣營（藍陣營，下同），有 $ k $ 個在 $ 1 $ 派系（鴨派系，下同）的方案數(shù)，然后隨便轉(zhuǎn)移，效率 $ O(nm^2) $

然后發(fā)現(xiàn)，如果沒有限制的話，陣營和派系的人數(shù)是互不影響的，接著就可以根據(jù)這個搞事情

記 $ g1[i][j] $ 表示前 $ i $ 個城市（沒有限制的），有 $ j $ 個在 $ 1 $ 陣營的方案數(shù)

$ g1[i][j]=g[i-1][j]+g[i-1][j-A[i]],A[i] $ 表示第 $ i $ 個城市的人數(shù)

記 $ g2[i][j] $ 表示前 $ i $ 個學(xué)校（沒有限制的），有 $ j $ 個在 $ 1 $ 派系的方案數(shù)

$ g1[i][j]=g[i-1][j]+g[i-1][j-B[i]],B[i] $ 表示第 $ i $ 個學(xué)校的人數(shù)

發(fā)現(xiàn) $ g1, g2 $ 是一個背包，可以壓成一位來搞

記 $ f[j][k] $ 同之前的定義，滾掉第一維，用于轉(zhuǎn)移有限制的城市和學(xué)校

對于城市，不同陣營的老師對派系的影響是不同的，因此記 $ h=f $ 記錄不同派系的轉(zhuǎn)移

先是對學(xué)校的轉(zhuǎn)移( $ ban[i] $ 為 $ i $ 學(xué)校不喜歡的導(dǎo)師)

$ f[j][k]=\left\{ \begin{matrix} f[j][k-B[v],ban[i]!=0 \\f[j][k],ban[i]!=1 \end{matrix} \right\} $

$ g[j][k]=\left\{ \begin{matrix} g[j][k-B[v],ban[i]!=2 \\g[j][k],ban[i]!=3 \end{matrix} \right\} $

對于城市的轉(zhuǎn)移

$ f[i][j]=f[i-A[l]][j]+g[i][j] $

考慮如何合并答案

枚舉 $ f[i][j] $，則

$ sum-i-c1 \leq x \leq c0-i $，$ x $ 為 $ 1 $ 陣營的人數(shù)
$ sum-j-d1 \leq y \leq d0-j $，$ y $ 為 $ 1 $ 派系的人數(shù)

則沒有限制的方案數(shù)為 $ \sum_{x=sum-i-c1}^{c0-i}\times \sum_{y=sum-j-d1}^{d0-i} $，前綴和優(yōu)化即可 $ O(1) $ 詢問

時間效率：$ O(km^2) $，時間主要在 $ f $ 上，但是跑不滿

代碼

?代碼寫丑了，常數(shù)比較大

1 #include<bits/stdc++.h> 2 #define LL long long 3 #define pb push_back 4 #define clr(s) memset(s,0,sizeof s) 5 #define _(d) while(d(isdigit(ch=getchar()))) 6 using namespace std; 7 int R(){ 8 int x;bool f=1;char ch;_(!)if(ch=='-')f=0;x=ch^48; 9 _()x=(x<<3)+(x<<1)+(ch^48);return f?x:-x;} 10 const int N=3e3+5,P=998244353; 11 int n,m,K,g1[N],g2[N],f[N][N],h[N][N],c0,c1,d0,d1,ban[N],Ban[N],A[N],B[N],fa[N],sum,ans,s0,s1; 12 vector<int>p[N]; 13 void clean(){ 14 clr(g1),clr(g2),clr(f),clr(h),clr(A),clr(B),clr(fa); 15 for(int i=1;i<=m;i++)p[i].clear(); 16 sum=s0=s1=ans=0; 17 } 18 int make(int x,int y){ 19 int l0=max(0,sum-x-c1),r0=c0-x; 20 int l1=max(0,sum-y-d1),r1=d0-y; 21 if(l0>r0||l1>r1)return 0; 22 return 1ll*(g1[r0]-(l0?g1[l0-1]:0)+P)*(g2[r1]-(l1?g2[l1-1]:0)+P)%P; 23 } 24 int main(){ 25 for(int T=R();T--;clean()){ 26 memset(ban,-1,sizeof ban); 27 memset(Ban,-1,sizeof Ban); 28 n=R(),m=R(),c0=R(),c1=R(),d0=R(),d1=R(); 29 for(int i=1;i<=n;i++){ 30 fa[i]=R(),B[i]=R(); 31 sum+=B[i],A[fa[i]]+=B[i]; 32 p[fa[i]].pb(i); 33 } 34 K=R(); 35 for(int i=1;i<=K;i++){ 36 int id=R(),fl=R(); 37 ban[id]=Ban[fa[id]]=fl; 38 } 39 g1[0]=g2[0]=1; 40 for(int i=1;i<=m;i++) 41 if(!~Ban[i]&&A[i]) 42 for(int j=min(c0,sum);j>=A[i];j--) 43 g1[j]=(g1[j-A[i]]+g1[j])%P; 44 for(int i=1;i<=c0;i++)g1[i]=(g1[i]+g1[i-1])%P; 45 for(int i=1;i<=n;i++) 46 if(!~ban[i]) 47 for(int j=min(d0,sum);j>=B[i];j--) 48 g2[j]=(g2[j-B[i]]+g2[j])%P; 49 for(int i=1;i<=d0;i++)g2[i]=(g2[i]+g2[i-1])%P; 50 f[0][0]=1; 51 for(int i=1;i<=m;i++) 52 if(~Ban[i]&&A[i]){ 53 for(int x=min(c0,s0);~x;x--) 54 for(int y=min(d0,s1);~y;y--) 55 h[x][y]=f[x][y]; 56 for(int v:p[i])if(ban[v]!=-1){ 57 s1+=B[v]; 58 int f00=ban[v]^0?1:0,f01=ban[v]^1?1:0,f10=ban[v]^2?1:0,f11=ban[v]^3?1:0; 59 for(int x=min(c0,s0);~x;x--) 60 for(int y=min(d0,s1);~y;y--) 61 if(y>=B[v]){ 62 f[x][y]=(f[x][y]*f01+f[x][y-B[v]]*f00)%P; 63 h[x][y]=(h[x][y]*f11+h[x][y-B[v]]*f10)%P; 64 } 65 else f[x][y]=f[x][y]*f01,h[x][y]=h[x][y]*f11; 66 } 67 s0+=A[i]; 68 for(int x=min(c0,s0);~x;x--) 69 for(int y=min(c0,s1);~y;y--) 70 if(x>=A[i])f[x][y]=(f[x-A[i]][y]+h[x][y])%P; 71 else f[x][y]=h[x][y]; 72 } 73 for(int i=0;i<=min(s0,c0);i++) 74 for(int j=0;j<=min(s1,d0);j++) 75 ans=(ans+1ll*f[i][j]*make(i,j))%P; 76 cout<<ans<<endl; 77 } 78 return 0; 79 } View Code

?

轉(zhuǎn)載于:https://www.cnblogs.com/chmwt/p/10677849.html

總結(jié)

以上是生活随笔為你收集整理的「十二省联考 2019」皮配——dp的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。