[十二省聯考2019]皮配
阿新 • • 發佈:2019-04-18
() 影響 https 相同 sum 容量 鏈接 define fine
題目鏈接
選一個派系和一個陣營可以唯一確定一名導師
因為每一個陣營裏的導師都分別來自不同派系,所以k=0時,對陣營的選擇是不影響對派系的選擇的
唯一的限制就是同城市的要在同一個陣營
所以以每個城市為物品,物品大小為該城市的人數,陣營人數為背包容量,做背包dp
再以每個學校為物品,物品大小為該學校的人數,派系人數為背包容量,做背包dp
只用一維記錄背包大小即可,因為總人數-背包裏的人數=在另一個陣營或派系的人數
然後合並答案即可
方案數是可以相互乘起來的,k很小,所以我們可以暴力做k!=0的情況,然後乘上符合要求的k==0的方案數
k!=0時,記\(f[x][t][i][j]\)為前x個學校,前一個學校選擇了t陣營,此時藍有i個人,鴨派有j個人的方案數
滾動第一維,否則空間會爆
將學校按城市排序,這樣相同城市的就會排在一起,轉移的時候如果和前一個學校同城就要選擇相同陣營
#include<bits/stdc++.h> #define rep(i,j,k) for(int i=j;i<=k;++i) using namespace std; typedef long long ll; typedef double db; char cch; inline int rd(){ int x=0,fl=1; cch=getchar(); while(cch>'9'||cch<'0'){ if(cch=='-') fl=-1; cch=getchar(); } while(cch>='0'&&cch<='9') x=(x<<3)+(x<<1)+cch-'0',cch=getchar(); return x*fl; } const int mod=998244353,N=3000; struct abc{ int ct,sum,ban; }p1[N],p2[N]; int f[N],g[N],ff[2][2][N][N],ok[N],bl[N],ct[N],city[N],sum[N],ban[N]; inline void inc(int &A,int B){//會比 %mod快一點點 A+=B; if(A>mod) A-=mod; if(A<0) A+=mod; } inline int mul(int A,int B){ return 1ll*A*B%mod; } inline int sub(int a,int b){ a-=b; if(a<0) a+=mod; return a; } inline int getg(int l,int r){ if(l>r) return 0; if(l==0) return g[r];//為了dp方便g[0]=1,實際上應該是0 return sub(g[r],g[l-1]); } inline int getf(int l,int r){ if(l>r) return 0; if(l==0) return f[r];//同理 return sub(f[r],f[l-1]); } inline bool cmp(const abc &a1,const abc &a2){ return a1.ct<a2.ct; } inline void sv(){ int n=rd(),c=rd(),c0=rd(),c1=rd(),d0=rd(),d1=rd(),ans=0,tot=0;//n所學校,c個城市,藍陣營 C0。 紅陣營 C1。 鴨派系 D0。 R派系D1。 memset(city,0,sizeof city); rep(i,1,n) ct[i]=rd(),sum[i]=rd(),tot+=sum[i],ban[i]=-1,city[ct[i]]+=sum[i];//city[i]表示第i個城市一共有多少人 int k=rd(),id; rep(i,1,k) id=rd(),ban[id]=rd(); int len1=0,len2=0; rep(i,1,n){ if(ban[i]!=-1) p1[++len1]=(abc){ct[i],sum[i],ban[i]};//有特殊要求的 else p2[++len2]=(abc){ct[i],sum[i],ban[i]};//沒有特殊要求的 } sort(p1+1,p1+len1+1,cmp);//按城市排序 rep(i,1,len1){ if(city[p1[i].ct]==-1) ok[i]=0;//陣營的轉移以城市為單位 else ok[i]=city[p1[i].ct]/*註意*/,city[p1[i].ct]=-1; } // memset(g,0,sizeof g),g[0]=1; rep(i,1,c) if(city[i]>0) for(int j=c0;j>=city[i];--j) inc(g[j],g[j-city[i]]); /*做前綴和*/rep(i,1,c0) inc(g[i],g[i-1]); // memset(f,0,sizeof f),f[0]=1; rep(i,1,len2) for(int j=d0;j>=p2[i].sum;--j) inc(f[j],f[j-p2[i].sum]); /*做前綴和*/rep(i,1,d0) inc(f[i],f[i-1]); // memset(ff,0,sizeof ff);//ff[x][t][i][j]為前x個學校,前一個學校選擇了t陣營,此時藍有i個人,鴨派有j個人的方案數,滾動第一維 ff[0][0][0][0]=1; int cnt=0,now=0; rep(i,1,len1){//對有要求的學校暴力求解 now^=1; int tmp=p1[i].sum,bn=p1[i].ban,d=ok[i],lst=cnt;//lst=之前的學校的總人數 cnt+=tmp; rep(t,0,1) rep(h,0,c0) rep(j,0,cnt) ff[now][t][h][j]=0;//這裏不可以用memset,用了會超時,因為一開始cnt很小,所以循環更快 rep(t,0,1){ int cs=-1;//choise if(i>1&&p1[i].ct==p1[i-1].ct) cs=t; for(int i=c0;i>=0;--i) for(int j=cnt;j>=0;--j){ if(cs!=1){//如果同城市的選擇了0陣營,或與上一個不同城, if(bn!=1&&i>=d&&j<=lst) inc(ff[now][0][i][j],ff[now^1][t]/*註意是t而不是0*/[i-d][j]);//沒有禁掉小R,可以加入R派 if(bn!=0&&i>=d&&j-tmp<=lst&&j>=tmp) inc(ff[now][0][i][j],ff[now^1][t][i-d][j-tmp]);//沒有禁掉Yazid,可以加入鴨派 } if(cs!=0){//如果同城市的選擇了1陣營,或與上一個不同城, if(bn!=3&&j<=lst) inc(ff[now][1][i][j],ff[now^1][t][i][j]); if(bn!=2&&j-tmp<=lst&&j>=tmp) inc(ff[now][1][i][j],ff[now^1][t][i][j-tmp]); } } } } // rep(t,0,1) rep(i,0,c0) rep(j,0,d0){ int v=ff[now][t][i][j]; if(!v)continue; int t1=c0-i,t2=max(0,tot-d1-j),t3=max(0,tot-c1-i),t4=d0-j;//符合人數要求的區間 inc(ans,mul(v,mul(getg(t3,t1),getf(t2,t4)))); } printf("%d\n",ans); } int main(){ int T=rd(); while(T--) sv(); } /* 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 */
[十二省聯考2019]皮配