AtCoder AGC029F Construction of a Tree (二分图匹配)

2024-02-15 15:32
文章标签 二分 匹配 tree atcoder construction agc029f

本文主要是介绍AtCoder AGC029F Construction of a Tree (二分图匹配)，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

题目链接

https://atcoder.jp/contests/agc029/tasks/agc029_f

题解

考虑如何才能构成一棵树：显然有一个必要条件是对于每个点\(u\)来说，整张图所有的边与除去\(u\)之外所有的点存在完美匹配（即考虑一张二分图左边是除了\(u\)之外点的集合右边是\(E_i\), \(u\)和\(E_i\)连边当且仅当\(u\in E_i\)，该图存在完美匹配）。用Hall定理来表达就是，设\(S\)为\(\{ E_1,E_2,..,E_n\}\)的任意一个子集，\(N(S)\)表示这些\(E_i\)对应的点集的并集，则\(|N(S)|\ge |S|+1\).
实际上这个条件也是充分条件。我们用一个构造算法来证明。直接从\(1\)号点开始BFS或者DFS，每次\(u\)选择一条出边到右边的一个点\(v\)，然后跳到右边的点的匹配点\(u'\). 若这两个点都没被访问过，则添加一条树边\((u,u')\). 上面的命题等价于这样搜索能够遍历所有的点。因为如果某一个时刻不能走到未走过的点，那么走过的左边点个数比右边点个数多\(1\)，左边点总个数比右边点总个数多\(1\)，故现在未遍历的右边的点的集合\(T\)满足\(N(T)\le T\)，这与上面的命题矛盾。而如果上面的命题不成立，显然无法搜出合法的方案。而这样搜索能遍历所有的点等价于原问题有解，故原命题等价于原问题有解。
时间复杂度\(O(\sum |E_i|\sqrt n)\).

注: 在这里由于时间复杂度的限制，我们只能用二分图匹配检验\(1\)号点是否满足去掉后有完美匹配，但是这并不代表所有点都有。我们证明了有解的充要条件是每个点去掉后都有完美匹配，也是\(1\)号点去掉后有完美匹配且BFS/DFS不会在遍历完所有点前终止，我们的算法是正确的。但是似乎数据里并没有这种\(1\)号点去掉后有完美匹配但实际上无解的情况，把下面BFS的代码中标注comment_1的那一行return 0;前面加一个assert(0);依然可以AC. 但实际上这种情况完全可能出现，Hack数据如下：

4
4 1 2 3 4
2 3 4
2 3 4

代码

BFS

