bzoj:3040: 最短路(road)
Description
N个点,M条边的有向图,求点1到点N的最短路(保证存在)。
1<=N<=1000000,1<=M<=10000000
Input
第一行两个整数N、M,表示点数和边数。
第二行六个整数T、rxa、rxc、rya、ryc、rp。
前T条边采用如下方式生成:
1.初始化x=y=z=0。
2.重复以下过程T次:
x=(x*rxa+rxc)%rp;
y=(y*rya+ryc)%rp;
a=min(x%n+1,y%n+1);
b=max(y%n+1,y%n+1);
则有一条从a到b的,长度为1e8-100*a的有向边。
后M-T条边采用读入方式:
接下来M-T行每行三个整数x,y,z,表示一条从x到y长度为z的有向边。
1<=x,y<=N,0<z,rxa,rxc,rya,ryc,rp<2^31
Output
一个整数,表示1~N的最短路。
Sample Input
3 3
0 1 2 3 5 7
1 2 1
1 3 3
2 3 1
0 1 2 3 5 7
1 2 1
1 3 3
2 3 1
Sample Output
2
HINT
【注释】
请采用高效的堆来优化Dijkstra算法。
![](https://images2015.cnblogs.com/blog/841250/201511/841250-20151127214823718-1112031862.png)
纯粹迪杰……
orz kpm 2S:
![](https://images2015.cnblogs.com/blog/841250/201511/841250-20151127215008718-810909465.png)
#include<cstdio> #include<algorithm> using namespace std; struct na{ int y,ne; long long z; na(){ ne=0; } }; struct dui{ long long v; int n; }; const long long INF=1e17; int n,c,m,l[1000001],r[1000001],x,y,z,pp,num=0,p[1000001],nu; na b[10000001]; dui d[1000001]; char cc; long long read(){ long long a=0; cc=getchar(); while(cc<'0'||cc>'9') cc=getchar(); while(cc>='0'&&cc<='9') a=a*10+cc-48,cc=getchar(); return a; } void in(int x,int y,int z){ num++; if (l[x]==0) l[x]=num;else b[r[x]].ne=num; b[num].y=y;b[num].z=z; r[x]=num; } void sw(int a,int b){ swap(p[d[a].n],p[d[b].n]); swap(d[a],d[b]); } void utz(int x){ while(x>1&&d[x].v<d[x>>1].v) sw(x,x>>1),x>>=1; } void dtz(int j){ j=j*2; while(j<=nu){ if (j<nu&&d[j].v>d[j+1].v) j++; sw(j>>1,j); j*=2; } sw(j>>1,nu); p[d[nu].n]=0; nu--; } int main(){ long long rxa,rxc,rya,ryc,rp,t; n=read();m=read(); t=read();rxa=read();rxc=read();rya=read();ryc=read();rp=read(); nu=n; x=y=z=0; for (int i=1;i<=t;i++){ long long a,b; x=(x*rxa+rxc)%rp; y=(y*rya+ryc)%rp; a=min(x%n+1,y%n+1); b=max(y%n+1,y%n+1); in(a,b,1e8-100*a); } for (int i=1;i<=m-t;i++){ x=read();y=read();z=read(); in(x,y,z); } for (int i=1;i<=n;i++) d[i].v=INF,d[i].n=i,p[i]=i; d[1].v=0; for (int i=1;i<=n;i++){ x=d[1].n;long long y=d[1].v; if (x==n) { printf("%lld\n",y); return 0; } dtz(1); for (int j=l[x];j;j=b[j].ne) if (p[b[j].y]&&d[p[b[j].y]].v>y+b[j].z){ d[p[b[j].y]].v=y+b[j].z; utz(p[b[j].y]); } } }