SPFA算法有两种优化算法 SLF 和 LLL
SLF:Small Label First 大概是双向队列?
如果要加入的元素v的距离dis[v]比当前的队首元素dis[q[head]]小的话,那么就把当前元素插到队首。
如果是手打队列的话,就先head-2,再q[head+1]=v;
LLL:Large Label Last
如果队列q的队首元素为x,队列中所有元素的平均距离为dist,如果dis[x]>dist,则将它移动到末尾,如果有一个元素的dis[x]<=dist,则将它出队。
贴个板:
(两个优化都有加)
(以bzoj 4152 为例 题面就不贴了)
但这道题跟SPFA有仇,SPFA应该是过不了的,要用dijkstra+堆优化
//这道题是有多看不起SPFA
//SLF && LLL
//然而有什么用呢 还是T掉了
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ms(x,y) memset(x,y,sizeof(x))
#include<queue>
#define ll long long
using namespace std;
const int N = 200000 + 10;
const int M = N*4;
inline ll read()
{
ll data=0,w=1; char ch=0;
while(ch!='-' && (ch<'0' || ch>'9')) ch=getchar();
if(ch=='-') w=-1,ch=getchar();
while(ch>='0' && ch<='9') data=data*10+ch-'0',ch=getchar();
return data*w;
}
int n;
deque<int> q;
struct Node{
ll x,y;
int i;
}a[N];
struct node{
int pre,v;
ll w;
}edge[M];
int num=0;
int head[N];
void addedge(int from,int to,ll w){
num++;
edge[num].pre=head[from];
edge[num].v=to;
edge[num].w=w;
head[from]=num;
}
bool cmp1(const Node &a,const Node &b){
return a.x<b.x;
}
bool cmp2(const Node &a,const Node &b){
return a.y<b.y;
}
void getmap(){
sort(a+1,a+n+1,cmp1);
for(register int i=2;i<=n;i++){
ll w=a[i].x-a[i-1].x;
addedge(a[i-1].i,a[i].i,w);
addedge(a[i].i,a[i-1].i,w);
}
sort(a+1,a+n+1,cmp2);
for(register int i=2;i<=n;i++){
ll w=a[i].y-a[i-1].y;
addedge(a[i-1].i,a[i].i,w);
addedge(a[i].i,a[i-1].i,w);
}
}
ll dis[N];
bool exist[N];
int tot;
ll sum=0;
void SPFA(){
ms(dis,127);ms(exist,0);
q.push_back(1),dis[1]=0,exist[1]=true,tot=1;
while(!q.empty()){
int u=q.front();q.pop_front();
if(dis[u]*tot>sum){
q.push_back(u);
continue;
}
exist[u]=false;tot--,sum -= dis[u];
for(int i=head[u];i;i=edge[i].pre){
int v=edge[i].v;
if(dis[v]>dis[u]+edge[i].w){
dis[v]=dis[u]+edge[i].w;
if(!exist[v]){
exist[v]=true;
if(!q.empty()&&dis[v]<dis[q.front()]) q.push_front( v ) ;
else q.push_back(v);
tot++,sum+=dis[v];
}
}
}
}
}
int main(){
scanf("%d",&n);
for(register int i=1;i<=n;i++){
a[i].x=read(),a[i].y=read();
a[i].i=i;
}
getmap();
/* for(int i=head[1];i;i=edge[i].pre){
int v=edge[i].v;
printf("%d %d\n",v,edge[i].w);
}*/
SPFA();
printf("%lld",dis[n]);
return 0;
}