题目询问 $ s $ 到 $ t $ 的最短距离，我们可以发现我们就没法跑记录二维是否用过传送的方法。
可以发现 $ s->t $ 的最短距离可以看做 $ s->x->y->t $ ，$ x,y $ 为随便选的两个点，设 $ dis1[i] $ 为 $ s $ 到 $ i $ 的最短距离 ，$ dis2[i] $ 为 $ t $ 到 $ i $ 的最短距离 ,那么答案就是 $ \min\limits_{1\leq x \leq n,1\leq y \leq n} dis1[x]+(x-y)^2+dis2[y] $ 。
建两个图，一个正图，一个反图，跑两遍最短路求出 $ dis1,dis2 $ ，然后 $ n^2 $ 枚举每一个点，我们就有了 $ O(n^2) $ 的做法，考虑优化。
观察答案的式子，我们可以给它化简成 \(dis1[x]+x^2-2xy+dis2[y]+y^2\) 的形式，因为我们 $ O(n) $ 枚举 $ x $ ，所以每次枚举 $ dis1[x]+x^2 $ 都是一个定值，我们只要求 $ dis2[y]+y^2-2xy $ 的最小值就好，把每一个 $ y $ 当成一条线段，那么它的斜率就是 $ -2x $，b就是 $ dis2[y]+y^2 $ ，我们就可以直接套一个李超跑过去。
点击查看代码
#include<bits/stdc++.h>
#define int long long
#define ls(x) x<<1
#define rs(x) x<<1|1
using namespace std;
constexpr int N=2e5+10;
int head1[N],cnt1,head2[N],cnt2,dis1[N],dis2[N],n,m,S,T,vis[N],s[N<<2],k[N],b[N];
struct edge{int next,to,w;}e1[N<<1],e2[N<<1];
inline int in(){int k=0,f=1;char c=getchar();while(c<'0'||c>'9'){if(c=='-') f=-1;c=getchar();}while(c>='0'&&c<='9') k=(k<<3)+(k<<1)+c-'0',c=getchar();return k*f;
}
void add1(int u,int v,int w){e1[++cnt1].next=head1[u],e1[cnt1].to=v,e1[cnt1].w=w,head1[u]=cnt1;
}
void add2(int u,int v,int w){e2[++cnt2].next=head2[u],e2[cnt2].to=v,e2[cnt2].w=w,head2[u]=cnt2;
}
struct node{int val,u;};
inline int operator>(node a,node b){return a.val>b.val;};
void dij1(int x){memset(dis1,0x3f,sizeof(dis1));memset(vis,0,sizeof(vis));dis1[x]=0;priority_queue<node,vector<node>,greater<node> > q;q.push({0,x});while(!q.empty()){node u=q.top();q.pop();if(vis[u.u]) continue;vis[u.u]=1;for(int i=head1[u.u];i;i=e1[i].next){int v=e1[i].to;if(dis1[v]>dis1[u.u]+e1[i].w){dis1[v]=dis1[u.u]+e1[i].w;if(!vis[v]) q.push({dis1[v],v});}}}
}
void dij2(int x){memset(dis2,0x3f,sizeof(dis2));memset(vis,0,sizeof(vis));dis2[x]=0;priority_queue<node,vector<node>,greater<node> > q;q.push({0,x});while(!q.empty()){node u=q.top();q.pop();if(vis[u.u]) continue;vis[u.u]=1;for(int i=head2[u.u];i;i=e2[i].next){int v=e2[i].to;if(dis2[v]>dis2[u.u]+e2[i].w){dis2[v]=dis2[u.u]+e2[i].w;if(!vis[v]) q.push({dis2[v],v});}}}
}
inline int calc(int id,int x){return k[id]*x+b[id];
}
void ins(int x,int l,int r,int u){if(!s[x]){s[x]=u;return;}if(l==r){if(calc(s[x],l)>calc(u,l))s[x]=u;return;}int mid=(l+r)>>1;if(calc(u,mid)<calc(s[x],mid)) swap(s[x],u);if(calc(u,l)<calc(s[x],l)) ins(ls(x),l,mid,u);if(calc(u,r)<calc(s[x],r)) ins(rs(x),mid+1,r,u);
}
int query(int x,int l,int r,int p){if(l==r) return calc(s[x],p);int res=calc(s[x],p);int mid=(l+r)>>1;if(p<=mid) res=min(res,query(ls(x),l,mid,p));else res=min(res,query(rs(x),mid+1,r,p));return res;
}
signed main(){// freopen("1.in","r",stdin);freopen("1.out","w",stdout);freopen("far.in","r",stdin);freopen("far.out","w",stdout);n=in(),m=in(),S=in(),T=in();for(int i=1;i<=m;++i){int u=in(),v=in(),w=in();add1(u,v,w);add2(v,u,w);}dij1(S);dij2(T);k[0]=0,b[0]=1e18;for(int i=1;i<=n;i++)k[i]=-2*i,b[i]=i*i+dis2[i];for(int i=1;i<=n;i++)ins(1,1,n,i);int ans=1e18;for(int i=1;i<=n;i++)ans=min(ans,query(1,1,n,i)+i*i+dis1[i]);printf("%lld\n",ans);return 0;
}
///////////////////////////////////////////////////
// ♪♪♪ //
///////////////////////////////////////////////////
//つ ◕_◕ つ
//༼ つ ◕_◕ ༽つ
注：此题卡spfa，所以请用dij