BZOJ4289 : PA2012 Tax
一个直观的想法是把每条边拆成两条有向边,同时每条有向边是新图中的一个点。对于两条边a->b与b->c,两点之间连有向边,费用为两条边费用的最大值。然后新建源点S与汇点T,由S向所有起点为1的边连边,T接受所有终点为n的边,那么答案就是S到T的最短路。
这样子的边数为$O(m^2)$,不能承受。
考虑枚举中转点x,将所有与它有关的边按费用从小到大排序。对于每条边,从以x为终点的点向以x为起点的点连边,费用为该边的费用。从以x为起点的点向下一条边连边,费用为两条边费用的差值,向上一条边连边,费用为0。
这样子建图,边数为$O(m)$,可以承受。
#include<cstdio> #include<algorithm> #include<queue> using namespace std; typedef long long ll; typedef pair<int,ll> PI; const int N=100010,M=200010; int n,m,i,j,x,y,z,cnt,t,S,T; int g[N],en[M<<1],st[M<<1],nxt[M<<1],ed; int G[M<<1],V[M*6],W[M*6],NXT[M*6],ED; ll d[M<<1]; priority_queue<PI,vector<PI>,greater<PI> >Q; struct P{int x,y,z;P(){}P(int _x,int _y,int _z){x=_x,y=_y,z=_z;}}a[M<<1],q[M]; inline bool cmp(const P&a,const P&b){return a.z<b.z;} inline void read(int&a){char c;while(!(((c=getchar())>='0')&&(c<='9')));a=c-'0';while(((c=getchar())>='0')&&(c<='9'))(a*=10)+=c-'0';} inline void addedge(int x,int y,int z){en[++ed]=y;st[ed]=z;nxt[ed]=g[x];g[x]=ed;} inline void add(int x,int y,int z){V[++ED]=y;W[ED]=z;NXT[ED]=G[x];G[x]=ED;} int main(){ read(n),read(m); while(m--){ read(x),read(y),read(z); a[++cnt]=P(x,y,z); a[++cnt]=P(y,x,z); addedge(x,cnt,cnt-1); addedge(y,cnt-1,cnt); } for(i=1;i<=n;i++){ for(t=0,j=g[i];j;j=nxt[j])q[++t]=P(en[j],st[j],a[en[j]].z); if(!t)continue; sort(q+1,q+t+1,cmp); for(j=1;j<=t;j++)add(q[j].x,q[j].y,q[j].z); for(j=1;j<t;j++){ add(q[j].y,q[j+1].y,q[j+1].z-q[j].z); add(q[j+1].y,q[j].y,0); } } S=cnt+1;T=S+1; for(i=1;i<=cnt;i++){ if(a[i].x==1)add(S,i,a[i].z); if(a[i].y==n)add(i,T,a[i].z); } for(i=1;i<=T;i++)d[i]=1LL<<60; Q.push(PI(d[S]=0,S)); while(!Q.empty()){ PI t=Q.top();Q.pop(); if(d[x=t.second]<t.first)continue; for(i=G[x];i;i=NXT[i])if(d[x]+W[i]<d[V[i]])Q.push(PI(d[V[i]]=d[x]+W[i],V[i])); } return printf("%lld",d[T]),0; }