[差分建图 最短路] BZOJ4289： PA2012. Tax

大概套路就是把边看作点

边边之间连边的话边数会被卡成 $O(n^2)$

把边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;
}