标签：IOI2018 int Kruskal Enum 狼人 DFS son ++ define
这题不就是Peaks(加强版)或者归程么。。这算是\(IOI2018\)撞上\(NOI2018\)的题了？
\(Kruskal\)重构树（具体是所有点按从小到大/从大到小的顺序，依次加入这些点的边），我们可以得到两棵树（和那两题不一样的是这题的权值在点上，不需要新建节点）。
对于询问\((S,T,L,R)\)，可以倍增找出\(S,T\)可以在哪棵子树中随便走。
那么只需要判断两棵子树是否有交就可以惹。
注意到子树的DFS序是连续的，我们可以在第一个子树的某个数据结构上，查第二个子树的DFS序区间。用主席树就可以了，主席树的下标是每个点在第二棵树上的DFS序。判断第一个子树在第二个子树的DFS序区间中是否有值即可。
写一写权当复习惹。
另外，这个题的交互有点逗逼啊...为什么是所有询问一起输出啊？
LOJ上15分，T掉了，不知道为什么。不管惹。
//656ms 84.03MB
#include <cstdio>
#include <cctype>
#include <algorithm>
#define BIT 17
#define gc() getchar()
#define MAXIN 300000
//#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
typedef long long LL;
const int N=2e5+5,M=8e5+5;
int n,root[N];
char IN[MAXIN],*SS=IN,*TT=IN;
struct Graph
{
int Enum,H[N],nxt[M],to[M];
inline void AE(int u,int v)
{
to[++Enum]=v, nxt[Enum]=H[u], H[u]=Enum;
to[++Enum]=u, nxt[Enum]=H[v], H[v]=Enum;
}
}G;
struct Tree
{
int fg,fa[N][BIT+1],F[N],H[N],nxt[N],L[N],R[N],Index;
int Find(int x)
{
return x==F[x]?x:F[x]=Find(F[x]);
}
inline void AE(int u,int v)
{
nxt[v]=H[u], H[u]=v;
}
inline int Jump(int x,int k)
{
if(fg)
{
for(int i=BIT; ~i; --i)
if(fa[x][i]&&fa[x][i]<=k) x=fa[x][i];
}
else
{
for(int i=BIT; ~i; --i)
if(fa[x][i]>=k) x=fa[x][i];
}
return x;
}
void DFS(int x,int dep)
{
L[x]=++Index;
for(int i=1; 1<<i<=dep; ++i) fa[x][i]=fa[fa[x][i-1]][i-1];
for(int v=H[x]; v; v=nxt[v]) fa[v][0]=x, DFS(v,dep+1);
R[x]=Index;
}
void Build(const int n)
{
for(int i=1; i<=n; ++i) F[i]=i;
if(fg)
{
for(int x=1,k=1; x<=n; ++x)
for(int i=G.H[x],v; i; i=G.nxt[i])
if((v=G.to[i])<x)
{
int r1=Find(x),r2=Find(v);//r1==x
if(r1!=r2)
{
F[r2]=r1, AE(r1,r2);
if(++k==n) break;
}
}
}
else
{
for(int x=n,k=1; x; --x)
for(int i=G.H[x],v; i; i=G.nxt[i])
if((v=G.to[i])>x)
{
int r1=Find(x),r2=Find(v);
if(r1!=r2)
{
F[r2]=r1, AE(r1,r2);
if(++k==n) break;
}
}
}
DFS(fg?n:1,0);
}
}T1,T2;
struct Segment_Tree
{
#define ls son[x][0]
#define rs son[x][1]
#define lson ls,son[y][0],l,m
#define rson rs,son[y][1],m+1,r
#define S N*19
int tot,sz[S],son[S][2];
#undef S
void Modify(int &x,int y,int l,int r,int p)
{
sz[x=++tot]=sz[y]+1;
if(l==r) return;
int m=l+r>>1;
p<=m?(rs=son[y][1],Modify(lson,p)):(ls=son[y][0],Modify(rson,p));
}
int Query(int x,int y,int l,int r,int L,int R)//y-x
{
if(sz[y]-sz[x]<=0) return 0;
if(L<=l && r<=R) return 1;
int m=l+r>>1;
if(L<=m)
if(m<R) return Query(lson,L,R)||Query(rson,L,R);
else return Query(lson,L,R);
return Query(rson,L,R);
}
}T;
inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-48,c=gc());
return now;
}
void DFS(int x)
{
static int Index=0;
++Index, T.Modify(root[Index],root[Index-1],1,n,T2.L[x]);
for(int v=T1.H[x]; v; v=T1.nxt[v]) DFS(v);
}
int main()
{
int n=read(),m=read(),Q=read(); ::n=n;
for(int i=1; i<=m; ++i) G.AE(read()+1,read()+1);
T1.Build(n), T2.fg=1, T2.Build(n), DFS(1);
for(int s,t; Q--; )
s=read()+1,t=read()+1,s=T1.Jump(s,read()+1),t=T2.Jump(t,read()+1),
putchar(T.Query(root[T1.L[s]-1],root[T1.R[s]],1,n,T2.L[t],T2.R[t])?'1':'0'), putchar('\n');
return 0;
}标签：IOI2018,int,Kruskal,Enum,狼人,DFS,son,++,define 来源： https://www.cnblogs.com/SovietPower/p/10685264.html
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。