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 ;
}