#include<bits/stdc++.h>
#define llong long long
#define mkpr make_pair
#define riterator reverse_iterator
#define pii pair<int,int>
using namespace std;inline int read()
{int x = 0,f = 1; char ch = getchar();for(;!isdigit(ch);ch=getchar()) {if(ch=='-') f = -1;}for(; isdigit(ch);ch=getchar()) {x = x*10+ch-48;}return x*f;
}const int INF = 1e7;namespace NetFlow
{const int N = 2e5+2;const int M = 4e5;struct Edge{int v,w,nxt,rev;} e[(M<<1)+3];int fe[N+3];int te[N+3];int dep[N+3];int que[N+3];int n,en,s,t;void addedge(int u,int v,int w){en++; e[en].v = v; e[en].w = w;e[en].nxt = fe[u]; fe[u] = en; e[en].rev = en+1;en++; e[en].v = u; e[en].w = 0;e[en].nxt = fe[v]; fe[v] = en; e[en].rev = en-1;}bool bfs(){for(int i=1; i<=n; i++) dep[i] = 0;int head = 1,tail = 1; que[1] = s; dep[s] = 1;while(head<=tail){int u = que[head]; head++;for(int i=fe[u]; i; i=e[i].nxt){int v = e[i].v;if(e[i].w>0 && dep[v]==0){dep[v] = dep[u]+1;if(v==t) return true;tail++; que[tail] = v;}}}return false;}int dfs(int u,int cur){if(u==t||cur==0) {return cur;}int rst = cur;for(int &i=te[u]; i; i=e[i].nxt){int v = e[i].v;if(e[i].w>0 && rst>0 && dep[v]==dep[u]+1){int flow = dfs(v,min(rst,e[i].w));if(flow>0){e[i].w -= flow; rst -= flow;e[e[i].rev].w += flow;if(rst==0) {return cur;}}}}if(rst==cur) {dep[u] = -2;}return cur-rst;}int dinic(int _n,int _s,int _t){n = _n,s = _s,t = _t;int ret = 0;while(bfs()){for(int i=1; i<=n; i++) te[i] = fe[i];memcpy(te,fe,sizeof(int)*(n+1));ret += dfs(s,INF);}return ret;}
}
using NetFlow::addedge;
using NetFlow::dinic;const int N = 1e5;
vector<int> adj[(N<<1)+3];
vector<pair<int,pii> > ans;
int mch[(N<<1)+3];
bool vis[(N<<1)+3];
int que[N+3];
int n;bool bfs()
{int hd = 1,tl = 1; que[1] = 1; vis[1] = true;while(hd<=tl){int u = que[hd]; hd++;for(int o=0; o<adj[u].size(); o++){int v = adj[u][o]; if(vis[v]) continue;if(vis[mch[v]]) continue;que[++tl] = mch[v]; vis[v] = vis[mch[v]] = true;ans.push_back(mkpr(v,mkpr(u,mch[v])));}}if(tl<n) {return false;}return true;
}int main()
{scanf("%d",&n);for(int i=1; i<=n; i++) addedge(1,i+2,1);for(int i=n+1; i<n+n; i++) addedge(i+2,2,1);for(int i=1; i<n; i++){int sz; scanf("%d",&sz);while(sz--){int x; scanf("%d",&x); adj[i+n].push_back(x); adj[x].push_back(i+n);if(x!=1) {addedge(x+2,i+n+2,1);}}}if(dinic(n+n+1,1,2)<n-1) {puts("-1"); return 0;}for(int u=3; u<=n+2; u++){for(int i=NetFlow::fe[u]; i; i=NetFlow::e[i].nxt){int v = NetFlow::e[i].v; if(v<=n+2) continue;if(NetFlow::e[i].w==0){mch[u-2] = v-2,mch[v-2] = u-2;break;}}}
//  printf("match: "); for(int i=1; i<=n+n-1; i++) printf("%d ",mch[i]); puts("");if(!bfs()) {puts("-1"); return 0;} //comment_1sort(ans.begin(),ans.end());for(int i=0; i<ans.size(); i++) printf("%d %d\n",ans[i].second.first,ans[i].second.second);return 0;
}

DFS

