标签：struct int luogu yyy fa P2597 ZJOI2012 灭绝 define
题面传送门
这道题形象地给我们展示了DAG下支配树地求法。
我们建立一棵树，每个点\(u\)的父亲\(fa_u\)表示如果\(fa_u\)灭绝了，\(u\)一定灭绝。容易发现答案就是子树节点个数-1
考虑这个东西怎么求。
我们从入度为\(0\)的点出发，对于一个点，它的父亲是在DAG上所有儿子在支配树上lca，因为只有这些都灭绝了他才会灭绝。
然后只要搞一个支持加点的倍增就好了。时间复杂度\(O((n+m)logn)\)
code:
#include<bits/stdc++.h>
#define I inline
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define abs(x) ((x)>0?(x):-(x))
#define re register
#define ll long long
#define db double
#define N 70000
#define K 150
#define mod 1000000007
#define eps (1e-9)
#define U unsigned int
#define it iterator
#define Gc() getchar()
#define Me(x,y) memset(x,y,sizeof(x))
#define d(x,y) (n*(x-1)+(y))
using namespace std;
int n,m,k,x,y,in[N+5],fa[N+5][20],lcas,d[N+5],siz[N+5],lg[N+5];
struct yyy{int to,z;}tmp;queue<int> Q;I void swap(int &x,int &y){x^=y^=x^=y;}
struct ljb{int head,h[N+5];yyy f[N+5<<4];I void add(int x,int y){f[++head]=(yyy){y,h[x]};h[x]=head;}}S1,S2,G;
I int lca(int x,int y){d[x]<d[y]&&(swap(x,y),0);while(d[x]^d[y]) x=fa[x][lg[d[x]-d[y]]];if(x==y) return x;for(int i=15;~i;i--)fa[x][i]&&fa[y][i]&&fa[x][i]^fa[y][i]&&(x=fa[x][i],y=fa[y][i]);return fa[x][0];}
I void dfs(int x){siz[x]=1;yyy tmp;for(int i=G.h[x];i;i=tmp.z) tmp=G.f[i],dfs(tmp.to),siz[x]+=siz[tmp.to];}
int main(){
freopen("1.in","r",stdin);
re int i;scanf("%d",&n);for(i=2;i<=n+1;i++) lg[i]=lg[i/2]+1;for(i=1;i<=n;i++) {scanf("%d",&x);while(x) S1.add(i,x),S2.add(x,i),scanf("%d",&x),in[i]++;}for(i=1;i<=n;i++) !in[i]&&(S1.add(i,n+1),S2.add(n+1,i),in[i]=1);Q.push(n+1);while(!Q.empty()){
x=Q.front();Q.pop();for(lcas=0,i=S1.h[x];i;i=tmp.z)tmp=S1.f[i],lcas=(lcas?lca(lcas,tmp.to):tmp.to);fa[x][0]=lcas;d[x]=d[lcas]+1;for(i=1;fa[x][i-1];i++) fa[x][i]=fa[fa[x][i-1]][i-1];G.add(lcas,x);for(i=S2.h[x];i;i=tmp.z)tmp=S2.f[i],in[tmp.to]--,!in[tmp.to]&&(Q.push(tmp.to),0);
}dfs(n+1);for(i=1;i<=n;i++) printf("%d\n",siz[i]-1);
}
标签：struct,int,luogu,yyy,fa,P2597,ZJOI2012,灭绝,define 来源： https://www.cnblogs.com/275307894a/p/15121419.html
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。