[最短路 主席树 Hash] Codeforces 464E #265 (Div. 1) E. The Classic Problem

21 篇文章 0 订阅
20 篇文章 0 订阅

很棒的数据结构练习题
因为边权太大了 我们不能直接储存
我们用主席树!
比较直接从高位到低位找第一个不同的位 可以Hash下
每次加法 我们找到这一位之后第一个0 那么0变为1 0之前一连串1都变成0 这里有个小trick可以不用打标记 我们先建一棵全0的线段树 清0直接找到对应节点接上去就好了

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<queue>
#include<map>
using namespace std;
typedef unsigned int uint;
typedef unsigned long long ull;
typedef pair<int,int> abcd;

inline char nc(){
  static char buf[100000],*p1=buf,*p2=buf;
  return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}

inline void read(int &x){
  char c=nc(),b=1;
  for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
  for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}

const int N=100205;
const int M=20000005;
const int P1=3,P2=5;

uint p1[N+5],p2[N+5];
int ls[M],rs[M],f[M]; ull H[M];
int rtw,ncnt;

inline void upd(int x,int l,int r){
  int mid=(l+r)>>1;
  f[x]=!f[ls[x]]?f[rs[x]]:f[ls[x]];
  uint x1=H[ls[x]]&((1LL<<32)-1),y1=H[ls[x]]>>32;
  uint x2=H[rs[x]]&((1LL<<32)-1),y2=H[rs[x]]>>32;
  uint xx=x1+x2*p1[mid-l+1],yy=y1+y2*p2[mid-l+1];
  H[x]=xx+((1LL<<32)*yy);
}

map<abcd,int> Pos;

inline void Build(int &x,int l,int r){
  x=++ncnt; int mid=(l+r)>>1;
  Pos[abcd(l,r)]=x;
  if (l==r){
    H[x]=0; f[x]=l; return;
  }
  Build(ls[x],l,mid); Build(rs[x],mid+1,r);
  upd(x,l,r);
}

inline void Pre(){
  p1[0]=p2[0]=1;
  for (int i=1;i<=N;i++) p1[i]=p1[i-1]*P1,p2[i]=p2[i-1]*P2;
  Build(rtw,1,N);
}

inline void set0(int &x,int y,int l,int r,int ql,int qr){
  if (ql<=l && r<=qr){
    x=Pos[abcd(l,r)];
    return;
  }
  x=++ncnt;
  ls[x]=ls[y],rs[x]=rs[y];
  int mid=(l+r)>>1;
  if (ql<=mid)
    set0(ls[x],ls[y],l,mid,ql,qr);
  if (qr>mid)
    set0(rs[x],rs[y],mid+1,r,ql,qr);
  upd(x,l,r);
}

inline void set(int &x,int y,int l,int r,int t){
  x=++ncnt;
  ls[x]=ls[y],rs[x]=rs[y];
  if (l==r){
    f[x]=0; H[x]=(1LL<<32)+1;
    return;
  }
  int mid=(l+r)>>1;
  if (t<=mid) set(ls[x],ls[y],l,mid,t);
  else set(rs[x],rs[y],mid+1,r,t);
  upd(x,l,r);
}

inline int query(int x,int l,int r,int ql,int qr){
  if (ql<=l && r<=qr)
    return f[x];
  int mid=(l+r)>>1;
  if (qr<=mid) return query(ls[x],l,mid,ql,qr);
  else if (ql>mid) return query(rs[x],mid+1,r,ql,qr);
  else{
    int d=query(ls[x],l,mid,ql,qr);
    if (d) return d; else return query(rs[x],mid+1,r,ql,qr);
  }
}

inline int Add(int x,int w){
  int t=query(x,1,N,w,N),ret,tmp;
  if (t>w)
    set0(tmp,x,1,N,w,t-1);
  else
    tmp=x;
  set(ret,tmp,1,N,t);
  return ret;
}

const int P=1000000007;
int ans=0,cur=1;
inline void getval(int x,int l,int r){
  if (l==r){
    ans=(ans+cur*(H[x]!=0))%P,cur=(cur<<1)%P;
    return;
  }
  int mid=(l+r)>>1;
  getval(ls[x],l,mid); getval(rs[x],mid+1,r);
}

inline bool cmp(int x,int y,int l,int r){
  if (l==r)
    return H[x]<H[y];
  int mid=(l+r)>>1;
  if (H[rs[x]]==H[rs[y]])
    return cmp(ls[x],ls[y],l,mid);
  else
    return cmp(rs[x],rs[y],mid+1,r);
}

struct info{
  int x,d;
  info(int x=0,int d=0):x(x),d(d) { }
  bool operator < (const info &B) const{
    return cmp(B.d,d,1,N);
  }
};

struct edge{
  int u,v,w,next;
}G[N<<1];
int head[N],inum;
inline void add(int u,int v,int w,int p){
  G[p].u=u; G[p].v=v; G[p].w=w; G[p].next=head[u]; head[u]=p;
}

int S,T;

priority_queue<info> Q;

int dis[N],pre[N];
int vst[N];
#define V G[p].v
inline void Dij(){
  Q.push(info(S,rtw)); vst[S]=1;
  dis[S]=rtw;
  while (!Q.empty()){
    int u=Q.top().x,d=Q.top().d; Q.pop();
    if (dis[u]!=d) continue;
    for (int p=head[u];p;p=G[p].next){
      int nw=Add(d,G[p].w);
      if (!vst[V] || cmp(nw,dis[V],1,N)){
    vst[V]=1; dis[V]=nw; pre[V]=u;
    Q.push(info(V,nw));
      }
    }
  }
}

int n,m;
int pnt,lst[N];

int main(){
  int iu,iv,iw;
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  read(n); read(m); Pre();
  for (int i=1;i<=m;i++)
    read(iu),read(iv),read(iw),iw++,add(iu,iv,iw,++inum),add(iv,iu,iw,++inum);
  read(S); read(T);
  Dij();
  if (!vst[T])
    printf("-1\n");
  else{
    getval(dis[T],1,N);
    printf("%d\n",ans);
    for (int t=T;t;t=pre[t]) lst[++pnt]=t;
    printf("%d\n",pnt);
    for (int i=pnt;i;i--) printf("%d ",lst[i]);
  }
  return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值