#include<bits/stdc++.h>
#define llong long long
#define mkpr make_pair
#define riterator reverse_iterator
#define pii pair<int,int>
using namespace std;inline int read()
{int x = 0,f = 1; char ch = getchar();for(;!isdigit(ch);ch=getchar()) {if(ch=='-') f = -1;}for(; isdigit(ch);ch=getchar()) {x = x*10+ch-48;}return x*f;
}const int INF = 1e7;namespace NetFlow
{const int N = 2e5+2;const int M = 4e5;struct Edge{int v,w,nxt,rev;} e[(M<<1)+3];int fe[N+3];int te[N+3];int dep[N+3];int que[N+3];int n,en,s,t;void addedge(int u,int v,int w){en++; e[en].v = v; e[en].w = w;e[en].nxt = fe[u]; fe[u] = en; e[en].rev = en+1;en++; e[en].v = u; e[en].w = 0;e[en].nxt = fe[v]; fe[v] = en; e[en].rev = en-1;}bool bfs(){for(int i=1; i<=n; i++) dep[i] = 0;int head = 1,tail = 1; que[1] = s; dep[s] = 1;while(head<=tail){int u = que[head]; head++;for(int i=fe[u]; i; i=e[i].nxt){int v = e[i].v;if(e[i].w>0 && dep[v]==0){dep[v] = dep[u]+1;if(v==t) return true;tail++; que[tail] = v;}}}return false;}int dfs(int u,int cur){if(u==t||cur==0) {return cur;}int rst = cur;for(int &i=te[u]; i; i=e[i].nxt){int v = e[i].v;if(e[i].w>0 && rst>0 && dep[v]==dep[u]+1){int flow = dfs(v,min(rst,e[i].w));if(flow>0){e[i].w -= flow; rst -= flow;e[e[i].rev].w += flow;if(rst==0) {return cur;}}}}if(rst==cur) {dep[u] = -2;}return cur-rst;}int dinic(int _n,int _s,int _t){n = _n,s = _s,t = _t;int ret = 0;while(bfs()){for(int i=1; i<=n; i++) te[i] = fe[i];memcpy(te,fe,sizeof(int)*(n+1));ret += dfs(s,INF);}return ret;}
}
using NetFlow::addedge;
using NetFlow::dinic;const int N = 1e5;
vector<int> adj[(N<<1)+3];
vector<pair<int,pii> > ans;
int mch[(N<<1)+3];
bool vis[(N<<1)+3];
int n;void dfs(int u)
{for(int o=0; o<adj[u].size(); o++){int v = adj[u][o]; if(vis[v]) continue;if(vis[mch[v]]) continue;vis[v] = vis[mch[v]] = true; ans.push_back(mkpr(v,mkpr(u,mch[v])));dfs(mch[v]);}
}int main()
{scanf("%d",&n);for(int i=1; i<=n; i++) addedge(1,i+2,1);for(int i=n+1; i<n+n; i++) addedge(i+2,2,1);for(int i=1; i<n; i++){int sz; scanf("%d",&sz);while(sz--){int x; scanf("%d",&x); adj[i+n].push_back(x); adj[x].push_back(i+n);if(x!=1) {addedge(x+2,i+n+2,1);}}}if(dinic(n+n+1,1,2)<n-1) {puts("-1"); return 0;}for(int u=3; u<=n+2; u++){for(int i=NetFlow::fe[u]; i; i=NetFlow::e[i].nxt){int v = NetFlow::e[i].v; if(v<=n+2) continue;if(NetFlow::e[i].w==0){mch[u-2] = v-2,mch[v-2] = u-2;break;}}}
//  printf("match: "); for(int i=1; i<=n+n-1; i++) printf("%d ",mch[i]); puts("");vis[1] = true; dfs(1);if(ans.size()<n-1) {puts("-1"); return 0;}sort(ans.begin(),ans.end());for(int i=0; i<ans.size(); i++) printf("%d %d\n",ans[i].second.first,ans[i].second.second);return 0;
}

这篇关于AtCoder AGC029F Construction of a Tree (二分图匹配)的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！


http://www.chinasem.cn/article/711787。

相关文章

Nginx中location实现多条件匹配的方法详解

Nginx中location实现多条件匹配的方法详解

《Nginx中location实现多条件匹配的方法详解》在Nginx中,location指令用于匹配请求的URI,虽然location本身是基于单一匹配规则的,但可以通过多种方式实现多个条件的匹配逻辑... 目录1. 概述2. 实现多条件匹配的方式2.1 使用多个 location 块2.2 使用正则表达式
阅读更多...
golang字符串匹配算法解读

golang字符串匹配算法解读

《golang字符串匹配算法解读》文章介绍了字符串匹配算法的原理,特别是Knuth-Morris-Pratt（KMP）算法,该算法通过构建模式串的前缀表来减少匹配时的不必要的字符比较,从而提高效率,在... 目录简介KMP实现代码总结简介字符串匹配算法主要用于在一个较长的文本串中查找一个较短的字符串（称为
阅读更多...
C++使用栈实现括号匹配的代码详解

