首页 > 其他 > 详细

Codeforces Round #257(Div.2) D Jzzhu and Cities --SPFA

时间:2014-07-22 22:48:43      阅读:369      评论:0      收藏:0      [点我收藏+]

题意:n个城市,中间有m条道路(双向),再给出k条铁路,铁路直接从点1到点v,现在要拆掉一些铁路,在保证不影响每个点的最短距离(距离1)不变的情况下,问最多能删除多少条铁路

分析:先求一次最短路,铁路的权值大于该点最短距离的显然可以删去,否则将该条边加入图中,再求最短路,记录每个点的前一个点,然后又枚举铁路,已经删去的就不用处理了,如果铁路权值大于该点最短距离又可以删去,权值相等时,该点的前一个点如果不为1,则这个点可以由其他路到达,这条铁路又可以删去。

由于本题中边比较多,最多可以有8x10^5条,所以用邻接链表会超时,最好用vector来存储边。

 

Time:296 ms Memory:30968 KB

bubuko.com,布布扣
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <cstdlib>
#include <algorithm>
#include <vector>
#include <queue>
#define lll __int64
using namespace std;
#define N 100006

struct Edge
{
    int v;
    lll w;
    Edge(int _v,lll _w)
    {
        v = _v;
        w = _w;
    }
    Edge(){}
};

struct Train
{
    int v;
    lll w;
}T[N];

vector<Edge> G[N];
lll dis[N];
int n,m;
int head[N],tot,pre[N];
int cut[N],vis[N];
queue<int> que;

void SPFA(int s)
{
    while(!que.empty())
        que.pop();
    memset(vis,0,sizeof(vis));
    vis[s] = 1;
    que.push(s);
    for(int i=0;i<=n;i++)
        dis[i] = 1000000000000000LL;
    dis[s] = 0;
    while(!que.empty())
    {
        int u = que.front();
        que.pop();
        vis[u] = 0;
        for(int i=0;i<G[u].size();i++)
        {
            int v = G[u][i].v;
            lll w = G[u][i].w;
            if(dis[v] > dis[u] + w)
            {
                dis[v] = dis[u] + w;
                if(!vis[v])
                {
                    vis[v] = 1;
                    que.push(v);
                }
            }
        }
    }
}

void SPFA2()
{
    while(!que.empty())
    {
        int u = que.front();
        que.pop();
        vis[u] = 0;
        for(int i=0;i<G[u].size();i++)
        {
            int v = G[u][i].v;
            lll w = G[u][i].w;
            if(dis[v] > dis[u] + w)
            {
                dis[v] = dis[u] + w;
                pre[v] = u;
                if(!vis[v])
                {
                    que.push(v);
                    vis[v] = 1;
                }
            }
            else if(dis[v] == dis[u] + w && pre[v] < u)
                pre[v] = u;
        }
    }
}

int main()
{
    int i,j,k;
    int u,v,y;
    lll w;
    while(scanf("%d%d%d",&n,&m,&k)!=EOF)
    {
        tot = 0;
        for(i=0;i<=n;i++)
            G[i].clear();
        memset(head,-1,sizeof(head));
        memset(cut,0,sizeof(cut));
        memset(pre,-1,sizeof(pre));
        for(i=0;i<m;i++)
        {
            scanf("%d%d%I64d",&u,&v,&w);
            G[u].push_back(Edge(v,w));
            G[v].push_back(Edge(u,w));
        }
        for(i=0;i<k;i++)
        {
            scanf("%d%I64d",&y,&w);
            T[i].v = y;
            T[i].w = w;
        }
        SPFA(1);
        int cnt = 0;
        for(i=0;i<k;i++)
        {
            int v = T[i].v;
            lll w = T[i].w;
            if(dis[v] <= w)
            {
                cut[i] = 1;
                cnt++;
            }
            else
            {
                G[1].push_back(Edge(v,w));
                G[v].push_back(Edge(1,w));
                pre[v] = 1;
                dis[v] = w;
                que.push(v);
                vis[v] = 1;
            }
        }
        SPFA2();
        for(i=0;i<k;i++)
        {
            if(cut[i])
                continue;
            int v = T[i].v;
            lll w = T[i].w;
            if((dis[v] == w && pre[v] != 1) || dis[v] < w)
                cnt++;
        }
        printf("%d\n",cnt);
    }
    return 0;
}
View Code

Codeforces Round #257(Div.2) D Jzzhu and Cities --SPFA,布布扣,bubuko.com

Codeforces Round #257(Div.2) D Jzzhu and Cities --SPFA

原文:http://www.cnblogs.com/whatbeg/p/3855921.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!