CodeForces 567E President and Roads(最短路 + tarjan求桥)

//最短路+桥
#include<cstdio>
#include<cstring>
#include<queue>
#include<algorithm>
using namespace std;
typedef long long LL;
const int M=100100,N=100010;
const LL INF=0x7fffffffffffff;
struct Edge{
    int from,to,cap,next;
    int id;
}edge[2*M];
int head[2][N];
int cnt;
struct road
{
    int u,v,w;
    int flag;
    LL ans;
}R[M];
void init(){
    cnt=0;
    memset(head, -1, sizeof(head));
}
void addedge(int u,int v,int w,int ii,int flag){
    edge[cnt].from=u;edge[cnt].to=v;edge[cnt].cap=w;edge[cnt].id=ii;
    edge[cnt].next=head[flag][u];head[flag][u]=cnt++;
}
struct HeapNode{
    LL dis;
    int id;
    bool operator<(const HeapNode &A)const{
        return dis>A.dis;
    }
};
LL dis[2][N];
int vis[N];
void Dijstra(int s,int flag,int n){
    priority_queue<HeapNode>Q;
    for(int i=1;i<=n;i++){
        dis[flag][i]=INF;
        vis[i]=0;
    }
    dis[flag][s]=0;
    Q.push((HeapNode){0,s});
    while(!Q.empty()){
        HeapNode k=Q.top();
        Q.pop();
        int u=k.id;
        if(vis[u])continue;
        vis[u]=1;
        for(int i=head[flag][u];i!=-1;i=edge[i].next){
            int v=edge[i].to;
            if(dis[flag][u]+edge[i].cap<dis[flag][v]){
                dis[flag][v]=dis[flag][u]+edge[i].cap;
                Q.push((HeapNode){dis[flag][v],v});
            }
        }
    }
}
void Reduce(int m,int t){
    LL mincost=dis[0][t];
    init();
    for(int i=0;i<m;i++){
        int u=R[i].u,v=R[i].v;
        if(dis[0][u]==INF||dis[1][v]==INF){
            R[i].flag=-2;continue;
        }
        if(dis[0][u]+dis[1][v]+(LL)R[i].w>mincost){
            R[i].flag=-1;
            R[i].ans=(LL)R[i].w+1-(mincost-dis[0][u]-dis[1][v]);
        }
        else if(dis[0][u]+dis[1][v]+(LL)R[i].w==mincost){
            addedge(u,v,R[i].w,i,0);
            addedge(v,u,R[i].w,i,0);
        }
    }
}
int dfn[N],low[N];
int indx;
int fa[N];//fa[v]记录的是来边的编号
void tarjan(int u){//重边求桥
    dfn[u]=low[u]=++indx;
    for(int i=head[0][u];i!=-1;i=edge[i].next){
        int v=edge[i].to;
        if(dfn[v]==-1){
            fa[v]=edge[i].id;
            tarjan(v);
            low[u]=min(low[u],low[v]);
            if(low[v]>dfn[u]){//子节点到不了u及以上
                R[edge[i].id].flag=1;
            }
        }
        else {
            if(fa[u]!=edge[i].id){
                low[u]=min(low[u],dfn[v]);
            }
        }
    }
}
void output(int m){
    for(int i=0;i<m;i++){
        if(R[i].flag==-1){
            if(R[i].ans<R[i].w&&R[i].ans>0)
                printf("CAN %lld\n",R[i].ans);
            else
                printf("NO\n");
        }
        else if(R[i].flag==1){
            printf("YES\n");
        }
        else {
            if(R[i].w>1&&R[i].flag!=-2)
                printf("CAN 1\n");
            else
                printf("NO\n");
        }
    }
}
int main(){
    int n,m,s,t;
    init();
    scanf("%d%d%d%d",&n,&m,&s,&t);
    for(int i=0;i<m;i++){
        scanf("%d%d%d",&R[i].u,&R[i].v,&R[i].w);
        addedge(R[i].u,R[i].v,R[i].w,i,0);
        addedge(R[i].v,R[i].u,R[i].w,i,1);
    }
    Dijstra(s,0,n);
    Dijstra(t,1,n);
    Reduce(m,t);
    memset(dfn, -1, sizeof(dfn));
    memset(low, 0, sizeof(low));
    indx=0;
    tarjan(s);
    output(m);
    return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值