http://uoj.ac/problem/168

没想到是网络流

官方题解地址

http://jiry-2.blog.uoj.ac/blog/1115

subtask2告诉我们度数为012的点对答案无影响

subtask3告诉我们原图\(|E| > 2|V| - 2\)时不是丛林的

证明一个结论

若对于原图所有的子图都满足\(|E| \le 2|V| - 2\)则是一个丛林

对子图大小施归纳法

\(n = 1\)是丛林

当\(n \geq 2\)时

设当前图度数最小的节点为\(u\)

仅讨论\(du[u] = 3\)的情况(0, 1, 2同subtask2)

记\(|E| = 2|V| - 2\)的子图为Dark图

Dark 交 Dark 得到的图也是Dark

证明的话考虑总的减去交的 每一部分都是满足的

然后设\(u\)在原图和\(a, b, c\)有边

在原图必不存在包含\(a, b, c\)的Dark子图(反证法)

然后就可以科学的构造\(n\)的时候的丛林了

如果有一个Dark子图包含\(a,b\)不包含\(c\)

那么一颗树选择\(u->a, u->b\)另一棵树选择\(u->c\)

然后就是问是否存在一个非空子图

满足\(|E| > 2|V| - 2\)

相当于点权值为\(-2\)边权值为\(1\)

这是一个最大权闭合子图的模型

但是是非空的

那么就枚举必选哪一个点

把这个点权值设成\(0\)跑网络流就好了

复杂度是\(\mathcal O(nm\sqrt n)\)的

观察到每一次枚举一个点只会改变一条边和他的反向边流量

所以退流就好了

学了一下退流的姿势

比如说要退边\(<u, v>\)的流

只需要\(t->v\)和\(u->s\)分别跑最大流就好了

感性YY一下

可能讲的好乱 看官方题解吧

还要卡常数qwq

#include <bits/stdc++.h>
#define rint register int
#define fo(i, n) for(int i = 1; i <= (n); i ++)
#define out(x) cerr << #x << " = " << x << "\n"
#define type(x) __typeof((x).begin())
#define foreach(it, x) for(type(x) it = (x).begin(); it != (x).end(); ++ it)
using namespace std;
// by piano
template<typename tp> inline void read(tp &x) {
x = 0;char c = getchar(); bool f = 0;
for(; c < '0' || c > '9'; f |= (c == '-'), c = getchar());
for(; c >= '0' && c <= '9'; x = (x << 3) + (x << 1) + c - '0', c = getchar());
if(f) x = -x;
}
namespace mf {
static const int N = 8888;
static const int M = 23333;
static const int inf = 1e9;
struct E {
int nxt, to, fl;
}e[M << 1];
int d[N], head[N], cur[N], e_cnt = 1;
int n, s, t;
inline void init(int _n, int _s, int _t) {
n = _n; e_cnt = 1;
fill(head, head + n + 1, 0);
}
inline void add(int u, int v, int fl) {
e[++ e_cnt] = (E) {head[u], v, fl}; head[u] = e_cnt;
}
inline void link(int u, int v, int fl) {
add(u, v, fl); add(v, u, 0);
}
inline bool bfs(void) {
queue<int> q;
fill(d, d + n + 1, -1);
d[s] = 0; q.push(s);
while(!q.empty()) {
int u = q.front(); q.pop();
for(rint i = head[u]; i; i = e[i].nxt) if(e[i].fl) {
int v = e[i].to;
if(!~ d[v]) {
d[v] = d[u] + 1, q.push(v);
if(v == t) return true;
}
}
}
return d[t] != -1;
}
inline int dfs(int u, int f) {
if(u == t || !f) return f;
rint used = 0, t;
for(rint &i = cur[u]; i; i = e[i].nxt) if(e[i].fl) {
int v = e[i].to;
if(d[v] != d[u] + 1) continue;
t = dfs(v, min(f - used, e[i].fl));
used += t; e[i].fl -= t; e[i ^ 1].fl += t;
if(used == f) break;
}
if(!used) d[u] = -1;
return used;
}
inline int doit(int _s, int _t, int f) {
int fl = 0; s = _s; t = _t;
while(bfs() && fl < f) {
memcpy(cur, head, sizeof (int) * (n + 2));
fl += dfs(s, f - fl);
}
return fl;
}
} const int inf = 1e9;
int n, m, s, t;
#define GG return puts("No"), 0;
main(void) {
read(n); read(m);
s = n + m + 1; t = s + 1;
mf::init(t, s, t);
for(rint i = 1; i <= n; i ++) {
mf::link(i, t, 2);
}
for(rint i = 1; i <= m; i ++) {
mf::link(s, i + n, 1);
}
for(rint i = 1; i <= m; i ++) {
int x, y; read(x); read(y);
mf::link(i + n, x, inf);
mf::link(i + n, y, inf);
}
int fl = mf::doit(s, t, inf);
if(m - fl > 0) GG;
for(rint i = 1; i <= n; i ++) {
fl -= mf::doit(i, s, mf::e[i * 2 + 1].fl);
mf::e[i * 2].fl = mf::e[i * 2 + 1].fl = 0;
fl += mf::doit(s, t, inf);
if(m - fl > 0) GG;
mf::e[i * 2].fl = 2;
}
cout << "Yes\n";
}

