日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

P3731 二分图匹配必经边

發布時間:2024/4/15 编程问答 29 豆豆
生活随笔 收集整理的這篇文章主要介紹了 P3731 二分图匹配必经边 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

題意經過一番轉換變成了 讓你在一個二分圖上刪一條邊使得二分圖的最大獨立集大小至少+1

二分圖的最大獨立集=點數-最小點覆蓋(最大匹配) 點數是固定不變的 所以我們要減少最大匹配數

則刪掉的哪一條邊必須是二分圖匹配里必須的一條邊

做法:先Dinic跑一次 然后Tarjan縮點 找到圖中端點不為S/T且滿流的邊 再判兩端點是否屬于一個SCC

如果屬于一個SCC則必定存在另一個邊可以替換掉這條匹配邊

#include<iostream> #include<cstdio> #include<cstdlib> #include<cstring> #include<cmath> #include<algorithm> #include<vector> #include<queue> using namespace std; #define MAXN 10100 #define MAXM 150150 #define inf 1e9 inline int read() {int x=0;bool t=false;char ch=getchar();while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();if(ch=='-')t=true,ch=getchar();while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();return t?-x:x; } struct Line{int v,next,w;}e[MAXM<<1]; int h[MAXN],cnt=2; inline void Add(int u,int v,int w) {e[cnt]=(Line){v,h[u],w};h[u]=cnt++;e[cnt]=(Line){u,h[v],0};h[v]=cnt++; } int S,T,level[MAXN]; bool bfs() {for(int i=S;i<=T;++i)level[i]=0;queue<int> Q;level[S]=1;Q.push(S);while(!Q.empty()){int u=Q.front();Q.pop();for(int i=h[u];i;i=e[i].next)if(e[i].w&&!level[e[i].v])level[e[i].v]=level[u]+1,Q.push(e[i].v);}return level[T]; } int cur[MAXN]; int dfs(int u,int flow) {if(u==T||!flow)return flow;int ret=0;for(int &i=cur[u];i;i=e[i].next){int v=e[i].v,d;if(e[i].w&&level[v]==level[u]+1){d=dfs(v,min(flow,e[i].w));ret+=d;flow-=d;e[i].w-=d;e[i^1].w+=d;if(!flow)break;}}if(!ret)level[u]=0;return ret; } int Dinic() {int ret=0;while(bfs()){for(int i=S;i<=T;++i)cur[i]=h[i];ret+=dfs(S,inf);}return ret; } vector<pair<int,int> > Ans; vector<int> E[MAXN]; int n,m,col[MAXN]; void pre(int u,int c){col[u]=c;for(int v:E[u])if(!col[v])pre(v,c^1);} int dfn[MAXN],low[MAXN],tim,G[MAXN],gr; int St[MAXN],top;bool ins[MAXN]; void Tarjan(int u) {dfn[u]=low[u]=++tim;St[++top]=u;ins[u]=true;for(int i=h[u];i;i=e[i].next){int v=e[i].v;if(!e[i].w)continue;if(!dfn[v])Tarjan(v),low[u]=min(low[u],low[v]);else if(ins[v])low[u]=min(low[u],dfn[v]);}if(dfn[u]==low[u]){++gr;int v;do{v=St[top--];G[v]=gr;ins[v]=false;}while(u!=v);} } int main() {n=read();m=read();for(int i=1;i<=m;++i){int u=read(),v=read();E[u].push_back(v);E[v].push_back(u);}for(int i=1;i<=n;++i)if(!col[i])pre(i,2);S=0;T=n+1;for(int i=1;i<=n;++i)if(col[i]==2){Add(S,i,1);for(int v:E[i])Add(i,v,1);}else Add(i,T,1);int ret=Dinic();for(int i=S;i<=T;++i)if(!dfn[i])Tarjan(i);for(int u=1;u<=n;++u)if(col[u]==2)for(int i=h[u];i;i=e[i].next){int v=e[i].v;if(e[i].w)continue;if(v==S||v==T)continue;if(G[u]!=G[v])Ans.push_back(u<v?make_pair(u,v):make_pair(v,u));}sort(Ans.begin(),Ans.end());printf("%d\n",(int)Ans.size());for(auto a:Ans)printf("%d %d\n",a.first,a.second);return 0; }

?

轉載于:https://www.cnblogs.com/Aragaki/p/10605852.html

總結

以上是生活随笔為你收集整理的P3731 二分图匹配必经边的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。