日韩av黄I国产麻豆传媒I国产91av视频在线观看I日韩一区二区三区在线看I美女国产在线I麻豆视频国产在线观看I成人黄色短片

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 >

P6113-[模板]一般图最大匹配【带花树】

發布時間:2023/12/3 66 豆豆
生活随笔 收集整理的這篇文章主要介紹了 P6113-[模板]一般图最大匹配【带花树】 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

正題

題目鏈接:https://www.luogu.com.cn/problem/P6113


題目大意

給出一張無向圖,求最大匹配。

1≤n≤103,1≤m≤5×1041\leq n\leq 10^3,1\leq m\leq 5\times 10^41n103,1m5×104


解題思路

帶花樹的模板,我也不會講/kel

所以看下面兩篇大佬的博客吧
yyb-帶花樹算法學習筆記
Bill Yang-帶花樹學習筆記

時間復雜度好像是O(n3)O(n^3)O(n3)


code

#include<cstdio> #include<cstring> #include<algorithm> #include<queue> using namespace std; const int N=1010,M=3e5+10; struct node{int to,next; }a[M]; int n,m,tot,tk,ans,ls[N],dfn[N]; int fa[N],pre[N],tag[N],match[N]; queue<int> q; void addl(int x,int y){a[++tot].to=y;a[tot].next=ls[x];ls[x]=tot;return; } int find(int x) {return (fa[x]==x)?x:(fa[x]=find(fa[x]));} int LCA(int x,int y){++tk;x=find(x);y=find(y);while(dfn[x]!=tk){dfn[x]=tk;x=find(pre[match[x]]);if(y)swap(x,y);}return x; } void Blossom(int x,int y,int lca){while(find(x)!=lca){pre[x]=y;y=match[x];if(tag[y]==2){tag[y]=1;q.push(y);}fa[x]=fa[y]=lca;x=pre[y];}return; } int Aug(int s){memset(tag,0,sizeof(tag));memset(pre,0,sizeof(pre));for(int i=1;i<=n;i++)fa[i]=i;while(!q.empty())q.pop();q.push(s);tag[s]=1;while(!q.empty()){int x=q.front();q.pop();for(int i=ls[x];i;i=a[i].next){int y=a[i].to;if(!tag[y]){tag[y]=2;pre[y]=x;if(!match[y]){for(int u=y,lst;u;u=lst)lst=match[pre[u]],match[u]=pre[u],match[pre[u]]=u;return 1;}tag[match[y]]=1;q.push(match[y]);}else if(tag[y]==1&&find(y)!=find(x)){int lca=LCA(x,y);Blossom(x,y,lca);Blossom(y,x,lca);}}}return 0; } int main() {scanf("%d%d",&n,&m);for(int i=1;i<=m;i++){int x,y;scanf("%d%d",&x,&y);addl(x,y);addl(y,x);}for(int i=1;i<=n;i++)if(!match[i])ans+=Aug(i);printf("%d\n",ans);for(int i=1;i<=n;i++)printf("%d ",match[i]);return 0; }

總結

以上是生活随笔為你收集整理的P6113-[模板]一般图最大匹配【带花树】的全部內容,希望文章能夠幫你解決所遇到的問題。

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