标签:
最短路+堆优化。
普通的堆还不行,自己用的是配对堆(貌似斐波那契堆也行?毕竟理论复杂度)
然后发现自己的配对堆比云神的不知快了多少。。。我照着他的模版打的喂。。
然后发现前T条边不理都能A。。。
数据啊数据233
#include <cctype> #include <algorithm> #include <cstring> #include <cstdlib> #include <cstdio> #include <iostream> #define rep(i, l, r) for(int i=l; i<=r; i++) #define clr(x, c) memset(x, c, sizeof(x)) #define travel(x) for(edge *p=fir[x]; p; p=p->n) #define ll long long #define maxn 1000009 #define maxm 10000009 using namespace std; const ll inf=(ll)(0x7fffffff)*(ll)(0x7fffffff); inline ll read() { ll x=0; char ch=getchar(); while (!isdigit(ch)) ch=getchar(); while (isdigit(ch)) x=x*10+ch-‘0‘, ch=getchar(); return x; } struct edge{int y, z; edge *n;} e[maxm], *fir[maxn], *pt=e; inline void Add(ll x, ll y, ll z){pt->y=y, pt->z=z, pt->n=fir[x], fir[x]=pt++;} struct node{int l, r, ch; node(){l=r=ch;}} h[maxn]; int n, m; ll d[maxn]; bool f[maxn]; int roof=0; int join(int v, int u) { if (d[v]<d[u]) swap(v, u); h[v].l=u, h[v].r=h[u].ch, h[h[u].ch].l=v; h[u].ch=v; return u; } void push(int v){if (!roof) roof=v; else roof=join(roof, v);} void update(int v) { if (v!=roof) { if (h[h[v].l].ch==v) h[h[v].l].ch=h[v].r; else h[h[v].l].r=h[v].r; if (h[v].r) h[h[v].r].l=h[v].l; h[v].l=h[v].r=0; roof=join(roof, v); } } int st[maxn], top; void pop() { if (!h[roof].ch) roof=0; else { top=0; int t=h[roof].ch; while (t) if (h[t].r) { int k=h[h[t].r].r; int v=h[t].r; h[t].l=h[t].r=h[v].l=h[v].r=0; st[++top]=join(t, v); t=k; } else { st[++top]=t; h[t].l=h[t].r=0; break; } roof=st[top]; rep(i, 1, top-1) roof=join(roof, st[i]); } } void Dijstra() { clr(f, 0); d[1]=0; rep(i, 2, n) d[i]=inf; push(1); f[1]=true; rep(i, 1, n) { int x=roof; pop(); f[x]=0; if (x==n) break; travel(x) if (d[p->y]>d[x]+(ll)(p->z)) { d[p->y]=d[x]+(ll)(p->z); if (!f[p->y]) push(p->y), f[p->y]=1; else update(p->y); } } } int main() { n=read(), m=read(); ll t=read(), rxa=read(), rxc=read(), rya=read(), ryc=read(), rp=read(), x, y, a, b; /* rep(i, 1, 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); Add(a, b, 1e8-100*a); } */ rep(i, 1, m-t){x=read(), y=read(); Add(x, y, read());} Dijstra(); printf("%lld\n", d[n]); return 0; }
标签:
原文地址:http://www.cnblogs.com/NanoApe/p/4479245.html