网路流 uoj 168 元旦老人与丛林的更多相关文章

  1. 2018.07.28 uoj#169. 【UR #11】元旦老人与数列(线段树)

    传送门 线段树好题. 维护区间加,区间取最大值,维护区间最小值,历史区间最小值. 同样先考虑不用维护历史区间最小值的情况,这个可以参考这道题的解法,维护区间最小和次小值可以解决前两个操作,然后使用历史 ...

  2. 【UOJ#169】元旦老人与数列

    论文题. 考虑到这题的维护和区间操作是反向的,也就是说无法像V那题快速的合并标记. 我们知道,一个区间的最小值和其他值是可以分开来维护的,因为如果一个区间被整体覆盖,那么最小值始终是最小值. 对于被覆 ...

  3. [UOJ #167]【UR #11】元旦老人与汉诺塔

    题目大意:给你一个有$n$个盘子的汉诺塔状态$S$,问有多少种不同的操作方法,使得可以在$m$步以内到达状态$T$.$n,m\leqslant100$ 题解:首先可以知道的是,一个状态最多可以转移到其 ...

  4. UOJ 52 元旦激光炮

    http://uoj.ac/problem/52 题意:每次可以得到3个序列中 思路:每次分别取出三个序列的K/3长度的位置,取最小的那个,然后每次减掉它,总复杂度是Nlog3N #include & ...

  5. UR11 A.元旦老人与汉诺塔

    题目:http://uoj.ac/contest/23/problem/167 如果我们拿个map来存状态的话.设当前状态是v,下一个状态是s.有f[i+1][s]+=f[i][v]. 初始f[0][ ...

  6. uoj169:元旦老人与数列

    题意:http://uoj.ac/problem/169 sol  :线段树..........蜜汁TLE了一个点,不管了..... 代码抄snowMyDream的,orz........... 线段 ...

  7. uoj167 元旦老人与汉诺塔(记忆化搜索)

    QwQ太懒了,题目直接复制uoj的了 QwQ这个题可以说是十分玄学的一道题了 首先可以暴搜,就是\(dfs\)然后模拟每个过程是哪个柱子向哪个柱子移动 不多解释了,不过实现起来还是有一点点难度的 直接 ...

  8. UOJ169. 【UR #11】元旦老人与数列

    传送门 考虑用 \(segment~tree~beats\) 那一套理论,维护区间最小值 \(mn\) 和严格次小值 \(se\) 那么可以直接 \(mlog^2n\) 维护前三个操作 考虑维护历史最 ...

  9. 【bzoj3894】文理分科 网路流

    [bzoj3894]文理分科 2015年3月25日3,4002 Description  文理分科是一件很纠结的事情!(虽然看到这个题目的人肯定都没有纠 结过)  小P所在的班级要进行文理分科.他的班 ...

随机推荐

  1. PMM--简介与部署

    一.PMM架构 Percona Monitoring and Management是percona一款开源的用于管理和监控MySQL 和MongoDB性能的开源平台,通过PMM客户端收集到的DB监控数 ...

  2. win10 查看本机的激活秘钥

    系统的注册表中,找到如下位置 计算机\HKEY_LOCAL_MACHINE\SOFTWARE\Microsoft\Windows NT\CurrentVersion\SoftwareProtectio ...

  3. 用JS实现输出两个数范围内的随机数

    const rs = require("readline-sync"); function roundNum(min = 0, max = 0) { if (!isNaN(min) ...

  4. Androidx初尝及其新旧包对照表

    x的最低实验条件 AndroidStudio 3.2.0+ gradle:gradle-4.6以上 本次实验条件: AndroidStudio 3.3 (强制要求最低gradle版本为gradle-4 ...

  5. Linux网卡上添加多个IP

    1.查看现有网卡ip addr show看到ens_xxx 2.进入网卡配置项cd /etc/sysconfig/network-scripts/vim ifcfg-ens192末尾添加2行:IPAD ...

  6. Java中重载(overloading)和重写(Overriding)的区别

    一:方法的重载 (1)方法重载指在类中定义方法名相同,参数不同的不同的多个方法(返回值类型可随意,不能以返回类型作为重载函数的区分标准). 参数不同表现: 1.参数的个数不同 2.参数的类型不同 3. ...

  7. 过滤器+用session验证是否登陆过

    过滤器: public class MyActionFilter : ActionFilterAttribute//继承ActionFilterAttribute类 { public override ...

  8. 2019.7月-前端面试总结(H5+C3+JS+ES6+Vue+浏览器)

    第二次面试 HTML HTML5中的新标签,举例一下 canvas绘画,本地离线存储localStorage,sessionStorage,video和audio元素,语义化元素,表单类型(date, ...

  9. Linux正则表达式,grep总结,sed用法

    原文: 1.sed   流编辑器,实现对文字的增删改替换查(过滤.取行),能同时处理多个文件多行的内容,可以不对原文件改动,把整个文件 输入到屏幕,可以把只匹配到模式的内容输入到屏幕上.还可以对原文件 ...

  10. C#-判断字符是否是全角半角

    C#字符串的全角是指用二个字节来表示的一个字符 C#字符串的半角是用一个字节来表示的一个字符 这样的话我们就可以用string.length 和System.text.Encoding.Default ...