http://acm.hdu.edu.cn/showproblem.php?pid=6203
给出q条链 由u和v确定 每条链上至少有一个节点是有故障的 问整个树图中至少有多少故障节点
对于每条链 求出lca 再按lca的深度降序排序
然后对于一条链 如果uv节点都没有被其他链覆盖过 那就将lca对应的整棵子树标记覆盖 答案加加 否则就说明当前这条链上的故障点已经可以和别的链合并了 可以忽略
至于为什么要按深度降序排序 我认为这样每次只需要判断一条链是不是已经通过其他链确定 如果升序排序 每一次要看lca到u和v这条链上有多少其它链上的lca被影响 很难写 (借鉴)同时 由于优先处理 LCA 深度大的点 不会出现点 U V 同时在同一个被禁止通行点 P 的子树内
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
using namespace std;
struct node0
{
int u;
int v;
int lca;
};
struct node1
{
int v;
int next;
};
node0 pre[50010];
node1 edge[60010];
int dp[30010][15];
int val[120010];
int first[30010],deep[30010],mp[30010],sum[30010];
int n,q,num;
bool cmp(node0 n1,node0 n2)
{
return deep[n1.lca]>deep[n2.lca];
}
void addedge(int u,int v)
{
edge[num].v=v;
edge[num].next=first[u];
first[u]=num++;
}
void dfs(int cur,int fa)
{
int i,v;
mp[cur]=++num,sum[cur]=1;
for(i=first[cur];i!=-1;i=edge[i].next)
{
v=edge[i].v;
if(v!=fa)
{
dp[v][0]=cur;
deep[v]=deep[cur]+1;
dfs(v,cur);
sum[cur]+=sum[v];
}
}
return;
}
void solve()
{
int i,j;
dp[1][0]=0;
deep[1]=1;
num=0;
dfs(1,0);
for(j=1;(1<<j)<=n;j++)
{
for(i=1;i<=n;i++)
{
dp[i][j]=dp[dp[i][j-1]][j-1];
}
}
return;
}
int getlca(int u,int v)
{
int i;
if(deep[u]<deep[v]) swap(u,v);
for(i=log2(n);i>=0;i--)
{
if(deep[dp[u][i]]>=deep[v])
{
u=dp[u][i];
}
}
if(u==v) return u;
for(i=log2(n);i>=0;i--)
{
if(dp[u][i]!=dp[v][i])
{
u=dp[u][i];
v=dp[v][i];
}
}
return dp[u][0];
}
void query(int tar,int &res,int l,int r,int cur)
{
int m;
res|=val[cur];
if(l==r) return;
m=(l+r)/2;
if(tar<=m) query(tar,res,l,m,2*cur);
else query(tar,res,m+1,r,2*cur+1);
}
void update(int pl,int pr,int l,int r,int cur)
{
int m;
if(pl<=l&&r<=pr)
{
val[cur]=1;
return;
}
m=(l+r)/2;
if(pl<=m) update(pl,pr,l,m,2*cur);
if(pr>m) update(pl,pr,m+1,r,2*cur+1);
}
int main()
{
int i,u,v,resu,resv,ans;
while(scanf("%d",&n)!=EOF)
{
n++;
memset(first,-1,sizeof(first));
num=0;
for(i=1;i<=n-1;i++)
{
scanf("%d%d",&u,&v);
u++,v++;
addedge(u,v);
addedge(v,u);
}
solve();
scanf("%d",&q);
for(i=1;i<=q;i++)
{
scanf("%d%d",&pre[i].u,&pre[i].v);
pre[i].u++,pre[i].v++;
pre[i].lca=getlca(pre[i].u,pre[i].v);
}
sort(pre+1,pre+q+1,cmp);
for(i=1;i<=4*n;i++) val[i]=0;
ans=0;
for(i=1;i<=q;i++)
{
resu=0,resv=0;
query(mp[pre[i].u],resu,1,n,1);
query(mp[pre[i].v],resv,1,n,1);
if(!resu&&!resv)
{
update(mp[pre[i].lca],mp[pre[i].lca]+sum[pre[i].lca]-1,1,n,1);
ans++;
}
}
printf("%d\n",ans);
}
return 0;
}