Warfare And Logistics UVALive

Warfare And Logistics

 UVALive - 4080 

题意：给n个点m条边。令c为每对节点的最短路长度之和。要求删除一条边后使得新的c值c‘最大，不连通的两点对短路视为L。

【如果每次删除一条边，要跑m次dijkstra，其实其中很多次都对最短路没有影响，因为删掉的边不在最短路里】

【因此，可以每次删除最短路树中的一条边，需要跑n次，复杂度降低到可接受程度】

#include <bits/stdc++.h>
#define LL long long
using namespace std;
const int inf=0x3f3f3f3f;
const int maxv=110;
const int maxe=2010;
int n,m,L;
struct Edge
{
    int u,v,w;
    int nex;
}e[maxe<<1];
int head[maxv];
int cnt=0;
void init()
{
    memset(head,-1,sizeof(head));
    cnt=0;
}
void add(int u,int v,int w)
{
    e[cnt].u=u;
    e[cnt].v=v;
    e[cnt].w=w;
    e[cnt].nex=head[u];
    head[u]=cnt++;
}

typedef pair<int,int> PII;
int dis[maxv],par[maxv];
void dijkstra(int s,bool c,int e1)
{
    priority_queue<PII,vector<PII>,greater<PII> > pq;
    for(int i=0;i<maxv;i++) dis[i]=L;
    dis[s]=0;
    pq.push(PII(0,s));
    par[s]=-1;
    while(!pq.empty())
    {
        PII t=pq.top();
        pq.pop();
        int u=t.second;
        if(dis[u]<t.first) continue;
        for(int i=head[u];i!=-1;i=e[i].nex)
        {
            if(i==e1||i==(e1^1)) continue;   //
            int v=e[i].v;
            if(dis[v]>dis[u]+e[i].w)
            {
                dis[v]=dis[u]+e[i].w;
                if(c) par[v]=i;
                pq.push(PII(dis[v],v));
            }
        }
    }
}
LL cal()
{
    LL temp=0;
    for(int i=0;i<n;i++)
        temp+=dis[i];
    return temp;
}
vector<int> vv;
int delta[maxe];

int main()
{
    while(scanf("%d%d%d",&n,&m,&L)!=EOF)
    {
        init();
        int u,v,w;
        for(int i=0;i<m;i++)
        {
            scanf("%d%d%d",&u,&v,&w);
            u--;v--;
            add(u,v,w);
            add(v,u,w);
        }
        LL ans=0;
        vv.clear();
        memset(delta,0,sizeof(delta));
        for(int i=0;i<n;i++)
        {
            dijkstra(i,1,-1);
            LL temp=cal();
            ans+=temp;
            for(int j=0;j<n;j++)
                if(par[j]!=-1)
            {
                dijkstra(i,0,par[j]);
                int id=par[j]&(~1);   //
                if(!delta[id]) vv.push_back(id);
                delta[id]+=cal()-temp;
            }

        }
        int maxdelta=0;
        for(int i=0;i<(int)vv.size();i++)
            maxdelta=max(maxdelta,delta[vv[i]]);
        printf("%lld %lld
",ans,ans+maxdelta);
    }
}