传送门
解析:
首先我们是不能直接在原图上直接跑最短路的,时间肯定爆炸。(废话)
我们注意一个性质,环外边最多只有300,第一种情况不管,反正只有一组询问。
那么这三百条边最多会连接600个节点。这些点我们称为关键点,这些边称为关键边。
可以发现,除了环上距离,能够更新两点之间的最短路的只有关键边。
那么考虑原来的环上相邻两个关键点之间的那些点。显然,当所求最短路的起始点不在这些点中时,这些点的更新只会拖慢求最短路的进度,就是拖延收敛时间。
那么我们考虑把这些点去掉,只保留关键点,保留关键点之间的关建边,同时关键点之间连权值为原来环上距离的边。
那么我们每次就只在简略后的图上面跑最短路,只需要把我们要求的 < s , t > <s,t> <s,t>加回原图,与原来相邻的关键点连距离为环上距离的边。显然,新图上的距离仍然符合原图的情况,求出来的最短路显然也是原图的最短路。
注意一点,我们必须在求完最短路后删去添加的点和边,不然复杂度显然会逐渐增大,直到最后无法承受。
删边其实很简单,我们只需要用前向星记录一下边,新增的边显然在前向星的末尾,直接删去就好了。
注意求的 < s , t > <s,t> <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。
最后一点,我们不能很方便的维护环上距离,所以直接断环成链,强行令 < 1 , n > <1,n> <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;
}