[ZJOI2009] 狼与羊的故事 - 最小割
给定一个\(N \times M\)方格矩阵,每个格子可在\(0,1,2\)中取值。要求在方格的边上进行划分,使得任意联通块内不同时包含\(1\)和\(2\)的格子。
若对方格矩阵的形式感到困惑,我们可以先考虑图上的形式。
假设图中有节点\(W_1\),\(Q_1\),\(S_1\),\(S_2\),且\(W_1\)与其它三者间各有一边,此外\(Q_1\)和\(S_1\)间还有一边。\(W_1\)是1值节点,\(S_1,S_2\)是2值节点。现在考虑进行划分。
我们发现在原问题中,\(0\)值点无论和哪个点划在一个连通块内都是合法的。换言之,\(0\)值节点在这里只起到传导连通的作用。亦即任意两个连通节点间可以有任意个\(0\)值节点。
根据这一性质,我们对图进行分层。建立网络,将\(W_i\)作为第一层(靠近源点),\(S_i\)作为最后一层(靠近汇点),\(Q_i\)放中间。将原图中的连接关系单向地对应到网络中:
- 由\(W_i\)向任意点的连接均转化为一条容量为\(1\)的弧
- 由\(S_i\)向任意点的连接均被忽略
- 由\(Q_i\)向非\(W_i\)点的连接均转化为一条容量为\(1\)的弧
容易发现,由于原图是无向图,我们在建立流网络时对其定向,从而避免对一个划分边的重复计算。而定向的依据则是依据\(W_i > Q_i > S_i\)的优先级顺序进行的。
这样建图后,求出的最大流就是流网络的最小割,也就是答案。
分析完图上的情况后,这里的情形就显得简单了。
我们考虑对任意格子\((i,j)\),其四周的格子即是上文中与之直接相连的节点。我们将图中的关系放到方格矩阵中处理就可以了。
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<vector>
#include<queue>
#define G g
using namespace std;
struct Item{int p,v,c;};
Item item(int _p,int _v,int _c){Item it;it.p=_p;it.v=_v;it.c=_c;return it;}
vector <Item> g[100005];
int n,m,t1,t2,t3,dis[100005],d[100005],vis[100005],s,t,t4,costs=0,tans,ans,inc[100005],cnt=0,k,inp[100005][4],edg[100005];
vector <int> oppo[100005];
int dinic_spfa(){
memset(dis,0xff,sizeof dis);
memset(d,0x3f,sizeof d);
memset(vis,0x00,sizeof vis);
memset(inc,0x00,sizeof inc);
queue <int> q;
q.push(s); dis[s]=0; d[s]=0;
while(!q.empty()){
int p=q.front(); q.pop(); vis[p]=0; inc[p]++;
for(int i=0;i<g[p].size();i++){
if(d[g[p][i].p]>d[p]+g[p][i].v&&g[p][i].c>0){
dis[g[p][i].p]=p; edg[g[p][i].p]=i;
d[g[p][i].p]=d[p]+g[p][i].v;
if(vis[g[p][i].p]==0) vis[g[p][i].p]=1, q.push(g[p][i].p);}}}
return dis[t]>0;}
int dinic_dfs(){
int p=t,a=0x7fffffff;
while(p-s) a=min(a,g[dis[p]][edg[p]].c),p=dis[p];
int lc=costs;
p=t;
while(p-s){
int tc,i=edg[p];
G[dis[p]][i].c-=a,
G[p][oppo[dis[p]][i]].c+=a,
costs+=a*G[dis[p]][i].v;
p=dis[p];}
return a;}
int dinic_main(int src,int dest){
s=src; t=dest;
while(dinic_spfa()) ans+=dinic_dfs();
return ans;}
void build(int w,int x,int y,int z){
oppo[x].push_back(g[w].size());
g[w].push_back(item(x,z,y));
oppo[w].push_back(g[x].size());
g[x].push_back(item(w,-z,0));}
int main(){
scanf("%d%d%d",&n,&m,&k);
for(int i=1;i<=m;i++) scanf("%d%d%d%d",&inp[i][0],&inp[i][1],&inp[i][2],&inp[i][3]),
build(inp[i][0],inp[i][1],inp[i][2],0);
dinic_main(1,n);
printf("%d ",ans);
for(int i=1;i<=m;i++)
build(inp[i][0],inp[i][1],0x7fffffff,inp[i][3]);
build(0,1,k,0);
dinic_main(0,n);
printf("%d\n",costs);
}
[ZJOI2009] 狼与羊的故事 - 最小割的更多相关文章
- BZOJ 1412: [ZJOI2009]狼和羊的故事( 最小割 )
显然是最小割...把狼的领地连S, 羊的领地连T, 然后中间再连边, 跑最大流就OK了 -------------------------------------------------------- ...
- 【BZOJ1412】[ZJOI2009]狼和羊的故事 最小割
[BZOJ1412][ZJOI2009]狼和羊的故事 Description “狼爱上羊啊爱的疯狂,谁让他们真爱了一场:狼爱上羊啊并不荒唐,他们说有爱就有方向......” Orez听到这首歌,心想: ...
- P2598 [ZJOI2009]狼和羊的故事(最小割)
P2598 [ZJOI2009]狼和羊的故事 说真的,要多练练网络流的题了,这么简单的网络流就看不出来... 题目要求我们要求将狼和羊分开,也就是最小割,(等等什么逻辑...头大....) 我们这样想 ...
- BZOJ1412[ZJOI2009]狼和羊的故事——最小割
题目描述 “狼爱上羊啊爱的疯狂,谁让他们真爱了一场:狼爱上羊啊并不荒唐,他们说有爱就有方向......” Orez听到这首歌,心想:狼和羊如此和谐,为什么不尝试羊狼合养呢?说干就干! Orez的羊狼圈 ...
- LG2598/BZOJ1412 「ZJOI2009」狼和羊的故事 最小割
问题描述 LG2598 BZOJ1412 题解 看到要把狼和羊两个物种分开 自然想到最小割. 发现\((x,y)\)可以向上下左右走以获得贡献,所以建边:\((x,y),(x-1,y)\),\((x, ...
- BZOJ1412 ZJOI2009 狼和羊的故事 【网络流-最小割】
BZOJ1412 ZJOI2009 狼和羊的故事 Description “狼爱上羊啊爱的疯狂,谁让他们真爱了一场:狼爱上羊啊并不荒唐,他们说有爱就有方向......” Orez听到这首歌,心想:狼和 ...
- BZOJ1412 [ZJOI2009]狼和羊的故事 【最小割】
1412: [ZJOI2009]狼和羊的故事 Time Limit: 10 Sec Memory Limit: 162 MB Submit: 3454 Solved: 1733 [Submit][ ...
- P2598 [ZJOI2009]狼和羊的故事(最小割)
P2598 [ZJOI2009]狼和羊的故事 题目描述 “狼爱上羊啊爱的疯狂,谁让他们真爱了一场:狼爱上羊啊并不荒唐,他们说有爱就有方向......” Orez听到这首歌,心想:狼和羊如此和谐,为什么 ...
- P2598 [ZJOI2009]狼和羊的故事(网络流)
P2598 [ZJOI2009]狼和羊的故事 源点和所有狼连 $inf$ 的边 所有羊和汇点连 $inf$ 的边 所有点向四周连 $1$ 的边 这样所有狼和羊之间的边都被割掉了 统计最小割就好辣 #i ...
随机推荐
- JavaScript数据类型typeof()和转换
javascript属于弱类型,值包含:数字,字符串和布尔值,c++与java属于强类型;int a="a",string a="a" 报错;var a ;原始 ...
- openlayers6结合geoserver实现地图矢量瓦片(附源码下载)
内容概览 1.基于openlayers6结合geoserver实现地图矢量瓦片2.源代码demo下载 效果图如下: 实现思路:利用Geoserver发布矢量切片服务,然后openlayers调用矢量瓦 ...
- 后端跨域的N种方法
简单来说,CORS是一种访问机制,英文全称是Cross-Origin Resource Sharing,即我们常说的跨域资源共享,通过在服务器端设置响应头,把发起跨域的原始域名添加到Access-Co ...
- c#winform自定义窗体,重绘标题栏,自定义控件学习
c#winform自定义窗体,重绘标题栏 虽然现在都在说winform窗体太丑了,但是我也能尽量让桌面应用程序漂亮那么一点点话不多说,先上图 重绘标题栏先将原生窗体设置成无边框,FormBoderSt ...
- 如何在Mac和Windows PC之间无线共享文件
有时候,我需要在Mac和PC之间无线共享文件.由于并非所有人都在使用macOS,因此无论是在办公室还是在家里,这种情况都会发生.尽管并非一帆风顺,但有一种无需任何第三方应用程序即可弥合差距的方法. 根 ...
- python——面向对象(1),基础
"""面向对象:抽象化编程思想.类,对象:用类来创建(实例化)对象.类:一系列特征和行为相同的事物总和, 1.属性:特征 2.行为:方法 定义类 :PEP 8要求标识符的 ...
- C# WPF简况(2/3)
微信公众号:Dotnet9,网站:Dotnet9,问题或建议:请网站留言, 如果对您有所帮助:欢迎赞赏. C# WPF简况(2/3) 阅读导航 本文背景 代码实现 本文参考 1.本文背景 承接上文(C ...
- amr格式转mp3和直接播放amr格式的文件-sunziren
原创文章,转载请注明出处! 前言: amr作为一种高压缩比的音频格式,受到很多客户的青睐.本文主要涉及两部分的内容,一是amr如何转为mp3格式,二是如何直接播放amr格式的文件. 1. 如何使用Ja ...
- PTA 邻接表存储图的广度优先遍历
试实现邻接表存储图的广度优先遍历. 函数接口定义: void BFS ( LGraph Graph, Vertex S, void (*Visit)(Vertex) ) 其中LGraph是邻接表存储的 ...
- 关于iScroll在安卓移动端/chrome模拟移动端上下滑动卡顿问题处理!!!!真实可靠!!!已解决!!!
滑动卡顿效果 安卓手机打开微信浏览网页,Chrome模拟手机浏览网页,都出现的问题滑动卡顿! 修改代码点: 1. <style type="text/css"> ...