C++使用栈实现括号匹配的代码详解

《C++使用栈实现括号匹配的代码详解》在编程中,括号匹配是一个常见问题,尤其是在处理数学表达式、编译器解析等任务时,栈是一种非常适合处理此类问题的数据结构,能够精确地管理括号的匹配问题,本文将通过C+... 目录引言问题描述代码讲解代码解析栈的状态表示测试总结引言在编程中，括号匹配是一个常见问题，尤其是在
阅读更多...
关于Gateway路由匹配规则解读

关于Gateway路由匹配规则解读

《关于Gateway路由匹配规则解读》本文详细介绍了SpringCloudGateway的路由匹配规则,包括基本概念、常用属性、实际应用以及注意事项,路由匹配规则决定了请求如何被转发到目标服务,是Ga... 目录Gateway路由匹配规则一、基本概念二、常用属性三、实际应用四、注意事项总结Gateway路由
阅读更多...
hdu2241（二分+合并数组）

hdu2241（二分+合并数组）

题意：判断是否存在a+b+c = x，a，b，c分别属于集合A,B,C 如果用暴力会超时，所以这里用到了数组合并，将b,c数组合并成d，d数组存的是b，c数组元素的和，然后对d数组进行二分就可以了 代码如下（附注释）： #include<iostream>#include<algorithm>#include<cstring>#include<stack>#include<que
阅读更多...
hdu2289（简单二分）

hdu2289（简单二分）

虽说是简单二分，但是我还是wa死了  题意：已知圆台的体积，求高度 首先要知道圆台体积怎么求：设上下底的半径分别为r1,r2,高为h,V = PI*(r1*r1+r1*r2+r2*r2)*h/3 然后以h进行二分 代码如下： #include<iostream>#include<algorithm>#include<cstring>#include<stack>#includ
阅读更多...
【Prometheus】PromQL向量匹配实现不同标签的向量数据进行运算

【Prometheus】PromQL向量匹配实现不同标签的向量数据进行运算

✨✨ 欢迎大家来到景天科技苑✨✨ 🎈🎈 养成好习惯，先赞后看哦~🎈🎈 🏆 作者简介：景天科技苑 🏆《头衔》：大厂架构师，华为云开发者社区专家博主，阿里云开发者社区专家博主，CSDN全栈领域优质创作者，掘金优秀博主，51CTO博客专家等。 🏆《博客》：Python全栈，前后端开发，小程序开发，人工智能，js逆向，App逆向，网络系统安全，数据分析，Django，fastapi
阅读更多...
poj 2976 分数规划二分贪心（部分对总体的贡献度） poj 3111

poj 2976 分数规划二分贪心（部分对总体的贡献度） poj 3111

poj 2976： 题意： 在n场考试中，每场考试共有b题，答对的题目有a题。 允许去掉k场考试，求能达到的最高正确率是多少。 解析： 假设已知准确率为x，则每场考试对于准确率的贡献值为： a - b * x，将贡献值大的排序排在前面舍弃掉后k个。 然后二分x就行了。 代码： #include <iostream>#include <cstdio>#incl
阅读更多...
poj 3104 二分答案

poj 3104 二分答案

题意： n件湿度为num的衣服，每秒钟自己可以蒸发掉1个湿度。 然而如果使用了暖炉，每秒可以烧掉k个湿度，但不计算蒸发了。 现在问这么多的衣服，怎么烧事件最短。 解析： 二分答案咯。 代码： #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <c
阅读更多...
poj 3258 二分最小值最大

poj 3258 二分最小值最大

题意： 有一些石头排成一条线，第一个和最后一个不能去掉。 其余的共可以去掉m块，要使去掉后石头间距的最小值最大。 解析： 二分石头，最小值最大。 代码： #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <c
阅读更多...

Copyright China编程 23002807@qq.com

沪ICP备2023033072号-2