P3731 题解
简要题意是找到一条边连接使得最大团大小增加。
在补图上最大团等于最大独立集。
所以问题转化为删掉一条边使得最大独立集增加,又因为团不超过两个,所以原图是二分图,也就是使得最大匹配减少。
考虑什么样的匹配边是可以被代替的,先跑一遍网络流求最大匹配,不难发现假若一条匹配在残量网络流上的一个环中,那么它就是可以被代替的,所以对残量网络进行缩点,两个端点不在同一个强连通分量中的匹配边可以作为答案。
#include <bits/stdc++.h>
using namespace std;
const int maxn = 1e4 + 114;
const int maxm = 2e6 + 114;
const int inf = 0x3f3f3f3f;
int maxflow, vis[maxn], tot = 1, cnt;
int s, t, n, m, Q, E;
vector<int> edge[maxn], Road[maxn];
int hd[maxn], road[maxn], dis[maxn];
map<int, int> f[maxn];
struct edge {
int next, to, w;
} e[maxm * 2];
void add(int u, int v, int w) {
e[++tot].to = v;
f[u][v] = tot;
e[tot].w = w;
e[tot].next = hd[u];
hd[u] = tot;
e[++tot].to = u;
e[tot].w = 0;
e[tot].next = hd[v];
hd[v] = tot;
}
bool bfs() {
memset(dis, 0, sizeof(dis));
dis[s] = 1;
queue<int>Q;
Q.push(s);
while (!Q.empty()) {
int u = Q.front();
Q.pop();
road[u] = hd[u];
for (int i = hd[u]; i; i = e[i].next) {
int v = e[i].to;
if (!dis[v] && e[i].w) {
dis[v] = dis[u] + 1;
Q.push(v);
}
}
}
return dis[t] != 0;
}
int dinic(int now, int res) {
if (now == t)
return res;
int tp = res;
for (int i = road[now]; i; i = e[i].next) {
int v = e[i].to;
road[now] = i;
if (dis[v] == dis[now] + 1 && e[i].w) {
int k = min(e[i].w, tp);
int del = dinic(v, k);
e[i].w -= del;
e[i ^ 1].w += del;
tp -= del;
if (!tp)
break;
}
}
return res - tp;
}
int col[maxn], Use[maxn];
void dfs(int u) {
if (Use[u] == 1)
return ;
Use[u] = 1;
for (int v : edge[u])
col[v] = col[u] ^ 1, dfs(v);
}
int dfsn[maxn];
int low[maxn];
stack<int> S;
int Vis[maxn], use[maxn];
int color[maxn], sum = 0;
int deep = 0;
void paint(int u) {
S.pop();
color[u] = sum;
Vis[u] = 0;
}
void tanjan(int u) {
dfsn[u] = ++deep;
low[u] = deep;
Vis[u] = 1;
use[u] = 1;
S.push(u);
for (int i = hd[u]; i; i = e[i].next) {
if (e[i].w == 0)
continue;
int v = e[i].to;
if (dfsn[v] == 0) {
tanjan(v);
low[u] = min(low[u], low[v]);
} else {
if (Vis[v] != 0) {
low[u] = min(low[u], low[v]);
}
}
}
if (dfsn[u] == low[u]) {
sum++;
while (S.top() != u) {
paint(S.top());
}
paint(u);
}
return ;
}
set< pair<int, int>> chifan;
int main() {
cin >> n >> m;
for (int i = 1; i <= m; i++) {
int u, v;
cin >> u >> v;
edge[u].push_back(v);
edge[v].push_back(u);
}
for (int i = 1; i <= n; i++) {
if (Use[i] == 0)
col[i] = 0, dfs(i);
}
for (int i = 1; i <= n; i++) {
if (col[i] == 0) {
for (int nxt : edge[i])
Road[i].push_back(nxt), add(i, nxt, 1);
}
}
s = n + 1, t = n + 2;
for (int i = 1; i <= n; i++)
if (col[i] == 0)
add(s, i, 1);
else
add(i, t, 1);
while (bfs() == true)
maxflow += dinic(s, inf);
for (int i = 1; i <= n + 2; i++) {
if (use[i] == 0)
tanjan(i);
}
for (int i = 1; i <= n; i++) {
for (int nxt : Road[i]) {
if (e[f[i][nxt]].w == 0 && color[i] != color[nxt]) {
chifan.insert(make_pair(min(i, nxt), max(i, nxt)));
}
}
}
cout << chifan.size() << '\n';
for (pair<int, int> OUT : chifan) {
cout << OUT.first << ' ' << OUT.second << '\n';
}
}
P3731 题解的更多相关文章
- 2016 华南师大ACM校赛 SCNUCPC 非官方题解
我要举报本次校赛出题人的消极出题!!! 官方题解请戳:http://3.scnuacm2015.sinaapp.com/?p=89(其实就是一堆代码没有题解) A. 树链剖分数据结构板题 题目大意:我 ...
- noip2016十连测题解
以下代码为了阅读方便,省去以下头文件: #include <iostream> #include <stdio.h> #include <math.h> #incl ...
- BZOJ-2561-最小生成树 题解(最小割)
2561: 最小生成树(题解) Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1628 Solved: 786 传送门:http://www.lyd ...
- Codeforces Round #353 (Div. 2) ABCDE 题解 python
Problems # Name A Infinite Sequence standard input/output 1 s, 256 MB x3509 B Restoring P ...
- 哈尔滨理工大学ACM全国邀请赛(网络同步赛)题解
题目链接 提交连接:http://acm-software.hrbust.edu.cn/problemset.php?page=5 1470-1482 只做出来四道比较水的题目,还需要加强中等题的训练 ...
- 2016ACM青岛区域赛题解
A.Relic Discovery_hdu5982 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Jav ...
- poj1399 hoj1037 Direct Visibility 题解 (宽搜)
http://poj.org/problem?id=1399 http://acm.hit.edu.cn/hoj/problem/view?id=1037 题意: 在一个最多200*200的minec ...
- 网络流n题 题解
学会了网络流,就经常闲的没事儿刷网络流--于是乎来一发题解. 1. COGS2093 花园的守护之神 题意:给定一个带权无向图,问至少删除多少条边才能使得s-t最短路的长度变长. 用Dijkstra或 ...
- CF100965C题解..
求方程 \[ \begin{array}\\ \sum_{i=1}^n x_i & \equiv & a_1 \pmod{p} \\ \sum_{i=1}^n x_i^2 & ...
- JSOI2016R3 瞎BB题解
题意请看absi大爷的blog http://absi2011.is-programmer.com/posts/200920.html http://absi2011.is-programmer.co ...
随机推荐
- CMS垃圾收集器小实验之CMSInitiatingOccupancyFraction参数
CMS垃圾收集器小实验之CMSInitiatingOccupancyFraction参数 背景 测试CMSInitiatingOccupancyFraction参数,测试结果和我的预期不符,所以花了一 ...
- Golang 爬虫01
目录 学习地址: 目录站: 爬虫概念: 工作流程: 百度贴吧爬虫实现: go实战代码 单进程 实现过程: 并发爬取 实现过程: 学习地址: https://www.bilibili.com/video ...
- console小知识
console.log(JSON.stringify(object,null,2));
- paramiko连接windows10详解,远程管理windows服务器
1.win10安装 OpenSSH 官网链接:https://docs.microsoft.com/zh-cn/windows-server/administration/openssh/openss ...
- sass 导入@import详解
@import Sass 拓展了 @import 的功能,允许其导入 SCSS 或 Sass 文件.被导入的文件将合并编译到同一个 CSS 文件中,另外,被导入的文件中所包含的变量或者混合指令 ( ...
- Windows平台git clone文件路径太长报错
问题描述 在Windows下拉取一些比较大的开源项目经常会提示文件路径太长(filename too long),然后死活都不成功 解决办法 1.配置git git config --system c ...
- wordpress博客系统报错
第一种,只显示nginx的默认网页 说明wordpress的网页配置文件没有被系统读取 我们就需要去查看nginx的配置文件/etc/nginx/conf.d/default.conf 首先,查看是不 ...
- proteus 器件名称被软件篡改bug的解决方案
proteus v7.8 器件名称被软件篡改bug 的解决方案 BUG描述 在做单片机实验时,发现从某一个时间保存的设计图文件开始,在添加新的电子元件时会出现部分旧元件的名称被捆绑替换为新元件的名称, ...
- webpack配置css预处理
webpack默认只支持js的打包,不支持其它类型,为了让它支持样式的打包就需要加载一些loader 打包css文件 在webpack中配置对应的loader 在入口js文件中通过import导入样式 ...
- 基于 Easysearch kNN 搭建即时图片搜索服务
环境准备 启动 Easysearch 服务: # Make sure your vm.max_map_count meets the requirement sudo sysctl -w vm.max ...