codeforces 544D Destroying Roads (最短路)

题目大意：n个点m条边的无向图（无重边），给出两个起点s1,s2，两个终点t1,t2，现在要尽可能多地删除图中的边，使得s1->t1路径长度不超过l1，s2->t2路径长度不超过l2。

分析：显然除了s1->t1的最短路径与s2->t2的最短路径保留外，其余边都删掉，删除的边能够达到最多。

以每个顶点为起点各跑一次最短路。求出各个顶点到其它顶点的最短距离。由于边长均为1，故路径长度即为边数。

判断s1->t1与s2->t2的最短距离是否分别小于l1、l2。若不是，则无解，输出-1.

然后先让结果为两条路径边数之和。

暴力枚举重复的路径，更新结果，时间复杂度为O(nm+n*n)

注意枚举的时候s1,s2,t1,t2的相对位置有四种情况，如下图所示：

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
#define N 3005
#define inf 0x3f3f3f3f
struct Edge
{
    int to,next,w;
}edge[N*N];

int cnt,n,head[N],low[N][N],w[N];
bool vis[N];

inline void add(int u,int v,int w)
{
    edge[cnt].to=v;
    edge[cnt].w=w;
    edge[cnt].next=head[u];
    head[u]=cnt++;
}

void spfa()
{
    int i,j;
    for(j=1; j<=n; ++j)
    {
        memset(vis,0,sizeof(vis));
        low[j][j]=0;
        vis[j]=1;
        queue<int> q;
        q.push(j);
        while(!q.empty())
        {
            int u=q.front();
            q.pop();
            vis[u]=0;
            for(i=head[u]; ~i; i=edge[i].next)
            {
                int v=edge[i].to;
                if(low[j][v]>low[j][u]+edge[i].w)
                {
                    low[j][v]=low[j][u]+edge[i].w;
                    if(!vis[v])
                    {
                        vis[v]=1;
                        q.push(v);
                    }
                }
            }
        }
    }
}

int main()
{
    int m,i,j,x,y,s1,t1,s2,t2,l1,l2;
    cin>>n>>m;
    memset(head,-1,sizeof(head));
    memset(low,inf,sizeof(low));
    for(i=1;i<=m;++i){
        scanf("%d%d",&x,&y);
        add(x,y,1);
        add(y,x,1);
    }
    spfa();
    scanf("%d%d%d%d%d%d",&s1,&t1,&l1,&s2,&t2,&l2);
    if(low[s1][t1]>l1||low[s2][t2]>l2) puts("-1");
    else{
        int ans=low[s1][t1]+low[s2][t2];
        for(i=1;i<=n;++i)
            for(j=1;j<=n;++j)
            {
                if(low[s1][i]+low[i][j]+low[j][t1]<=l1&&low[s2][i]+low[i][j]+low[j][t2]<=l2)
                    ans=min(ans,low[s1][i]+low[i][j]+low[j][t1]+low[s2][i]+low[i][j]+low[j][t2]);
                if(low[s1][i]+low[i][j]+low[j][t1]<=l1&&low[t2][i]+low[i][j]+low[j][s2]<=l2)
                    ans=min(ans,low[s1][i]+low[i][j]+low[j][t1]+low[t2][i]+low[i][j]+low[j][s2]);
                if(low[t1][i]+low[i][j]+low[j][s1]<=l1&&low[s2][i]+low[i][j]+low[j][t2]<=l2)
                    ans=min(ans,low[t1][i]+low[i][j]+low[j][s1]+low[s2][i]+low[j][t2]);
                if(low[t1][i]+low[i][j]+low[j][s1]<=l1&&low[t2][i]+low[i][j]+low[j][s2]<=l2)
                    ans=min(ans,low[t1][i]+low[i][j]+low[j][s1]+low[t2][i]+low[j][s2]);
            }
        printf("%d\n",m-ans);
    }
    return 0;
}