2018.10.23【校内模拟】“新”的家园(最短路)(玄学建图)

传送门


解析:

首先我们是不能直接在原图上直接跑最短路的,时间肯定爆炸。(废话)

我们注意一个性质,环外边最多只有300,第一种情况不管,反正只有一组询问。

那么这三百条边最多会连接600个节点。这些点我们称为关键点,这些边称为关键边。
可以发现,除了环上距离,能够更新两点之间的最短路的只有关键边。

那么考虑原来的环上相邻两个关键点之间的那些点。显然,当所求最短路的起始点不在这些点中时,这些点的更新只会拖慢求最短路的进度,就是拖延收敛时间。

那么我们考虑把这些点去掉,只保留关键点,保留关键点之间的关建边,同时关键点之间连权值为原来环上距离的边。

那么我们每次就只在简略后的图上面跑最短路,只需要把我们要求的 &lt; s , t &gt; &lt;s,t&gt; <s,t>加回原图,与原来相邻的关键点连距离为环上距离的边。显然,新图上的距离仍然符合原图的情况,求出来的最短路显然也是原图的最短路。

注意一点,我们必须在求完最短路后删去添加的点和边,不然复杂度显然会逐渐增大,直到最后无法承受。

删边其实很简单,我们只需要用前向星记录一下边,新增的边显然在前向星的末尾,直接删去就好了。

注意求的 &lt; s , t &gt; &lt;s,t&gt; <s,t>有可能本来就是不需要经过关键点的最短路,所以加边的时候需要直接加一条 s , t s,t s,t的环上距离的边。

还有就是其实点数很小,所以我们可以承受跑那么多次最短路,无所畏惧。

而且这是一个 n , m n,m n,m同阶的稀疏图,所以 S P F A SPFA SPFA不会被卡,且表现优于 D i j k s t r a Dijkstra Dijkstra

最后一点,我们不能很方便的维护环上距离,所以直接断环成链,强行令 &lt; 1 , n &gt; &lt;1,n&gt; <1,n>这条边为特殊边。维护所有节点到 1 1 1的链上距离,其他的直接差分就好了。


代码:

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define re register
#define gc getchar
#define pc putchar
#define cs const

inline int getint(){
	re int num;
	re char c;
	while(!isdigit(c=gc()));num=c^48;
	while(isdigit(c=gc()))num=(num<<1)+(num<<3)+(c^48);
	return num;
}

inline void outint(int a){
	static char ch[23];
	if(a==0)pc('0');
	while(a)ch[++ch[0]]=a-a/10*10,a/=10;
	while(ch[0])pc(ch[ch[0]--]^48);
}

cs int N=500005,M=2000006,INF=0x3f3f3f3f;
bool iskey[N];
int head[N];
int last[N],nxt[M<<1],to[M<<1],ecnt;
int w[M<<1];
inline void addedge(int u,int v,int val){
	assert(val>=0);
	nxt[++ecnt]=last[u],last[u]=ecnt,to[ecnt]=v,w[ecnt]=val;
	nxt[++ecnt]=last[v],last[v]=ecnt,to[ecnt]=u,w[ecnt]=val;
}

int top,sta[N];
int Ecnt;
inline void clear(){
	ecnt=Ecnt;
	while(top){
		last[sta[top]]=head[sta[top]];
		--top;
	}
}

int pt[N],pcnt;
int sum[N];
inline void ins(int s){
	if(iskey[s])return ;
	for(int re i=1;i<=pcnt;++i){
		if(pt[i]<s&&pt[i+1]>s){
			sta[++top]=s;
			sta[++top]=pt[i];
			sta[++top]=pt[i+1];
			addedge(s,pt[i],sum[s]-sum[pt[i]]);
			addedge(s,pt[i+1],sum[pt[i+1]]-sum[s]);
			return ;
		}
	}
}

bool vis[N];
int dist[N];
queue<int > q;
inline int SPFA(int s,int t){
	if(s==t)return 0;
	for(int re i=1;i<=pcnt;++i)dist[pt[i]]=INF;
	dist[s]=0;dist[t]=INF;
	q.push(s);
	while(!q.empty()){
		int u=q.front();
		q.pop();
		vis[u]=false;
		for(int re e=last[u],v=to[e];e;v=to[e=nxt[e]]){
			if(dist[v]>dist[u]+w[e]){
				dist[v]=dist[u]+w[e];
				if(vis[v])continue;
				vis[v]=true;
				q.push(v);
			}
		}
	}
	return dist[t];
}

int n,m,Q;
signed main(){
	memset(sum,0x3f,sizeof sum);
	n=getint();
	m=getint();
	Q=getint();
	for(int re i=1;i<=m;++i){
		int u=getint(),v=getint(),val=getint();
		if(u>v)swap(u,v);
		if(u+1!=v){
			iskey[u]=iskey[v]=true;
			addedge(u,v,val);
		}
		else sum[v]=min(sum[v],val);
	}
	sum[0]=sum[1]=0;
	for(int re i=1;i<=n;++i)sum[i]+=sum[i-1];
	for(int re i=1;i<=n;++i){
		static int last=0;
		if(iskey[i]){
			pt[++pcnt]=i;
			if(!last){
				last=i;
				continue;
			}
			addedge(last,i,sum[i]-sum[last]);
			last=i;
		}
	}
	Ecnt=ecnt;
	memcpy(head,last,sizeof last);
	while(Q--){
		clear();
		int s=getint(),t=getint();
		if(s>t)swap(s,t);
		addedge(s,t,sum[t]-sum[s]);
		sta[++top]=s;sta[++top]=t;
		ins(s);ins(t);
		outint(SPFA(s,t));
		pc('\n');
	}
	return 0;
}
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值