博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
hdu3081 Marriage Match II
阅读量:5330 次
发布时间:2019-06-14

本文共 3118 字,大约阅读时间需要 10 分钟。

新年第一篇,又花了一早上,真是蠢啊!

二分+网络流
之前对于讨论哪些人是朋友的时候复杂度过高
直接n3的暴力虽然看起来复杂度高,其实并不是每次都成立

#include
using namespace std;const int N = 205;const int INF = 0x3f3f3f3f;#define sz(X) ((int)X.size())int A[N], B[N];int mp[N][N];int f[N];int find(int x) { return f[x] == x? x : f[x] = find(f[x]); }/***********Dinic*********/int s, t;struct Edge { int from, to, cap, flow, nx; Edge(int a=0, int b=0, int c=0, int d=0, int e=0):from(a), to(b), cap(c), flow(d),nx(e){}}E[N*N];int head[N], tot;bool vis[N];int d[N];int cur[N];void init() { tot = 0; memset(head,-1,sizeof(head));}void add(int from, int to, int cap) { E[tot] = Edge(from,to,cap,0,head[from]); head[from] = tot++; E[tot] = Edge(to,from,0,0,head[to]); head[to] = tot++;}bool bfs() { memset(vis, 0, sizeof(vis)); queue
Q; Q.push(s); d[s] = 0; vis[s] = 1; while(!Q.empty()) { int x = Q.front(); Q.pop(); for(int i = head[x]; ~i; i = E[i].nx) { Edge& e = E[i]; if(!vis[e.to] && e.cap > e.flow) { vis[e.to] = 1; d[e.to] = d[x] + 1; Q.push(e.to); } } } return vis[t];}int dfs(int x, int a) { if(x == t || a == 0) return a; int flow = 0, f; for(int& i = cur[x]; ~i; i = E[i].nx) { Edge& e = E[i]; if(d[x]+1 == d[e.to] && (f = dfs(e.to, min(a, e.cap-e.flow))) > 0) { e.flow += f; E[i^1].flow -= f; flow += f; a -= f; if(a == 0) break; } } return flow;}int Maxflow(int start, int ed) { s = start; t = ed; int flow = 0; while(bfs()) { for(int i = s; i <= t; ++i) cur[i] = head[i]; flow += dfs(s, INF); } return flow;}int main() { int n,m,l; int _; scanf("%d",&_); while(_--) { memset(mp,0,sizeof(mp)); scanf("%d %d %d",&n,&m,&l); for(int i = 0; i <= n; ++i) f[i] = i; for(int i = 0; i < m; ++i) { int a, b; scanf("%d %d",&a,&b); mp[a][b] = 1; } for(int i = 0; i < l; ++i) { int c,d; scanf("%d %d",&c,&d); int t1 = find(c); int t2 = find(d); if(t1 != t2) f[t1] = t2; } for(int i = 1; i <= n; ++i) { for(int j = 1; j <= n; ++j) { if(find(i) == find(j)) { for(int k = 1; k <= n; ++k) { if(mp[i][k]) mp[j][k] = 1; } } } } int s = 0; int t = 2*n+1; int l = 0; int r = n; int res; while(l <= r) { int mid = (l + r)>>1; init(); for(int i = 1; i <= n; ++i) add(s,i,mid); for(int i = 1; i <= n; ++i) add(i+n,t,mid); for(int i = 1; i <= n; ++i) for(int j = 1; j <= n; ++j) { if(mp[i][j]) { add(i,j+n, 1); } } int tt = Maxflow(s, t); if(tt != mid*n) r = mid-1; else { l = mid+1; res = mid; } } printf("%d\n",res); } return 0;}

转载于:https://www.cnblogs.com/Basasuya/p/8433719.html

你可能感兴趣的文章
【转】redo与undo
查看>>
wpf样式绑定 行为绑定 事件关联 路由事件实例
查看>>
Oracle事务
查看>>
String类中的equals方法总结(转载)
查看>>
内存地址对齐
查看>>
创新课程管理系统数据库设计心得
查看>>
Could not resolve view with name '***' in servlet with name 'dispatcher'
查看>>
[转载] redis 的两种持久化方式及原理
查看>>
管道,数据共享,进程池
查看>>
SDUTOJ3754_黑白棋(纯模拟)
查看>>
php中的isset和empty的用法区别
查看>>
把word文档中的所有图片导出
查看>>
ubuntu 18.04取消自动锁屏以及设置键盘快捷锁屏
查看>>
arcgis api 4.x for js 结合 Echarts4 实现散点图效果(附源码下载)
查看>>
YTU 2625: B 构造函数和析构函数
查看>>
apache自带压力测试工具ab的使用及解析
查看>>
加固linux
查看>>
Hyper-V虚拟机上安装一个图形界面的Linux系统
查看>>
js千分位处理
查看>>
字符串类型的相互转换
查看>>