大概套路就是把边看作点
边边之间连边的话边数会被卡成 O(n2)
把边sort一遍,可以用前缀和优化建图
但是这样常数有点大…反正我是卡完内存卡时间
然后发现可以差分…对于一个点,一条无向边拆成两条有向边,那么进去的边与对应的出去的边连0的边。把边sort一边后,这条边向它后面的边连一条边权大小为他们的差值的边,向它前面的边连边权大小为0的边。
然后跑最短路就好了
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#define fi first
#define se second
using namespace std;
const int N=4000010;
typedef long long ll;
int n,m,cnt,G[100010];
struct edge{
int s,t,w,nx;
}E[400010];
inline char nc(){
static char buf[100000],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline void read(int &x){
char c=nc(); x=0;
for(;c>'9'||c<'0';c=nc());for(;c>='0'&&c<='9';x=x*10+c-'0',c=nc());
}
inline void addedge(int x,int y,int z){
E[++cnt].s=x; E[cnt].t=y; E[cnt].nx=G[x]; E[cnt].w=z; G[x]=cnt;
E[++cnt].s=y; E[cnt].t=x; E[cnt].nx=G[y]; E[cnt].w=z; G[y]=cnt;
}
vector<int> e;
inline bool cmp(const int &a,const int &b){
return E[a].w<E[b].w;
}
int t;
typedef pair<ll,int> PAR;
namespace Gr{
int cnt,G[N];
struct edge{
int t,w,nx;
}E[N];
inline void addedge(int x,int y,int z){
E[++cnt].t=y; E[cnt].nx=G[x]; E[cnt].w=z; G[x]=cnt;
}
priority_queue<PAR> Q;
ll dis[N]; int vis[N];
ll spfa(){
for(int i=1;i<=t;i++) dis[i]=1LL<<60;
for(int i=::G[1];i;i=::E[i].nx){
dis[i]=::E[i].w; Q.push(PAR(-dis[i],i));
}
while(!Q.empty()){
int x=Q.top().se; ll y=-Q.top().fi; Q.pop();
if(y!=dis[x]) continue;
for(int i=G[x];i;i=E[i].nx)
if(dis[E[i].t]>dis[x]+E[i].w){
dis[E[i].t]=dis[x]+E[i].w;
Q.push(PAR(-dis[E[i].t],E[i].t));
}
}
ll ret=1LL<<60;
for(int i=::G[n];i;i=::E[i].nx)
ret=min(ret,::E[i].w+dis[i^1]);
return ret;
}
}
int main(){
read(n); read(m); cnt=1;
for(int i=1,x,y,z;i<=m;i++)
read(x),read(y),read(z),addedge(x,y,z);
t=cnt;
for(int x=1;x<=n;x++){
e.clear();
for(int i=G[x];i;i=E[i].nx)
e.push_back(i);
if(!e.size()) continue;
sort(e.begin(),e.end(),cmp);
for(int i=0;i<e.size();i++) Gr::addedge(e[i]^1,e[i],E[e[i]].w);
for(int i=0;i<e.size()-1;i++) Gr::addedge(e[i],e[i+1],E[e[i+1]].w-E[e[i]].w);
for(int i=1;i<e.size();i++) Gr::addedge(e[i],e[i-1],0);
}
printf("%lld\n",Gr::spfa());
//for(int i=2;i<=cnt;i++) cout<<E[i].s<<' '<<E[i].t<<' '<<Gr::dis[i]<<endl;
return 0;
}