bzoj 4289: PA2012 Tax 最短路

        QAQ手写堆好慢。。。

       把一条边拆成两条，然后把边看成新的图中的点。然后考虑原图中的点x的入边和出边，显然每一条入边在原图中对应一条出边，在新图的这两个新点上连边为原图边权；然后由于一个点的点权为它出边和入边的较大值，考虑对x的出边排序后，相邻两条小的向大的连边权为原图边权差值，大的向小的连边权为0。

       然后处理一下源点和汇点即可。。

       不是很清楚SPFA能不能过，直接手写堆dijkstra。。

AC代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 400005
#define M 1600005
#define ll long long
using namespace std;

int n,m,cnt=1,tot,sta,gol,tp,pnt[M],edg[M],nxt[M],q[N],val[N],id[N]; ll d[N];
struct node{ int x,y,z; }a[N];
int read(){
	int x=0; char cr=getchar();
	while (cr<'0' || cr>'9') cr=getchar();
	while (cr>='0' && cr<='9'){ x=x*10+cr-'0'; cr=getchar(); }
	return x;
}
struct graph{
	int fst[N];
	void add(int x,int y,int z){
		pnt[++tot]=y; edg[tot]=z; nxt[tot]=fst[x]; fst[x]=tot;
	}
}g1,g2;
bool cmp(const node &u,const node &v){ return u.z<v.z; }
void godn(int x){
	int y=q[x],z=x<<1; ll t=d[y];
	if (z<cnt && d[q[z|1]]<d[q[z]]) z|=1;
	while (z<=cnt && d[q[z]]<t){
		q[x]=q[z]; id[q[x]]=x;
		x=z; z<<=1; if (z<cnt && d[q[z|1]]<d[q[z]]) z|=1;
	}
	q[x]=y; id[y]=x;
}
void goup(int x){
	int y=q[x]; ll t=d[y];
	for (; x>1 && t<d[q[x>>1]]; x>>=1){
		q[x]=q[x>>1]; id[q[x]]=x;
	}
	q[x]=y; id[y]=x;
}
int main(){
	n=read(); m=read(); int i,j,x,y,z,p;
	for (i=1; i<=m; i++){
		x=read(); y=read(); z=read();
		val[++cnt]=z; val[++cnt]=z;
		g1.add(x,cnt,cnt-1); g1.add(y,cnt-1,cnt);
	}
	for (i=2; i<n; i++){
		for (p=g1.fst[i],tp=0; p; p=nxt[p])
			a[++tp]=(node){pnt[p],edg[p],val[pnt[p]]};
		if (!tp) continue;
		sort(a+1,a+tp+1,cmp);
		for (j=1; j<=tp; j++) g2.add(a[j].x,a[j].y,a[j].z);
		for (j=1; j<tp; j++){
			g2.add(a[j].y,a[j+1].y,a[j+1].z-a[j].z);
			g2.add(a[j+1].y,a[j].y,0);
		}
	}
	gol=++cnt;
	for (p=g1.fst[1]; p; p=nxt[p])
		g2.add(1,edg[p],val[edg[p]]);
	for (p=g1.fst[n]; p; p=nxt[p])
		g2.add(pnt[p],gol,val[edg[p]]);
	for (i=1; i<=cnt; i++) q[i]=id[i]=i;
	memset(d,0x3f,sizeof(d)); d[1]=0;
	while (cnt){
		x=q[1]; q[1]=q[cnt--]; id[q[1]]=1; godn(1);
		for (p=g2.fst[x]; p; p=nxt[p]){
			y=pnt[p];
			if (d[x]+edg[p]<d[y]){
				d[y]=d[x]+edg[p]; goup(id[y]);
			}
		}
	}
	printf("%lld\n",d[gol]);
	return 0;
}

by lych

2016.5.31