codeforces464E The Classic Problem -- 最短路+主席树+Hash

题目大意：

求$s$到$t$的最短路，边权为$2^{x_i},x_i\leq10^5$。

看Manchery的题解就可以了。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<cmath>
#include<queue>
#include<map>
using namespace std;
#define N 100010
#define P 1000000007
#define INF 1e6
typedef unsigned long long ull;
struct Edge{
    int t,nx,w;
}e[N<<1];
struct Node{
    ull w;
    int l,r,f;
}c[N*280];
struct Interval{
    int l,r;
    Interval(int l=0,int r=0):l(l),r(r){}
    bool operator < (Interval a)const{
        return l<a.l||(l==a.l&&r<a.r);
    }
};
map<Interval,int>Map;
ull p[N+100];
int i,j,k,n,m,h[N],x,y,z,d[N],Num,t,M,Pre[N],Cnt,g[N],B;
bool b[N];
inline bool Cmp(int x,int y,int l,int r){
    if(l==r)return c[x].w>c[y].w;
    int Mid=l+r>>1;
    if(c[c[x].r].w==c[c[y].r].w)return Cmp(c[x].l,c[y].l,l,Mid);
    return Cmp(c[x].r,c[y].r,Mid+1,r);
}
struct A{
    int Rt,f;
    A(int Rt=0,int f=0):Rt(Rt),f(f){}
    bool operator < (A a)const{
        return Cmp(Rt,a.Rt,1,M);
    }
}Tmp;
priority_queue<A>Q;
inline int Max(int x,int y){
    return x<y?y:x;
}
inline int Min(int x,int y){
    return x<y?x:y;
}
inline void Add(int x,int y,int z){
    e[++t].t=y;e[t].nx=h[x];h[x]=t;e[t].w=z;
}
inline void Up(int x){
    c[x].w=c[c[x].l].w*7+c[c[x].r].w*11;
    c[x].f=c[c[x].l].f<INF?c[c[x].l].f:c[c[x].r].f;
}
inline void Build(int& x,int l,int r){
    x=++Num;Map[Interval(l,r)]=x;
    if(l==r){
        c[x].f=l;
        return;
    }
    int Mid=l+r>>1;
    Build(c[x].l,l,Mid);
    Build(c[x].r,Mid+1,r);
    Up(x);
}
inline int Update2(int x,int l,int r,int y){
    int z=++Num;c[z]=c[x];
    if(l==r){
        c[z].w=p[l];
        c[z].f=INF;
        return z;
    }
    int Mid=l+r>>1;
    if(y<=Mid)c[z].l=Update2(c[x].l,l,Mid,y);else c[z].r=Update2(c[x].r,Mid+1,r,y);
    Up(z);
    return z;
}
inline int Update3(int x,int l,int r,int L,int R){
    if(l>=L&&r<=R)return Map[Interval(l,r)];
    int Mid=l+r>>1,y=++Num;c[y]=c[x];
    if(Mid>=L)c[y].l=Update3(c[x].l,l,Mid,L,R);
    if(Mid<R)c[y].r=Update3(c[x].r,Mid+1,r,L,R);
    Up(y);
    return y;
}
inline int Find(int x,int l,int r,int y){
    if(l==r)return c[x].f;
    int Mid=l+r>>1;
    if(y<=Mid)return Min(Find(c[x].l,l,Mid,y),c[c[x].r].f);
    return Find(c[x].r,Mid+1,r,y);
}
inline int Update1(int x,int y){
    int z=Find(x,1,M,y);
    int a1=Update2(x,1,M,z);
    if(z>y)a1=Update3(a1,1,M,y,z-1);
    return a1;
}
inline int Getval(int x,int l,int r){
    if(l==r){
        int ret=c[x].w?B:0;
        return B=(B<<1)%P,ret;
    }
    int Mid=l+r>>1;
    return (Getval(c[x].l,l,Mid)+Getval(c[x].r,Mid+1,r))%P;
}
int main(){
    scanf("%d%d",&n,&m);
    for(i=1;i<=m;i++)scanf("%d%d%d",&x,&y,&z),Add(x,y,++z),Add(y,x,z),M=Max(M,z);
    for(i=1;1<<i<n;i++);M+=i;
    for(i=p[0]=1;i<=M;i++)p[i]=p[i-1]*5;
    scanf("%d%d",&x,&y);
    if(x==y){
        printf("0\n1\n%d\n",x);
        return 0;
    }
    Build(d[x],1,M);
    Q.push(A(d[x],x));
    while(!Q.empty()){
        Tmp=Q.top();Q.pop();
        if(d[Tmp.f]!=Tmp.Rt)continue;
        for(int i=h[Tmp.f];i;i=e[i].nx){
            z=Update1(d[Tmp.f],e[i].w);
            if(!b[e[i].t]||Cmp(d[e[i].t],z,1,M)){
                b[e[i].t]=1;d[e[i].t]=z;Pre[e[i].t]=Tmp.f;
                Q.push(A(z,e[i].t));
            }
        }
    }
    if(!b[y])printf("-1\n");else{
        B=1;
        printf("%d\n",Getval(d[y],1,M));
        for(i=y;i!=x;i=Pre[i])g[++Cnt]=i;
        printf("%d\n%d ",Cnt+1,x);
        for(i=Cnt;i;i--)printf("%d ",g[i]);
    }
    return 0;
}