AcWing 3797. 最大化最短路(排序优化+最短路)

传送门

思路:

首先从起点跟终点分别跑一遍最短路,假设 d i s 1 [ i ] dis1[i] dis1[i]表示从 1 1 1出发到 i i i的最短路, d i s 2 [ i ] dis2[i] dis2[i]表示从 n n n出发到 i i i的最短路。
假设中间的两个特殊点为 a , b a,b a,b,在 a , b a,b a,b连边后新的最短路为 d i s 1 [ a ] + d i s 2 [ b ] + 1 dis1[a]+dis2[b]+1 dis1[a]+dis2[b]+1
如果暴力枚举特殊点的话,时间复杂度就会变成 O ( k 2 ) O(k^2) O(k2),显然不可取。
比赛时就卡这了,一直没想到优化的方法,赛后看了巨巨的代码才懂。
正常应该有两种类型的最短路:
1 − > a − > b − > n 1->a->b->n 1>a>b>n
1 − > b − > a − > n 1->b->a->n 1>b>a>n
假设最短路为第一种情况,那么会有:
d i s 1 [ a ] + d i s 2 [ b ] + 1 < d i s 1 [ b ] + 1 + d i s 2 [ a ] dis1[a]+dis2[b]+1<dis1[b]+1+dis2[a] dis1[a]+dis2[b]+1<dis1[b]+1+dis2[a]
移项与消掉重复项后得:
d i s 1 [ a ] − d i s 2 [ a ] < d i s 1 [ b ] − d i s 2 [ b ] dis1[a]-dis2[a]<dis1[b]-dis2[b] dis1[a]dis2[a]<dis1[b]dis2[b]
所以可以按这个排序后,枚举 b b b并且维护最大的dis1[a],取最大值就好了,时间复杂度降为 O ( k l o g k ) O(klogk) O(klogk).
要注意的是,向图中加边不会让最短路增加,所以说如果最大值大于原先的最短路,答案还是原先的最短路。

代码:

const int maxn=2e5+10,N=maxn*2,inf=0x3f3f3f3f;

int n,m,k,a[maxn],dis1[maxn],dis2[maxn],st[maxn];

int h[N], w[N], e[N], ne[N], idx;
int dist[N];

struct node{
	int id,ds,dt;
}b[maxn];

bool cmp(node a,node b){
	return a.ds-a.dt<b.ds-b.dt;
}

void add(int a, int b, int c)
{
    e[idx] = b, w[idx] = c, ne[idx] = h[a], h[a] = idx ++ ;
}

void spfa(int s,int dist[])
{
    memset(st,0,sizeof st);
    dist[s] = 0;
    queue<int> q;
    q.push(s);
    st[s] = true;

    while (q.size())
    {
        int t = q.front();
        q.pop();

        st[t] = false;

        for (int i = h[t]; i != -1; i = ne[i])
        {
            int j = e[i];
            if (dist[j] > dist[t] + w[i])
            {
                dist[j] = dist[t] + w[i];
                if (!st[j])
                {
                    q.push(j);
                    st[j] = true;
                }
            }
        }
    }
}

int main(){
	memset(h,-1,sizeof h);
	memset(dis1,0x3f,sizeof dis1);
	memset(dis2,0x3f,sizeof dis2);
	n=read,m=read,k=read;
	rep(i,1,k) a[i]=read;
	rep(i,1,m){
		int x=read,y=read;
		add(x,y,1);add(y,x,1);
	}
	spfa(1,dis1);spfa(n,dis2);
	rep(i,1,k){
		b[i]={a[i],dis1[a[i]],dis2[a[i]]};
	}
	sort(b+1,b+1+k,cmp);
	int res=0;
	int las=b[1].ds;
	rep(i,2,k){
		res=max(res,las+1+b[i].dt);
		las=max(las,b[i].ds);
	}
	if(res>dis1[n]) res=dis1[n];
	write(res);
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

豆沙睡不醒

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值