[Luogu4899][IOI2018] werewolf 狼人

luogu

sol

\(\mbox{IOI2018}\)的出题人有没有看过\(\mbox{NOI2018}\)的题目呀。。。

\(\mbox{Kruskal}\)重构树+二维数点。

题目相当于是问你从\(S\)出发只经过\(\ge L\)的点,和从\(T\)出发只经过\(\le R\)的点有没有交集。

分别建两棵最大/最小\(\mbox{Kruskal}\)重构树,这样从\(S\)出发,从\(T\)出发能到达的点就分别是两棵树上的一个子树。以点在两棵树上的\(dfs\)序作为两维坐标,每次询问就是问你某个矩形区域里面有没有点。

code

#include<cstdio>
#include<algorithm>
#include<vector>
using namespace std;
int gi(){
int x=0,w=1;char ch=getchar();
while ((ch<'0'||ch>'9')&&ch!='-') ch=getchar();
if (ch=='-') w=0,ch=getchar();
while (ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar();
return w?x:-x;
}
const int N = 4e5+5;
struct president_tree{int ls,rs,sz;}t[N*35];
int n,m,q,val[N],rt[N],tot;
vector<int>E[N];
struct Kruskal_rebuild_tree{
int ty;
int fa[N],nxt[N],hd[N],pa[20][N],dfn[N],low[N],tim;
int find(int x){return x==fa[x]?x:fa[x]=find(fa[x]);}
void dfs(int u){
dfn[u]=++tim;
for (int j=1;j<20;++j) pa[j][u]=pa[j-1][pa[j-1][u]];
for (int v=hd[u];v;v=nxt[v]) dfs(v);
low[u]=tim;
}
void work(){
for (int i=1;i<=n;++i) fa[i]=i;
if (ty){
for (int u=n;u;--u)
for (int v:E[u])
if (u<v){
v=find(v);if (u==v) continue;
nxt[v]=hd[u],hd[u]=v,fa[v]=pa[0][v]=u;
}
dfs(1);
}else{
for (int u=1;u<=n;++u)
for (int v:E[u])
if (u>v){
v=find(v);if (u==v) continue;
nxt[v]=hd[u],hd[u]=v,fa[v]=pa[0][v]=u;
}
dfs(n);
}
}
int jump(int u,int k){
for (int j=19;~j;--j)
if (pa[j][u]&&((ty&&pa[j][u]>=k)||(!ty&&pa[j][u]<=k)))
u=pa[j][u];
return u;
}
}A,B;
void modify(int &x,int l,int r,int p){
t[++tot]=t[x];++t[x=tot].sz;
if (l==r) return;int mid=l+r>>1;
if (p<=mid) modify(t[x].ls,l,mid,p);
else modify(t[x].rs,mid+1,r,p);
}
int query(int x,int y,int l,int r,int ql,int qr){
if (l>=ql&&r<=qr) return t[x].sz-t[y].sz;
int mid=l+r>>1,s=0;
if (ql<=mid) s+=query(t[x].ls,t[y].ls,l,mid,ql,qr);
if (qr>mid) s+=query(t[x].rs,t[y].rs,mid+1,r,ql,qr);
return s;
}
int main(){
n=gi();m=gi();q=gi();
for (int i=1;i<=m;++i){
int u=gi()+1,v=gi()+1;
E[u].push_back(v);E[v].push_back(u);
}
A.ty=1;A.work();B.work();
for (int i=1;i<=n;++i) val[A.dfn[i]]=B.dfn[i];
for (int i=1;i<=n;++i) modify(rt[i]=rt[i-1],1,n,val[i]);
while (q--){
int s=gi()+1,e=gi()+1,l=gi()+1,r=gi()+1;
s=A.jump(s,l);e=B.jump(e,r);
puts(query(rt[A.low[s]],rt[A.dfn[s]-1],1,n,B.dfn[e],B.low[e])?"1":"0");
}
return 0;
}
上一篇:设置tableview 右侧 索引字体的大小


下一篇:对C# .Net4.5异步机制测试