bzoj3551

3551: [ONTAK2010]Peaks加强版

Time Limit: 20 Sec  Memory Limit: 128 MB
Submit: 877  Solved: 297
[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

Source

题解：

　　过了这么久才又开始写博客（省选被虐成渣，现在还要准备半期考试 only shit）好吧，算了毕竟太渣，不多说了，写题解QAQ

　　看了题目，一脸懵逼，怎么写？（不会） 怎么办？（看题解）

　　首先看看这道题，它是一个无向图，根本不知道怎么来写主席树，操蛋。

　　仔细看看题，我们可以这样想，我要走一个边权比x小的边，并且一直走下去，我们不就可以用最小生成树吗？每次走最小的边，那么一定是最优的，所以我们就用kruskal来写

　　但是要注意的是，在找到一个新的边时，新建一个节点，权值为边的权值，把边的两端连向这个点，不难证明这是一课树。。。。。

　　然后每次询问就从当前点往上面倍增地跳，直到不能走（点权大于x） 那么那个点的子树中的点就是我可以到达的点

　　于是题目就变成求一个子树第k大的问题。。。用lca加和主席树就A了

 #include<iostream>
 #include<cstring>
 #include<cstdio>
 #include<algorithm>
 #include<cmath>
 #define maxn 300005
 #define maxnode 2000005
 using namespace std;
 int n,m,q,lastans,top;
 int size,tot,ext;
 int pre[maxn],now[maxn],v[maxn],fa[maxn],sum[maxnode],son[maxnode][],pp[maxn];
 int bin[],deep[maxn],f[maxn][],mx[maxn][];
 int h[maxn],list[maxn],qz[maxn],root[maxn],st[maxn],ed[maxn];
 bool vis[maxn];
 struct date{int u,v,val;
 }a[];
 int read()
 {
     int x=; char ch; bool bo=;
     while (ch=getchar(),ch<''||ch>'') if (bo=='-') bo=;
     while (x=x*+ch-'',ch=getchar(),ch>=''&&ch<='');
     if (bo) return -x; return x;
 }
 bool cmp_val(date a, date b)
 {
     return a.val<b.val;
 }
 int find(int x)
 {
     if (fa[x]!=x) fa[x]=find(fa[x]);
     return fa[x];
 }
 void insert(int x,int y){tot++; pre[tot]=now[x]; now[x]=tot; v[tot]=y;
 }
 void ins(int l,int r,int x,int &y,int val)
 {
     int t;
     if (!y) y=++tot; sum[y]=sum[x]+;
     if (l==r) return;
     int mid=(l+r)>>;
     if (val<=mid) t=,r=mid; else t=,l=mid+;
     son[y][t^]=son[x][t^];
     ins(l,r,son[x][t],son[y][t],val);
 }
 void dfs(int x)
 {
     vis[x]=; pp[++top]=x;
     for (int i=; i<=; i++)
         if (bin[i]<=deep[x])
         {
             f[x][i]=f[f[x][i-]][i-];
             mx[x][i]=max(mx[x][i-],mx[f[x][i-]][i-]);
         }
     for (int p=now[x]; p; p=pre[p])
     {
         int son=v[p];
         deep[son]=deep[x]+;
         f[son][]=x;
         mx[son][]=qz[x];
         dfs(son);
     }
     if (x>n) pp[++top]=x;
 }
 void build()
 {
     ext=n;
     sort(a+,a+m+,cmp_val);
     for (int i=; i<=m; i++)
     {
         int q=find(a[i].u),p=find(a[i].v);
         if (q!=p)
         {
             ext++;
             fa[q]=fa[p]=ext; qz[ext]=a[i].val;
             insert(ext,p);  insert(ext,q);
             if (ext==*n-) break;
         }
     }
     for (int i=; i<=n; i++)
     {
         if (!vis[i]) dfs(find(i));
     }
     for (int i=; i<=top; i++)
     {
         int t=pp[i];
         if (t<=n) ins(,n,root[i-],root[i],h[t]);
         else
         {
             root[i]=root[i-];
             if (!st[t]) st[t]=i; else ed[t]=i;
         }
     }
 }
 int valfind(int x,int val)
 {
     for(int i=;i>=;i--)
         if(deep[x]>=bin[i]&&mx[x][i]<=val)x=f[x][i];
     return x;
 }
 int query(int l,int r,int x,int y,int val)
 {
     if (l==r) return l;
     int mid=(l+r)>>;
     int kk=sum[son[y][]]-sum[son[x][]],t;
     //cout<<"         "<<kk<<endl;
     if (kk>=val) return query(l,mid,son[x][],son[y][],val);
     else return query(mid+,r,son[x][],son[y][],val-kk);
 }
 void solve()
 {
     while (q--)
     {
         int x=read(),val=read(),k=read();
         if (lastans!=-) x^=lastans,val^=lastans,k^=lastans;
         int t=valfind(x,val);
         int a=root[st[t]],b=root[ed[t]];
         if (sum[b]-sum[a]<k) lastans=-;
         else lastans=list[query(,n,a,b,sum[b]-sum[a]-k+)];
         printf("%d\n",lastans);
     }

 }
 void prework()
 {
     bin[]=;
     for (int i=; i<; i++) bin[i]=bin[i-]<<;
     for (int i=; i<=*n; i++) fa[i]=i;
     sort(list+,list++n);
     for (int i=; i<=n; i++) h[i]=lower_bound(list+,list++n,h[i])-list;
     for (int i=; i<=m; i++)
     {
         a[i].u=read(); a[i].v=read(); a[i].val=read();
     }
     build();
 }
 int main()
 {
     n=read(); m=read(); q=read();
     for (int i=; i<=n; i++) list[i]=h[i]=read();
     prework();
     solve();
     return ;
 }