2019 Multi-University Training Contest 1 E Path(最短路+最小割)
题意
链接:https://vjudge.net/problem/HDU-6582
给定一个有向图,可以有重边,每条边上有一个权值表示删掉这条边的代价,问最少花费多少代价能使从s到t节点的最短路径增大?1≤n,m≤10000
思路
容易想到应该是删最短路上的边,最短路可能不止一条,所以使原图1到n的所有最短路不连通即可,这就是最小割呀!选出权值和最小的边使得图不连通,这里是使最短路图不连通。
所以做法就是先建两个图,一个是u->v的有向边,另一个是v->u的有向边,从1跑一下Dijkstra,从n跑一下Dijkstra,我们枚举每条边(u,v),如果这条边是最短路上的边,那么应满足dis1[u]+val(u,v)+dis2[v]=dis1[n],即1到u的最短距离+u到v的距离+v到n的最短距离=1到n的最短距离。
最小割=最大流,跑一下Dinic即可啦。
代码
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define int ll
const int N=1e6+5;
const ll inf=1e16;
struct node
{
int p,w;
node(int a,int b)
{
p=a;
w=b;
}
friend bool operator<(node a,node b) //权值小的先出队
{
if(a.w!=b.w) return a.w>b.w;
return a.p>b.p;
}
};
vector <node> eg1[N],eg2[N];
int dis1[N],n,dis2[N];
void add1(int u,int v,int w)
{
eg1[u].push_back(node(v,w));
}
void add2(int u,int v,int w)
{
eg2[u].push_back(node(v,w));
}
void Dijkstra1(int now)
{
for(int i=0; i<=n; i++) dis1[i]=inf;
dis1[now]=0;
priority_queue <node> pq;
pq.push(node(now,dis1[now]));
while(!pq.empty())
{
node f=pq.top();
pq.pop();
for(int i=0; i<eg1[f.p].size(); i++)
{
node t=eg1[f.p][i];
if(dis1[t.p]>t.w+f.w)
{
dis1[t.p]=t.w+f.w;
pq.push(node(t.p,dis1[t.p]));
}
}
}
}
void Dijkstra2(int now)
{
for(int i=0; i<=n; i++) dis2[i]=inf;
dis2[now]=0;
priority_queue <node> pq;
pq.push(node(now,dis2[now]));
while(!pq.empty())
{
node f=pq.top();
pq.pop();
for(int i=0; i<eg2[f.p].size(); i++)
{
node t=eg2[f.p][i];
if(dis2[t.p]>t.w+f.w)
{
dis2[t.p]=t.w+f.w;
pq.push(node(t.p,dis2[t.p]));
}
}
}
}
/***************************************/
int x,y,z,maxflow,deep[N];//deep深度
struct Edge
{
int next,to,dis;
} edge[N];
int num_edge=-1,head[N],cur[N];//cur用于复制head
queue <int> q;
void add_edge(int from,int to,int dis,bool flag)
{
edge[++num_edge].next=head[from];
edge[num_edge].to=to;
if (flag) edge[num_edge].dis=dis;//反图的边权为 0
head[from]=num_edge;
}
//bfs用来分层
bool bfs(int s,int t)
{
memset(deep,0x7f,sizeof(deep));
while (!q.empty()) q.pop();
for (int i=1; i<=n; i++) cur[i]=head[i];
deep[s]=0;
q.push(s);
while (!q.empty())
{
int now=q.front();
q.pop();
for (int i=head[now]; i!=-1; i=edge[i].next)
{
if (deep[edge[i].to]>inf && edge[i].dis)//dis在此处用来做标记 是正图还是返图
{
deep[edge[i].to]=deep[now]+1;
q.push(edge[i].to);
}
}
}
if (deep[t]<inf) return true;
else return false;
}
//dfs找增加的流的量
int dfs(int now,int t,int limit)//limit为源点到这个点的路径上的最小边权
{
if (!limit || now==t) return limit;
int flow=0,f;
for (int i=cur[now]; i!=-1; i=edge[i].next)
{
cur[now]=i;
if (deep[edge[i].to]==deep[now]+1 && (f=dfs(edge[i].to,t,min(limit,edge[i].dis))))
{
flow+=f;
limit-=f;
edge[i].dis-=f;
edge[i^1].dis+=f;
if (!limit) break;
}
}
return flow;
}
void Dinic(int s,int t)
{
while (bfs(s,t))
maxflow+=dfs(s,t,inf);
}
/***********************************************/
signed main()
{
int u,v,m,w,t;
scanf("%lld",&t);
while(t--)
{
scanf("%lld%lld",&n,&m);
for(int i=0; i<=n; i++) eg1[i].clear(),eg2[i].clear();
for(int i=0; i<m; i++)
{
scanf("%lld%lld%lld",&u,&v,&w);
add1(u,v,w);
add2(v,u,w);
}
Dijkstra1(1);
Dijkstra2(n);
memset(head,-1,sizeof(head));
for(int i=1; i<=n; i++)
{
int sz=eg1[i].size();
for(int j=0; j<sz; j++)
{
if(dis1[i]+eg1[i][j].w+dis2[eg1[i][j].p]==dis1[n])
{
add_edge(i,eg1[i][j].p,eg1[i][j].w,1);
add_edge(eg1[i][j].p,i,eg1[i][j].w,0);
}
}
}
maxflow=0;
Dinic(1,n);
printf("%lld\n",maxflow);
}
return 0;
}
2019 Multi-University Training Contest 1 E Path(最短路+最小割)的更多相关文章
- 2019 Multi-University Training Contest 1 Path(最短路+最小割)
题意:给你n个点 m条边 现在你能够堵住一些路 问怎样能让花费最少且让1~n走的路比最短路的长度要长 思路:先跑一边最短路 建一个最短路图 然后我们跑一边最大流求一下最小割即可 #include &l ...
- [2019杭电多校第一场][hdu6582]Path(最短路&&最小割)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6582 题意:删掉边使得1到n的最短路改变,删掉边的代价为该边的边权.求最小代价. 比赛时一片浆糊,赛后 ...
- 2019HDU多校Path——最短路最小割
题目 给出一个 $n$ 个顶点 $m$ 条边的图,要求阻塞一些边,使得从 $1$ 到 $n$ 的最短路变长,求阻塞的边长度和的最小值,不必保证阻塞后可达. 分析 很显然,要阻塞的边肯定在最短路图上,先 ...
- HDU - 6582 Path (最短路+最小割)
题意:给定一个n个点m条边的有向图,每条边有个长度,可以花费等同于其长度的代价将其破坏掉,求最小的花费使得从1到n的最短路变长. 解法:先用dijkstra求出以1为源点的最短路,并建立最短路图(只保 ...
- 2019 Nowcoder Multi-University Training Contest 4 E Explorer
线段树分治. 把size看成时间,相当于时间 $l$ 加入这条边,时间 $r+1$ 删除这条边. 注意把左右端点的关系. #include <bits/stdc++.h> ; int X[ ...
- 2019 Nowcoder Multi-University Training Contest 1 H-XOR
由于每个元素贡献是线性的,那么等价于求每个元素出现在多少个异或和为$0$的子集内.因为是任意元素可以去异或,那么自然想到线性基.先对整个集合A求一遍线性基,设为$R$,假设$R$中元素个数为$r$,那 ...
- [2019杭电多校第二场][hdu6598]Harmonious Army(最小割)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6598 题意是说一个军队有n人,你可以给他们每个人安排战士或者法师的职业,有m对人有组合技,组合技的信息 ...
- [最短路,最大流最小割定理] 2019 Multi-University Training Contest 1 Path
题目:http://acm.hdu.edu.cn/showproblem.php?pid=6582 Path Time Limit: 2000/1000 MS (Java/Others) Mem ...
- 2019 Multi-University Training Contest 1
2019 Multi-University Training Contest 1 A. Blank upsolved by F0_0H 题意 给序列染色,使得 \([l_i,r_i]\) 区间内恰出现 ...
随机推荐
- android 在基类activity中注册BroadcastReceiver,子activity类实现响应
android app 一般都会定义自己的BaseActivity, 如果各子Activity都需要接收广播但对广播的处理又不同时,可以考虑在BaseActivity中注册BroadcastRecei ...
- USB HOST与 USB OTG的区别及工作原理
在SmartQ 7上面,同时存在USB HOST与 USB OTG两个接口,我想问一下,这两个接口有什么区别么?我怎么认为HOST属于是多余呢? 麻烦高手解答,感激不尽!!! 转自USB HOST与 ...
- MATLAB实例:对称双随机矩阵
MATLAB实例:对称双随机矩阵 作者:凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ 双随机矩阵(doubly stochastic matrix):元素属 ...
- java8-06-四大函数式接口
在使用lamdba表达式需要函数式接口的支持 java8已经提供了很多函数式接口 在java.util.function包下
- 201871010112-梁丽珍《面向对象程序设计(java)》第十三周学习总结
项目 内容 这个作业属于哪个课程 https://www.cnblogs.com/nwnu-daizh/ 这个作业的要求在哪里 https://www.cnblogs.com/nwnu-daizh/p ...
- DRF--认证和权限
前戏 大家都知道http协议是无状态的,每次发送请求他们怎么知道我们是不是登录过呢?我们可以在用户登录之后给用户一个“暗号”,下次请求的时候带着这个“暗号”来.我们拿自己存的和携带过来的进行对比,如果 ...
- 无限可能,Elasticsearch(一)
这个世界已然被数据淹没.多年来,我们系统间流转和产生的大量数据已让我们不知所措. 现有的技术都集中在如何解决数据仓库存储以及如何结构化这些数据. 这些看上去都挺美好,直到你实际需要基于这些数据实时做决 ...
- golang数据结构之栈
stack.go package stack import ( "errors" "fmt" ) type Stack struct { MaxTop int ...
- linux--解决anaconda升级pip问题
Anaconda指的是一个开源的Python发行版本,其包含了conda.Python等180多个科学包及其依赖项. 在用pip install命令安装东西时,有时会提示如下错误:升级pip You ...
- 关于unittest单元测试框架中常用的几种用例加载方法
unittest模块是Python自带的一个单元测试模块,我们可以用来做单元测试.unittest模块包含了如下几个子模块: 测试用例:TestCase 测试集:TestSuite 加载用例:Test ...