codeforces 567 E. President and Roads 【 最短路 桥 】
给出一个有向图,从起点走到终点(必须走最短路),问一条边是否一定会被经过,如果不经过它,可以减小它的多少边权使得经过它(边权不能减少到0)
正反向建图,分别求出起点到每个点的最短距离,终点到每个点的最短距离(用这个可以算出减小的边权)
再将在最短路径上的边重新建图。求出里面的桥,就是必须经过的边
wa了一上午------呜呜呜呜
先wa 19 是因为求桥的时候是无向图,数组开小了一半
然后 wa 46 ,是因为dis[]数组初始化为 1 << 30 -1 ,应该再开大点 ,开成 1 << 50 -1
我写成 1 LL * 50 -----一直wa 19------
5555555555555555555--------sad-------
- #include <cstdio>
- #include <cstring>
- #include <cstdlib>
- #include <cmath>
- #include <vector>
- #include <map>
- #include <set>
- #include <stack>
- #include <queue>
- #include <iostream>
- #include <algorithm>
- using namespace std;
- #define lp (p << 1)
- #define rp (p << 1|1)
- #define getmid(l,r) (l + (r - l) / 2)
- #define MP(a,b) make_pair(a,b)
- typedef long long ll;
- typedef unsigned long long ull;
- typedef pair<int,int> pii;
- const int INF = ( << ) - ;
- const int maxn = ;
- int N,M;
- ll path;
- int first1[maxn],nxt1[ * maxn],ecnt1;
- int first2[maxn],nxt2[ *maxn],ecnt2;
- int first[maxn],nxt[*maxn],ecnt;
- ll dis1[maxn],dis2[maxn];
- int bg[*maxn],vis[*maxn],low[maxn],dfn[maxn];
- int ans[*maxn],res[maxn];
- int tot,num;
- struct edge{
- int v,u,cost;
- int tag;
- int ge;
- int id;
- friend bool operator < (edge a,edge b){
- return a.cost > b.cost;
- }
- };
- edge e1[*maxn],e2[*maxn],e[*maxn];
- void init(){
- ecnt1 = ecnt2 = ecnt = ;
- memset(first1,-,sizeof(first1));
- memset(first2,-,sizeof(first2));
- memset(first,-,sizeof(first));
- }
- void Add_edge1(int u,int v,int c){
- nxt1[++ecnt1] = first1[u];
- e1[ecnt1].u = u;
- e1[ecnt1].v = v;
- e1[ecnt1].cost = c;
- e1[ecnt1].tag = ;
- e1[ecnt1].ge = ;
- first1[u] = ecnt1;
- }
- void Add_edge2(int u,int v,int c){
- nxt2[++ecnt2] = first2[u];
- e2[ecnt2].v = v;
- e2[ecnt2].cost = c;
- e2[ecnt2].tag = ;
- e2[ecnt2].ge = ;
- first2[u] = ecnt2;
- }
- void Add_edge(int u,int v,int c){
- nxt[ecnt] = first[u];
- e[ecnt].v = v;
- e[ecnt].u = u;
- e[ecnt].id = c;
- first[u] = ecnt++;
- nxt[ecnt] = first[v];
- e[ecnt].v = u;
- e[ecnt].u = v;
- e[ecnt].id = c;
- first[v] = ecnt++;
- }
- struct cmp{
- bool operator ()(pii a,pii b){
- return a.first > b.first;
- }
- };
- void Dijstra1(int s){
- priority_queue<pii,vector<pii >,cmp> PQ;
- dis1[s] = ;
- PQ.push(MP(dis1[s],s));
- int cnt = ;
- while(!PQ.empty()){
- pii x = PQ.top(); PQ.pop();
- if(dis1[x.second] < x.first) continue;
- for(int i = first1[x.second]; i != -; i = nxt1[i]){
- int v = e1[i].v;
- if(dis1[v] > dis1[x.second] + e1[i].cost){
- dis1[v] = dis1[x.second] + e1[i].cost;
- PQ.push(MP(dis1[v],v));
- }
- }
- }
- }
- void Dijstra2(int s){
- priority_queue<pii,vector<pii >,cmp> PQ;
- dis2[s] = ;
- PQ.push(MP(dis2[s],s));
- int cnt = ;
- while(!PQ.empty()){
- pii x = PQ.top(); PQ.pop();
- if(dis2[x.second] < x.first) continue;
- for(int i = first2[x.second]; i != -; i = nxt2[i]){
- int v = e2[i].v;
- if(dis2[v] > dis2[x.second] + e2[i].cost){
- dis2[v] = dis2[x.second] + e2[i].cost;
- PQ.push(MP(dis2[v],v));
- }
- }
- }
- }
- void Dfs(int p,int pre){
- dfn[p] = low[p] = ++tot;
- for(int i = first[p]; ~i; i = nxt[i]){
- int v = e[i].v;
- if(vis[i]) continue;
- vis[i] = vis[i ^ ] = true;
- if(!dfn[v]){
- Dfs(v,p);
- low[p] = min(low[p],low[v]);
- if(low[v] > dfn[p]) {
- bg[i] = bg[i ^ ] = true;
- ans[e[i].id] = ;
- }
- }
- else low[p] = min(low[p],dfn[v]);
- }
- }
- void Tarjan(){
- memset(dfn,,sizeof(dfn));
- memset(low,,sizeof(low));
- memset(bg,false,sizeof(bg));
- memset(vis,false,sizeof(vis));
- tot = ;
- for(int i = ; i <= N; ++i) if( dfn[i] == ) Dfs(i,-);
- }
- void solve(){
- memset(res,-,sizeof(res));
- for(int i = ;i <= M;i++){
- int u = e1[i].u;
- int v = e1[i].v;
- if(ans[i] == ) {
- res[i] = ;
- continue;
- }
- ll need = path - dis1[u] - dis2[v];
- if(need > ){
- res[i] = e1[i].cost - need + ;
- }
- }
- for(int i = ;i <= M;i++){
- if(res[i] == ) puts("YES");
- else if(res[i] == -) puts("NO");
- else printf("CAN %d\n",res[i]);
- }
- }
- int main(){
- int a,b,c,s,t;
- num = ;
- while(scanf("%d%d%d%d",&N,&M,&s,&t) != EOF){
- num++;
- init();
- int x;
- for(int i = ; i <= M; ++i){
- scanf("%d%d%d",&a,&b,&c);
- Add_edge1(a,b,c);
- Add_edge2(b,a,c);
- }
- for(int i = ;i <= N;i++) dis1[i] = dis2[i] = 1ll << ;
- Dijstra1(s);
- Dijstra2(t);
- path = dis1[t];
- for(int u = ;u <= N;u++){
- for(int i = first1[u];~i;i = nxt1[i]){
- int v = e1[i].v;
- if(dis1[u] + dis2[v] + e1[i].cost == path) {
- e1[i].tag = ;
- e2[i].tag = ;
- Add_edge(u,v,i);
- }
- }
- }
- memset(ans,,sizeof(ans));
- Tarjan();
- solve();
- }
- return ;
- }
codeforces 567 E. President and Roads 【 最短路 桥 】的更多相关文章
- Codeforces Round #Pi (Div. 2) E. President and Roads 最短路+桥
题目链接: http://codeforces.com/contest/567/problem/E 题意: 给你一个带重边的图,求三类边: 在最短路构成的DAG图中,哪些边是必须经过的: 其他的(包括 ...
- 【CodeForces 567E】President and Roads(最短路)
Description Berland has n cities, the capital is located in city s, and the historic home town of th ...
- Codeforces.567E.President and Roads(最短路 Dijkstra)
题目链接 \(Description\) 给定一张有向图,求哪些边一定在最短路上.对于不一定在最短路上的边,输出最少需要将其边权改变多少,才能使其一定在最短路上(边权必须为正,若仍不行输出NO). \ ...
- cf567E. President and Roads(最短路计数)
题意 题目链接 给出一张有向图,以及起点终点,判断每条边的状态: 是否一定在最短路上,是的话输出'YES' 如果不在最短路上,最少减去多少权值会使其在最短路上,如果减去后的权值\(< 1\),输 ...
- Codeforces Round #Pi (Div. 2) E. President and Roads tarjan+最短路
E. President and RoadsTime Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://codeforces.com/contest/567 ...
- Codeforces GYM 100876 J - Buying roads 题解
Codeforces GYM 100876 J - Buying roads 题解 才不是因为有了图床来测试一下呢,哼( 题意 给你\(N\)个点,\(M\)条带权边的无向图,选出\(K\)条边,使得 ...
- Codeforces Round #302 (Div. 2) D. Destroying Roads 最短路
题目链接: 题目 D. Destroying Roads time limit per test 2 seconds memory limit per test 256 megabytes input ...
- Codeforces Gym 100338C Important Roads 最短路+Tarjan找桥
原题链接:http://codeforces.com/gym/100338/attachments/download/2136/20062007-winter-petrozavodsk-camp-an ...
- Codeforces 543B Destroying Roads(最短路)
题意: 给定一个n个点(n<=3000)所有边长为1的图,求最多可以删掉多少条边后,图满足s1到t1的距离小于l1,s2到t2的距离小于l2. Solution: 首先可以分两种情况讨论: 1: ...
随机推荐
- [IOI2011]Race $O(nlog^{2}n)$ 做法
这个应该还不是正解,明天看一下正解到底是什么... Code: #include<bits/stdc++.h> #define setIO(s) freopen(s".in&qu ...
- AppScan用来进行漏洞扫描
https://www.cnblogs.com/mawenqiangios/p/8573525.html
- ASP组件AspJpeg(加水印)生成缩略图等使用方法
ASP组件AspJpeg(加水印)生成缩略图等使用方法 作者: 字体:[增加 减小] 类型:转载 时间:2012-12-17我要评论 ASPJPEG是一款功能相当强大的图象处理组件,用它可以轻松地做出 ...
- IOS - xib(Interface Builder,view) - can't change view size(view不能改变大小问题)
很多时候,我们自定义tableview.collectionview的cell,也有时候我们要自定义窗口xib,但创建xib后,其height.width不可修改. 这时问题就来了,怎么才能使我们的自 ...
- 简单JavaScript小程序
<!DOCTYPE html><html> <head> <meta charset="UTF-8"> ...
- 用Python图像处理
前几天弄了下django的图片上传,上传之后还需要做些简单的处理,python中PIL模块就是专门用来做这个事情的. 于是照葫芦画瓢做了几个常用图片操作,在这里记录下,以便备用. 这里有个字体文件,大 ...
- nodejs-NPM基本使用
搜索模块 npm search express 更新模块 npm update express 卸载模块 npm uninstall express 安装模块 npm install express ...
- jdk环境变量设置理解
1.系统变量→新建 JAVA_HOME 变量 . 变量值填写jdk的安装目录(本人是 E:\Java\jdk1.7.0) 2.系统变量→寻找 Path 变量→编辑 在变量值最后输入 %JAVA_HOM ...
- 关于心理的二十五种倾向(查理·芒格)-2
5)避免不一致倾向避免不一致倾向实际上就是人天生就害怕改变.相同是由于人类大脑的生理机制决定的.由于这样的倾向能够带来节省运算空间和能量的优点.这样的抗改变模式的形成,可能的原因例如以下:A) 迅速作 ...
- 16进制颜色字符串转为UIColor
//16进制颜色(html颜色值)字符串转为UIColor +(UIColor *) hexStringToColor: (NSString *) stringToConvert { NSS ...