Luogu4897 【模板】最小割树
没事干写一发模板。
- #include<iostream>
- #include<cstdio>
- #include<cmath>
- #include<cstdlib>
- #include<cstring>
- #include<algorithm>
- #include<map>
- using namespace std;
- int read()
- {
- int x=,f=;char c=getchar();
- while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
- while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
- return x*f;
- }
- #define N 510
- #define M 1600
- #define inf 1000000000
- int n,m,p[N],v[N],u[N],tot=,t=-;
- bool flag[N];
- int d[N],cur[N],q[N];
- struct data{int to,nxt,cap,flow;
- }edge[M<<];
- map<int,bool> f;
- namespace tree
- {
- int p[N],fa[N],deep[N],len[N],t=;
- struct data{int to,nxt,len;}edge[M<<];
- void addedge(int x,int y,int z){t++;edge[t].to=y,edge[t].nxt=p[x],edge[t].len=z,p[x]=t;}
- void dfs(int k)
- {
- for (int i=p[k];i;i=edge[i].nxt)
- if (edge[i].to!=fa[k])
- {
- deep[edge[i].to]=deep[k]+;
- fa[edge[i].to]=k;
- len[edge[i].to]=edge[i].len;
- dfs(edge[i].to);
- }
- }
- void solve(int q)
- {
- while (q--)
- {
- int x=read(),y=read(),ans=;
- if (deep[x]<deep[y]) swap(x,y);
- while (deep[x]>deep[y]) ans=min(ans,len[x]),x=fa[x];
- while (x!=y) ans=min(ans,len[x]),ans=min(ans,len[y]),x=fa[x],y=fa[y];
- printf("%d\n",ans);
- }
- }
- }
- void addedge(int x,int y,int z)
- {
- t++;edge[t].to=y,edge[t].nxt=p[x],edge[t].cap=z,edge[t].flow=,p[x]=t;
- }
- bool bfs(int S,int T)
- {
- memset(d,,sizeof(d));d[S]=;
- int head=,tail=;q[]=S;
- do
- {
- int x=q[++head];
- for (int i=p[x];~i;i=edge[i].nxt)
- if (d[edge[i].to]==-&&edge[i].flow<edge[i].cap)
- {
- d[edge[i].to]=d[x]+;
- q[++tail]=edge[i].to;
- }
- }while (head<tail);
- return ~d[T];
- }
- int work(int k,int T,int f)
- {
- if (k==T) return f;
- int used=;
- for (int i=cur[k];~i;i=edge[i].nxt)
- if (d[k]+==d[edge[i].to])
- {
- int w=work(edge[i].to,T,min(f-used,edge[i].cap-edge[i].flow));
- edge[i].flow+=w,edge[i^].flow-=w;
- if (edge[i].flow<edge[i].cap) cur[k]=i;
- used+=w;if (used==f) return f;
- }
- if (used==) d[k]=-;
- return used;
- }
- void dinic(int S,int T)
- {
- for (int i=;i<=t;i++) edge[i].flow=;
- int ans=;
- while (bfs(S,T))
- {
- memcpy(cur,p,sizeof(p));
- ans+=work(S,T,inf);
- }
- tree::addedge(S,T,ans);
- }
- void dfs(int k)
- {
- flag[k]=;
- for (int i=p[k];~i;i=edge[i].nxt)
- if (!flag[edge[i].to]&&edge[i].flow<edge[i].cap)
- dfs(edge[i].to);
- }
- void solve(int l,int r)
- {
- if (l>=r) return;
- dinic(v[l],v[r]);
- memset(flag,,sizeof(flag));
- dfs(v[l]);
- int cnt=l-;
- for (int i=l;i<=r;i++)
- if (flag[v[i]]) u[++cnt]=v[i];
- cnt=r+;
- for (int i=l;i<=r;i++)
- if (!flag[v[i]]) u[--cnt]=v[i];
- for (int i=l;i<=r;i++) v[i]=u[i];
- solve(l,cnt-);
- solve(cnt,r);
- }
- int main()
- {
- #ifndef ONLINE_JUDGE
- freopen("mincuttree.in","r",stdin);
- freopen("mincuttree.out","w",stdout);
- const char LL[]="%I64d\n";
- #else
- const char LL[]="%lld\n";
- #endif
- n=read(),m=read();
- memset(p,,sizeof(p));
- for (int i=;i<=m;i++)
- {
- int x=read(),y=read(),z=read();
- addedge(x,y,z),addedge(y,x,z);
- }
- for (int i=;i<=n;i++) v[i]=i;
- solve(,n);
- tree::dfs();
- int q=read();
- tree::solve(q);
- return ;
- }
Luogu4897 【模板】最小割树的更多相关文章
- [模板]最小割树(Gomory-Hu Tree)(luogu4897)
给定一个\(n\)个点\(m\)条边的无向连通图,多次询问两点之间的最小割 两点间的最小割是这样定义的:原图的每条边有一个割断它的代价,你需要用最小的代价使得这两个点不连通 Input 第一行两个数\ ...
- 洛谷.4897.[模板]最小割树(Dinic)
题目链接 最小割树模板.具体见:https://www.cnblogs.com/SovietPower/p/9734013.html. ISAP不知为啥T成0分了.. Dinic: //1566ms ...
- CQOI2016 不同的最小割 (最小割树模板)(等价流树的Gusfield构造算法)
题目 最小割树模板 算法详解及证明见: 2016年国家队候选队员论文 <浅谈无向图最小割问题的一些算法及应用--绍兴一中 王文涛> 3.2节 CODE #include <bits/ ...
- 【模板】最小割树(Gomory-Hu Tree)
传送门 Description 给定一个\(n\)个点\(m\)条边的无向连通图,多次询问两点之间的最小割 两点间的最小割是这样定义的:原图的每条边有一个割断它的代价,你需要用最小的代价使得这两个点不 ...
- [洛谷]P3729 曼哈顿计划EX(最小割树/等价流树)
题目大意:给出一张n个点m条边的无向图,每个点有点权,q次询问,每次给出k,要求选出若干个点点权之和不小于k,求一个最大的值x,使得选出的点中任意两点之间至少有x条互不相交的链.(n<=550, ...
- BZOJ4519[Cqoi2016]不同的最小割——最小割树+map
题目描述 学过图论的同学都知道最小割的概念:对于一个图,某个对图中结点的划分将图中所有结点分成 两个部分,如果结点s,t不在同一个部分中,则称这个划分是关于s,t的割.对于带权图来说,将 所有顶点处在 ...
- BZOJ2229[Zjoi2011]最小割——最小割树
题目描述 小白在图论课上学到了一个新的概念——最小割,下课后小白在笔记本上写下了如下这段话: “对于一个图,某个对图中结点的划分将图中所有结点分成两个部分,如果结点s,t不在同一个部分中,则称这个划分 ...
- 最小割树(Gomory-Hu Tree)求无向图最小割详解 附 BZOJ2229,BZOJ4519题解
最小割树(Gomory-Hu Tree) 前置知识 Gomory-Hu Tree是用来解决无向图最小割的问题的,所以我们需要了解无向图最小割的定义 和有向图类似,无向图上两点(x,y)的割定义为一个边 ...
- [学习笔记]最小割树(Gomory-Hu Tree)
最小割树(\(\mathcal{Gomory-Hu Tree}\))简明指南 对于单源最短路径,我们有\(SPFA\)和\(Dijkstra\),对于多源最短路径,我们有\(Floyd\):对于两点间 ...
- 最小割树(Gomory-Hu Tree)
当我们遇到这样的问题: 给定一个 \(n\) 个点 \(m\) 条边的无向连通图,多次询问两点之间的最小割 我们通常要用到最小割树. 博客 建树 分治.记录当前点集,然后随便找俩点当 \(s\) 和 ...
随机推荐
- BZOJ1010_玩具装箱toy_KEY
题目传送门 这道题可以很快想到暴力DP的做法: f[i]=min(f[i],f[j]+(C[i]-C[j]+i-j--L)^); 但是数据范围有50000,这就需要用斜率优化了. 我们设S[i]=C[ ...
- 4361: isn
4361: isn https://lydsy.com/JudgeOnline/problem.php?id=4361 分析: dp+容斥. 首先计算出每个长度有多少种子序列是非降的.这一步可以$n^ ...
- 海思NB-IOT的SDK里面写KV键值对值掉电保存参数
1. 写函数,ICCID_CTRL_KEY 这个是键,有些是海思已经使用的,所以自己定义的时候要错开,函数的第二个参数是数据指针,第三个参数是数据长度 #define ICCID_CTRL_KEY ...
- 强化学习读书笔记 - 11 - off-policy的近似方法
强化学习读书笔记 - 11 - off-policy的近似方法 学习笔记: Reinforcement Learning: An Introduction, Richard S. Sutton and ...
- itchat个人练习 语音与文本图灵测试例程
背景介绍 itchat是一个开源的微信个人号接口,使用python调用微信从未如此简单. 使用不到三十行的代码,你就可以完成一个能够处理所有信息的微信机器人. 官方文档参考https://itchat ...
- uvaoj 156Ananagrams(map和vector组合使用)
https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem& ...
- this指向问题(箭头函数没有this)
全局环境下,指向windows console.log(this.document === document); // true 函数上下文调用 function f1(){ return this; ...
- div不设置高度背景颜色或外边框不能显示的解决方法
在使用div+css进行网页布局时,如果外部div有背景颜色或者边框,而不设置其高度,在浏览时出现最外层Div的背景颜色和边框不起作用的问题. 大体结构<div class="oute ...
- 如何在 Debian 9 下安装 LEMP 和 WHMCS 7.5
WHMCS 7.5 发布了,它开始支持 PHP 7.2,这里就写个简单的教程记录一下安装方式. 1.准备工作 首先,我们需要按照 在Debian 9 / Debian 8 下使用源安装方式安装 LEM ...
- 75.[LeetCode] Sort Colors
Given an array with n objects colored red, white or blue, sort them in-place so that objects of the ...