CF1422D Returning Home -- 最短路+贪心连边

题面

CodeForces - 1422D 

题解

关于题目翻译，去百度吧，看得懂的。

先读懂题：“instant-movement locations”（简称IML）的作用是把你瞬间（即耗时为零）从同一行或同一列召到这个位置，

考虑把起点、终点和每个IML之间连边跑最短路，由于IML的特殊性质，任何一个点与一个IML之间的距离为x坐标之差与y坐标之差的最小值；

由于你可以从同一行或同一列瞬间到达某个IML，但不能从某个IML位置出发瞬间到达同一行或同一列的某个位置，所以从IML连出去的时候得把它看成普通点。

暴力连边是的，时间空间都会炸；

考虑如下的情况：

如果把它们两两之间连边，从a到c的距离是，但如果只连a-b、b-c，从a到c的距离仍是，

所以我们只需要连每个IML上方或下方y坐标离它最近的点，横向同理，只连x坐标最近的两个点。

这样最多连4m条边，堆优化Dij，线段树优化

代码

我用的是线段树优化的Dij，其实大可不必...

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<queue>
#include<vector>
#define ll long long
#define MAXN 100005
#define MOD 998344353ll
using namespace std;
inline ll read(){
	ll x=0;bool f=1;char s=getchar();
	while((s<'0'||s>'9')&&s>0){if(s=='-')f^=1;s=getchar();}
	while(s>='0'&&s<='9')x=(x<<1)+(x<<3)+s-'0',s=getchar();
	return f?x:-x;
}
int n,m,f[MAXN*3],p,st[MAXN];
bool vis[MAXN];
ll d[MAXN];

struct P{ll x,y;}s[MAXN];
struct itn{
	int v;ll w;
	itn(){}
	itn(int V,ll W){v=V,w=W;}
};
inline bool cmp(P a,P b){
	if(a.x!=b.x)return a.x<b.x;
	else return a.y<b.y;
}
inline bool cmp_(int a,int b){return s[a].y<s[b].y;}
vector<itn>G[MAXN];

inline ll ads(ll x){return x>0?x:-x;}
inline ll getd(P a,P b){return min(ads(a.x-b.x),ads(a.y-b.y));}
inline ll getd_(P a,P b){return ads(a.x-b.x)+ads(a.y-b.y);}
inline int MIN(int a,int b){return d[a]<d[b]?a:b;}

inline void change(int x,int d){//四行zkw
	for(f[p+x]=d,x=(x+p)>>1;x>0;x>>=1)
		f[x]=MIN(f[x<<1],f[x<<1|1]);
}

int main()
{
	n=read(),m=read();
	s[0].x=read(),s[0].y=read(),s[m+1].x=read(),s[m+1].y=read();
	for(int i=1;i<=m;i++)s[i].x=read(),s[i].y=read();
	sort(s+1,s+1+m,cmp);
	for(int i=1;i<=m;i++)st[i]=i;
	sort(st+1,st+1+m,cmp_);
	G[0].push_back(itn(m+1,getd_(s[0],s[m+1])));
	for(int i=1;i<m;i++){
		G[i].push_back(itn(i+1,getd(s[i],s[i+1])));
		G[i+1].push_back(itn(i,getd(s[i],s[i+1])));
	}
	for(int i=1;i<m;i++){
		G[st[i]].push_back(itn(st[i+1],getd(s[st[i]],s[st[i+1]])));
		G[st[i+1]].push_back(itn(st[i],getd(s[st[i]],s[st[i+1]])));
	}
	for(int i=1;i<=m;i++){
		G[0].push_back(itn(i,getd(s[0],s[i])));
		G[i].push_back(itn(m+1,getd_(s[i],s[m+1])));
	}
	memset(d,0x7f,sizeof(d));
	for(p=1;p<m+3;p<<=1);
	for(int i=0;i<G[0].size();i++){
		int v=G[0][i].v;ll w=G[0][i].w;
		d[v]=min(d[v],w),change(v,v);
	}
	for(int T=1;T<=m;T++){
		int u=f[1];//不用查询，直接问祖
		if(u==0)break;
		change(u,0),vis[u]=1;
		for(int i=0;i<G[u].size();i++)
			if(!vis[G[u][i].v]){
				int v=G[u][i].v;ll w=G[u][i].w;
				d[v]=min(d[v],d[u]+w),change(v,v);
			}
	}
	printf("%lld\n",d[m+1]);
	return 0;
}