bzoj 3551: [ONTAK2010]Peaks加强版
3551: [ONTAK2010]Peaks加强版
Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 3889 Solved: 1143
[Submit][Status][Discuss]
Description
【题目描述】同3545
Input
第一行三个数N，M，Q。
第二行N个数，第i个数为h_i
接下来M行，每行3个数a b c，表示从a到b有一条困难值为c的双向路径。
接下来Q行，每行三个数v x k，表示一组询问。v=v xor lastans，x=x xor lastans，k=k xor lastans。如果lastans=-1则不变。
Output
同3545
Sample Input
Sample Output
HINT
【数据范围】同3545
这题是kruskal重构树，这题一定要在线，而3545可以离线，离线的话我们可以排序以后用线段树合并维护连通性。
构树以后，由于答案都在子树里面，而子树在dfs序里面是一个连续的区间，所以我们维护一个主席树。
1 #include<bits/stdc++.h> 2 using namespace std; 3 int const N=200000+10; 4 int const M=500000+10; 5 #define mid (l+r)/2 6 struct edge{ 7 int to,nt; 8 }e[N]; 9 struct E{ 10 int x,y,z; 11 bool operator < (const E &rhs) const{ 12 return z<rhs.z; 13 } 14 }t[M]; 15 int a[N],b[N],n,m,q,c[N],f[N],val[N],h[N],cnt,sum[N*20],rt[N],tin[N],tout[N],T,id[N],lch[N*20],rch[N*20]; 16 int d[N][20],anc[N][20]; 17 void read(int &x){ 18 x=0; char c=getchar(); 19 while (!isdigit(c)) c=getchar(); 20 while (isdigit(c)) x=x*10+(c^48),c=getchar(); 21 } 22 int find(int x){ 23 return x==f[x]? x: f[x]=find(f[x]); 24 } 25 void add(int a,int b){ 26 e[++cnt].to=b; 27 e[cnt].nt=h[a]; 28 h[a]=cnt; 29 } 30 void kruscal(){ 31 int id=n; 32 for(int i=1;i<2*n;i++) f[i]=i; 33 for(int i=1;i<=m;i++){ 34 int fx=find(t[i].x); 35 int fy=find(t[i].y); 36 if(fx!=fy){ 37 id++; 38 f[fx]=f[fy]=id; 39 val[id]=t[i].z; 40 add(id,fx); 41 add(id,fy); 42 d[fx][0]=d[fy][0]=t[i].z; 43 anc[fx][0]=anc[fy][0]=id; 44 if(id==2*n-1) break; 45 } 46 } 47 } 48 void dfs(int x){ 49 tin[x]=++T; 50 id[T]=x; 51 for(int i=h[x];i;i=e[i].nt){ 52 int v=e[i].to; 53 dfs(v); 54 } 55 tout[x]=T; 56 } 57 void insert(int now,int &x,int l,int r,int p){ 58 x=++cnt; 59 lch[x]=lch[now];rch[x]=rch[now];sum[x]=sum[now]; 60 if(l==r) { 61 sum[x]++; 62 return; 63 } 64 if(p<=mid) insert(lch[now],lch[x],l,mid,p); 65 else insert(rch[now],rch[x],mid+1,r,p); 66 sum[x]=sum[lch[x]]+sum[rch[x]]; 67 } 68 int query(int x,int y,int l,int r,int k){ 69 if(l==r) return b[r]; 70 int num=sum[lch[y]]-sum[lch[x]]; 71 if(num>=k) return query(lch[x],lch[y],l,mid,k); 72 else return query(rch[x],rch[y],mid+1,r,k-num); 73 } 74 int main(){ 75 //scanf("%d%d%d",&n,&m,&q); 76 read(n);read(m); read(q); 77 for(int i=1;i<=n;i++) 78 read(a[i]),b[i]=a[i]; 79 sort(b+1,b+n+1); 80 int rk=unique(b+1,b+n+1)-b-1; 81 for(int i=1;i<=n;i++) 82 c[i]=lower_bound(b+1,b+rk+1,a[i])-b; 83 int last=0; 84 for(int i=1;i<=m;i++) 85 // scanf("%d%d%d",&t[i].x,&t[i].y,&t[i].z); 86 read(t[i].x),read(t[i].y),read(t[i].z); 87 sort(t+1,t+m+1); 88 kruscal(); 89 cnt=0; 90 for(int i=2*n-1;i>=1;i--) 91 if(!tin[i]) { 92 anc[i][0]=i; 93 dfs(i); 94 } 95 /* 96 for(int i=2*n-1;i>=1;i--){ 97 cout<<i<<"*** "<<val[i]<<" ** "; 98 for(int j=h[i];j;j=e[j].nt){ 99 cout<<e[j].to<<" "; 100 } 101 cout<<endl; 102 } */ 103 for(int j=1;j<20;j++) 104 for(int i=1;i<2*n;i++){ 105 d[i][j]=max(d[i][j-1],d[anc[i][j-1]][j-1]); 106 anc[i][j]=anc[anc[i][j-1]][j-1]; 107 } 108 for(int i=1;i<2*n;i++){ 109 if(id[i]<=n) insert(rt[i-1],rt[i],1,rk,c[id[i]]); 110 else rt[i]=rt[i-1]; 111 } 112 while (q--){ 113 int v,x,k; 114 // scanf("%d%d%d",&v,&x,&k); 115 read(v); read(x);read(k); 116 if(last!=-1){ 117 v^=last; 118 x^=last; 119 k^=last; 120 } 121 int t=v; 122 for(int i=19;i>=0;i--) 123 if(d[t][i]<=x) t=anc[t][i]; 124 int l=tin[t]; 125 int r=tout[t]; 126 int num=sum[rt[r]]-sum[rt[l-1]]; 127 if(num<k) printf("%d\n",last=-1); 128 else { 129 k=num-k+1; 130 //cout<<"nima"<<k<<endl; 131 printf("%d\n",last=query(rt[l-1],rt[r],1,rk,k)); 132 } 133 //last=0; 134 } 135 return 0; 136 137 }
