2016-2017 ACM-ICPC, NEERC, Northern Subregional Contest G Gangsters in Central City (LCA)
题意:
有一棵树,水源在根节点1,房子在叶子节点。有若干操作,操作可能是歹徒占领或者离开一个房子。我们不想给歹徒供水,可以通过切断边实现(如果某个叶子节点到根节点的路径上有一条边被切掉,那么就不能供水了。)对于每次操作后,问不给所有歹徒供水最少要切多少条边,并且问在切满足前面最小的情况下,最少使得多少个良民受影响。初始没有歹徒。
思路:
我们先考虑第一个问题。容易知道,假设与根相连的有k条边,那么最多只需要切k次,就切断了所有房子的水源。
也就是说,从最少切的次数角度的考虑,切与水源相连的边一定是最优的。
我们可以考虑把树根1去掉,这样得到k棵子树
然后可以预处理出,对于每个叶子节点,其非根的最远祖先是谁,也就是k棵子树的根节点都是谁。
那么对于每次出现歹徒,假设其非根的最远祖先是x,只需要切1-x这条边即可。
现在考虑第二个问题,在保证问题一最小的情况下,一个比较直观的想法是,我们尽量往低了切,也就是尽量往原理根的边上切,这样才能使收到影响的良民比较少。
容易想到,一个子树上该切的点是,所有被歹徒占领的坏点的LCA。这样可以使得受到影响的良民最少。
因为我们要得到受到影响的良民的数量,所以用siz[i]维护以i为根的子树的叶子数量,以及歹徒的数量。
这道题的关键结论是,:树上多个点的LCA,就是DFS序最小的和DFS序最大的这两个点的LCA"
这道题的关键结论是,:树上多个点的LCA,就是DFS序最小的和DFS序最大的这两个点的LCA"
这道题的关键结论是,:树上多个点的LCA,就是DFS序最小的和DFS序最大的这两个点的LCA"
因此这题就是写个LCA就可以了,切掉的坏点的dfs序可以用个set维护下。
/* ***********************************************
Author :111qqz
Created Time :2017年10月12日 星期四 17时06分57秒
File Name :G.cpp
************************************************ */
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <string>
#include <cmath>
#include <cstdlib>
#include <ctime>
#define PB push_back
#define fst first
#define sec second
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define ms(a,x) memset(a,x,sizeof(a))
typedef long long LL;
#define pi pair < int ,int >
#define MP make_pair
using namespace std;
const double eps = 1E-8;
const int dx4[4]={1,0,0,-1};
const int dy4[4]={0,-1,1,0};
const int inf = 0x3f3f3f3f;
const int N=1E5+7;
int n,q;
int val[N];
vector < int > edge[N];
int in[N];
int E[2*N],R[2*N],dis[N],depth[2*N];
int p;
int fa[N];
int dp[2*N][20];
int siz[N]; //siz[i]表示以i为根的子树的叶子数量。
int FA[N]; //将树根去掉之后,每棵子树最上面能到达的顶点。
set<int>T[N];//T[i]表示与根节点直接连接的节点i 切掉的点的dfs序 id
int a[N];
void dfs( int u,int dep,int d,int pre)
{
fa[u] = pre;
p++;
E[p] = u;
depth[p] = dep;
R[u] = p;
dis[u] = d;
int SIZ = edge[u].size();
if (SIZ==0) siz[u] = 1;
for ( int i = 0 ; i < SIZ ; i++)
{
int v = edge[u][i];
if (v==pre) continue;
if (u==1) FA[v] = v; //把树根去掉,原树分裂乘若干个子树
else FA[v] = FA[u];
dfs(v,dep+1,d+1,u);
siz[u] += siz[v];
p++;
E[p] = u;
depth[p] = dep;
}
}
int _min( int l,int r)
{
if (depth[l] < depth[r]) return l;
return r;
}
void rmq_init()
{
for ( int i = 1 ; i <= 2*n+2 ; i++) dp[i][0] = i;
for ( int j = 1 ; (1<<j) <= 2*n+2 ; j++)
for ( int i = 1 ; i + (1<<j)-1 <= 2*n+2 ; i++)
dp[i][j] = _min(dp[i][j-1],dp[i+(1<<(j-1))][j-1]);
}
int rmq_min(int l,int r)
{
if (l>r) swap(l,r);
int k = 0 ;
while (1<<(k+1)<=r-l+1) k++;
return _min(dp[l][k],dp[r-(1<<k)+1][k]);
}
int ans1=0,ans2=0;
int main()
{
#ifndef ONLINE_JUDGE
//freopen("./in.txt","r",stdin);
#endif
freopen("gangsters.in","r",stdin);
freopen("gangsters.out","w",stdout);
ms(a,0);
ms(siz,0);
scanf("%d %d",&n,&q);
for ( int i = 2 ; i <= n ; i++)
{
int x;
scanf("%d",&x);
// cout<<"x:"<<x<<endl;
edge[x].PB(i);
}
p = 0;
dfs(1,0,0,-1);
rmq_init();
while (q--)
{
char opt[3];
int x;
scanf("%s%d",opt,&x);
// cout<<opt<<" "<<x<<endl;
if (opt[0]=='+')
{
int fx = FA[x];
if (T[fx].size()==0) ++ans1;
T[fx].insert(R[x]);
ans2-=a[fx];//因为切掉了,先去掉之前的影响
int fi = *T[fx].begin();
int la = *T[fx].rbegin();
int LCA = E[rmq_min(fi,la)]; //尽量往低了切,才可能影响小。
// printf( "(LCA(%d %d)=%d,siz[LCA]=%d \n",R[fi],R[la],LCA,siz[LCA]);
a[fx] = siz[LCA] - T[fx].size();
ans2 += a[fx];//再考虑新的影响
// cout<<"ans2:"<<ans2<<endl;
}
else
{
int fx = FA[x];
if (T[fx].size()==1) ans1--;
T[fx].erase(R[x]);
ans2-=a[fx];
if (T[fx].size()!=0)
{
int fi = *T[fx].begin();
int la = *T[fx].rbegin();
int LCA = E[rmq_min(fi,la)];
a[fx] = siz[LCA] - T[fx].size();
// printf( "(LCA(%d %d)=%d\n",R[fi],R[la],LCA);
}
else a[fx] = 0;
ans2 +=a[fx];
}
printf("%d %d\n",ans1,ans2);
}
#ifndef ONLINE_JUDGE
fclose(stdin);
#endif
return 0;
}