「luogu - P4126」「ahoi 2009」最小割
也许题不错,反正有点降智…
先给结论,在
E_N=E \\
c(x,y)=w(x,y)
\]
的流网络中:
- 可行边:在增广完的 induced subgraph 中,不存在 \(u\) 到 \(v\) 的路径;
- 必要边:在增广完的 induced subgraph 中,可以从 \(S\) 到 \(u\) 且可以从 \(v\) 到 \(T\)。
先看可行边。不存在 \((u,v)\) 有两个条件,一个是 \(c_f(u,v)=0\),另一个是与之并联(特指以 \(u\) 为起点,\(v\) 为终点的)的线路中存在 \(c_f(u',v')=0\)。第一个的理解是,如果它没满流,则与之串联的 arcs 中存在比它的容量更小的边,根据最小割串联割最小的原则成立;第二个就是,如果你不把并联的砍了,你的划分压根不合法,何谈可行与否。
那么关于可行边的判断,把图缩点即可。
再看必要边。这个类比可行边即可,不赘述。
#include<bits/stdc++.h>
using namespace std;
int n,m,S,T,co[4100],dfsnt,colnt,inst[4100],sta[4100],top,dfn[4100],low[4100],vi[4100],rec[60100];
vector<pair<int,int>> arc;
template<typename T> struct network {
const int n;
struct edge {
int to,r; T w;
// friend bool operator<(const edge& one,const edge& ano) { return one.to<ano.to || (one.to==ano.to && one.r<ano.r); }
};
vector<vector<edge>> e;
vector<int> lev,iter;
network(int n):n(n),e(n+1),lev(n+1),iter(n+1) {}
void link(const int x,const int y,T w,int ID) {
assert(1<=x && x<=n && 1<=y && y<=n);
rec[ID]=int(e[x].size());
e[x].push_back((edge){y,int(e[y].size())+(x==y),w});
e[y].push_back((edge){x,int(e[x].size())-1,0});
}
bool BFS(const int s,const int t) {
queue<int> q; lev.assign(n+1,0);
for(q.push(s),lev[s]=1; q.size(); q.pop()) {
for(int now=q.front(),i=iter[now]=0,y; i<int(e[now].size()); ++i) {
if(now==t) return 1;
if(!lev[y=e[now][i].to] && e[now][i].w) lev[y]=lev[now]+1,q.push(y);
}
}
return lev[t];
}
T DFS(const int now,T f,const int t) {
if(now==t) return f;
T res=0,tt;
for(int& i=iter[now],y; i<int(e[now].size()); ++i) {
if(lev[y=e[now][i].to]==lev[now]+1 && e[now][i].w && (tt=DFS(y,min(f,e[now][i].w),t))) {
e[now][i].w-=tt; e[y][e[now][i].r].w+=tt; res+=tt; f-=tt;
if(!f) break;
}
}
if(!res) lev[now]=0;
return res;
}
T get(const int s,const int t) {
T res=0;
while(BFS(s,t)) res+=DFS(s,numeric_limits<T>::max(),t);
return res;
}
};
template<typename T=int> T rd() {
T x=0; char IO=getchar(); bool f=0;
while(IO<'0' || IO>'9') f|=IO=='-',IO=getchar();
while(IO>='0' && IO<='9') x=x*10+(IO&15),IO=getchar();
return f?-x:x;
}
void DFS(const int now,const network<int>& G) {
dfn[now]=low[now]=++dfsnt;
inst[sta[++top]=now]=1;
for(int i=0,y; i<int(G.e[now].size()); ++i) {
if(!G.e[now][i].w) continue;
if(!dfn[y=G.e[now][i].to]) DFS(y,G),low[now]=min(low[now],low[y]);
else if(inst[y]) low[now]=min(low[now],dfn[y]);
}
if(dfn[now]==low[now]) {
++colnt;
while(sta[top]!=now) co[sta[top]]=colnt,inst[sta[top]]=0,top--;
top--; co[now]=colnt; inst[now]=0;
}
}
void DFS_network(const int now,const network<int>& G,const int c) {
vi[now]=c;
for(int i=0,y; i<int(G.e[now].size()); ++i) {
if(!vi[y=G.e[now][i].to] && (c-1?G.e[y][G.e[now][i].r].w:G.e[now][i].w)) DFS_network(y,G,c);
}
}
signed main() {
freopen("mincut.in","r",stdin);
freopen("mincut.out","w",stdout);
n=rd(); m=rd(); S=rd(); T=rd();
network<int> G(n);
for(int i=0,x,y; i<m; ++i) x=rd(),y=rd(),G.link(x,y,rd(),i),arc.emplace_back(x,y);
for(int i=(G.get(S,T),1); i<=n; ++i) if(!dfn[i]) DFS(i,G);
DFS_network(S,G,1); DFS_network(T,G,2);
// for(int i=1; i<=n; ++i) printf(" --- %d ",vi[i]); puts("\n");
// for(int i=1; i<=n; ++i) printf(" --- %d ",co[i]); puts("");
// for(int i=1; i<=n; ++i) sort(G.e[i].begin(),G.e[i].end());
// for(int now=1; now<=n; ++now) {
// printf(" --- current node = %d:",now);
// for(int i=0; i<int(G.e[now].size()); ++i) printf(" %d",G.e[now][i].to);
// puts("");
// }
for(int i=0; i<m; ++i) {
// printf(" (%d %d)\n",arc[i].first,lower_bound(G.e[arc[i].first].begin(),G.e[arc[i].first].end(),(network<int>::edge){arc[i].second,0,0})->to);
// if(lower_bound(G.e[arc[i].first].begin(),G.e[arc[i].first].end(),(network<int>::edge){arc[i].second,0,0})->w) puts("0 0");
// printf(" (%d %d)[%d %d]\n",arc[i].first,G.e[arc[i].first][rec[i]].to,co[arc[i].first],co[G.e[arc[i].first][rec[i]].to]);
if(G.e[arc[i].first][rec[i]].w) puts("0 0");
else printf("%d %d\n",co[arc[i].first]!=co[arc[i].second],vi[arc[i].first]==1 && vi[arc[i].second]==2);
}
return 0;
}
「luogu - P4126」「ahoi 2009」最小割的更多相关文章
- 「Luogu P2060 [HNOI2006]马步距离」
一道神奇的BFS 前置芝士 BFS(DFS):这次真的不是我懒,我也不知道DFS怎么写. STL中的set或者map. 具体做法 数据范围非常大,直接BFS肯定是一片黑色(指TLE,MLE),直接贪心 ...
- 「SHOI2007」「Codevs2341」 善意的投票(最小割
2341 善意的投票 2007年省队选拔赛上海市队选拔赛 时间限制: 5 s 空间限制: 128000 KB 题目等级 : 大师 Master 题目描述 Description 幼儿园里有n个小朋 ...
- BZOJ1001/LG4001 「ICPC Beijing2006」狼抓兔子 平面图最小割转对偶图最短路
问题描述 BZOJ1001 LG4001 题解 平面图最小割=对偶图最短路 假设起点和终点间有和其他边都不相交的一条虚边. 如图,平面图的若干条边将一个平面划分为若干个图形,每个图形就是对偶图中的一个 ...
- 【Luogu】P2057善意的投票(最小割转最大流)
题目链接 也算水题一道吧,不过Round1感性理解一下就xjb建了个图,40 Round2仔细分析了一会,理性建了个图,90 然后分析了半天……改大数组就A了…… 从S到所有值为1的点连一条inf的边 ...
- 【Luogu】P3355骑士共存问题(最小割)
题目链接 像题面那样把棋盘染成红黄点.发现骑士迈一步能到达的点的颜色一定是跟他所在的格子的颜色不同的.于是(woc哪来的于是?这个性质有这么明显吗?)从源点向所有红点连边,从所有黄点向汇点连边,红点向 ...
- BZOJ 4823 Luogu P3756 老C的方块 染色+最小割
题面太长了请各位自行品尝—>老C的方块 分析: 我们要解决掉所有使人弃疗的组合,还要保证花费最小,容易想到最小割(当然你要是想费用流的话,我们就没办法定义流量了) 我们来分析一下那些令人弃疗的组 ...
- Luogu SP839 OPTM - Optimal Marks(按位最小割)
这道题和 BZOJ 2400 是一道题,不多讲了 CODE #include <cstdio> #include <cstring> #include <vector&g ...
- loj#2020 「AHOI / HNOI2017」礼物 ntt
loj#2020 「AHOI / HNOI2017」礼物 链接 bzoj没\(letex\),差评 loj luogu 思路 最小化\(\sum\limits_1^n(a_i-b_i)^2\) 设改变 ...
- Loj #2495. 「AHOI / HNOI2018」转盘
Loj #2495. 「AHOI / HNOI2018」转盘 题目描述 一次小 G 和小 H 原本准备去聚餐,但由于太麻烦了于是题面简化如下: 一个转盘上有摆成一圈的 \(n\) 个物品(编号 \(1 ...
- Loj #2494. 「AHOI / HNOI2018」寻宝游戏
Loj #2494. 「AHOI / HNOI2018」寻宝游戏 题目描述 某大学每年都会有一次 Mystery Hunt 的活动,玩家需要根据设置的线索解谜,找到宝藏的位置,前一年获胜的队伍可以获得 ...
随机推荐
- Galaxy Project | 一些尝试与思考
很久都没有更新推文了,脑壳羞涩,快码不出字的节奏! 最近在尝试内部 Galaxy 一些新工具的开发和 Galaxy 核心版本的升级测试,发现一些问题,简单记录和聊一下吧. 一些尝试 对于在线的 web ...
- CKS 考试题整理 (11)-沙箱运行容器gVisor
Context 该 cluster使用 containerd作为CRI运行时.containerd的默认运行时处理程序是runc. containerd已准备好支持额外的运行时处理程序runsc (g ...
- 【技术积累】C语言中基础知识【三】
什么是C语言[了解即可] C语言是一种通用的高级编程语言,由美国贝尔实验室的Dennis Ritchie在20世纪70年代早期开发出来的.它在计算机科学和软件工程领域中被广泛使用. C语言具有以下特点 ...
- php处理emoji表情 存数据库
PHP 处理emoji表情 存数据库 直接过滤掉 1 function filter_emoji($str) { 2 $regex = '/(\\\u[ed][0-9a-f]{3})/i'; 3 $s ...
- Unity中的RegisterPlugins:深入解析与实用案例
Unity中的RegisterPlugins:深入解析与实用案例 在Unity游戏开发中,我们经常需要使用第三方插件来实现一些特定的功能.为了让这些插件能够在Unity中正常工作,我们需要对它们进行注 ...
- 【从0开始编写webserver·基础篇#03】TinyWeb源码阅读,还是得看看靠谱的项目
[前言] 之前通过看书.看视频和博客拼凑了一个webserver,然后有一段时间没有继续整这个项目 现在在去看之前的代码,真的是相当之简陋,而且代码设计得很混乱,我认为没有必要继续在屎堆上修改了,于是 ...
- python学习笔记:第七章面向对象
与java类似,python作为一种面向对象的编程语言,也可以创建自定义的对象和类. 它的特性主要有:继承,封装,多态,方法,属性,超类 1.变量的作用域 c = 50 #全局变量, 作用域为整个模块 ...
- SSE图像算法优化系列三十一:RGB2HSL/RGB2HSV及HSL2RGB/HSV2RGB的指令集优化-上。
RGB和HSL/HSV颜色空间的相互转换在我们的图像处理中是有着非常广泛的应用的,无论是是图像调节,还是做一些肤色算法,HSL/HSV颜色空间都非常有用,他提供了RGB颜色空间不具有的一些独特的特性, ...
- 2023年ccpc大学生程序设计竞赛-crf
第一次参加这种大型线下程序设计比赛,心情自然是很激动,但比赛中并没有想象中那么顺利,甚至可以说有些惊险,比赛开始的时候,我们三人随便看看题,顺便等着跟榜,对于签到题我们自然是相信clk可以很快地独立完 ...
- 一对多数据的多选筛选([EF,ABP]但是sql思路通用)
一对多数据的多选筛选([EF,ABP]但是sql思路通用) 标题起的并不是特别准确,但是现在又没有更加准确的描述,暂时这样了 业务需求: 数据格式为一条数据,对应多条二级数据,暂时表示为主信息为&qu ...