Description

Input

Output

Sample Input

3

1 1 1

2 1 1

1 1 1

Sample Output

4

2 1

Data Constraint

對(duì)于30%的數(shù)據(jù)，N<=30

對(duì)于100%的數(shù)據(jù)，N<=80

Solution

• 這題挺熟悉，典型的最大匹配

• 可是它的第二問(wèn)是求最大匹配的交集，這使我一下子沒(méi)了思緒

• 后來(lái)我才猛然想到：這可以用費(fèi)用流來(lái)做！

• 先求完美匹配的值，于是設(shè)源點(diǎn)和匯點(diǎn)；

• 源點(diǎn)向每個(gè)男孩各連一條容量為 1 、費(fèi)用為 0 的邊；

• 每個(gè)女孩各向匯點(diǎn)連一條容量為 1 、費(fèi)用為 0 的邊；

• 最后每個(gè)男孩各向每個(gè)女孩連一條容量為 1 、費(fèi)用為 1 的邊。

• 之后跑一遍最大費(fèi)用最大流就可以啦！

• 而所謂的交集，就可以枚舉最大匹配的每一條邊，

• 將其刪除之后重新求最大匹配，若最大匹配發(fā)生變化則說(shuō)明這條邊在交集中！

• 注意：枚舉的邊必須是求解時(shí)使用過(guò)的邊，否則沒(méi)有必要枚舉。

Code

#include<cstdio> #include<cstring> using namespace std; const int N=200; struct data {int x,y; }a[N*N]; int n,t,ans,tot,num; int f[N][N],c[N][N],b[N][N]; int que[N*4],dis[N],pre[N]; bool bz[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; } inline bool spfa() {memset(dis,128,sizeof(dis));int l=que[1]=dis[0]=0,r=1;while(l<r){int now=que[++l];bz[now]=false;for(int i=0;i<=t;i++)if(f[now][i] && dis[now]+c[now][i]>dis[i]){dis[i]=dis[now]+c[now][i];pre[i]=now;if(!bz[i]) bz[que[++r]=i]=true;}}return dis[t]>0; } inline int work() {int sum=0;for(int i=t;i;i=pre[i]){f[pre[i]][i]--;f[i][pre[i]]++;sum+=c[pre[i]][i];}return sum; } int main() {n=read(),t=n*2+1;for(int i=1;i<=n;i++){for(int j=1;j<=n;j++){f[i][n+j]=1;c[i][n+j]=read();c[n+j][i]=-c[i][n+j];}f[0][i]=f[n+i][t]=1;}memcpy(b,f,sizeof(b));while(spfa()) ans+=work();printf("%d\n",ans);for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)if(!f[i][n+j]){a[++tot].x=i;a[tot].y=n+j;}for(int i=1;i<=tot;i++){memcpy(f,b,sizeof(f));int cnt=f[a[i].x][a[i].y]=0;while(spfa()) cnt+=work();if(cnt<ans) printf("%d %d\n",a[i].x,a[i].y-n);}return 0; }

總結(jié)

以上是生活随笔為你收集整理的JZOJ 3739. 【TJOI2014】匹配的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

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