P2597 [ZJOI2012]灾难(倍增LCA+拓扑排序)
传送门
据大佬说这玩意儿好像叫灾难树还是灭绝树?
我们先按建图,设点$u$的食物有$x[1]...x[k]$,即在图中这些点都有一条指向它的边
以样例来说,对于人,羊和牛都有一条指向它的边,然而不管是羊的灭绝还是牛的灭绝都不会导致人的灭绝,会导致人灭绝的只有草,也就是羊和牛的LCA
(严格来说这个应该不叫LCA,不过反正这是一个有向无环图,大家能理解就好)
于是我们把人挂到草的下面,表示只有草的灭绝会使人灭绝。
于是我们按拓扑序把每一个生物都依次加入灭绝树,然后一遍dfs,每一个生物的灾难值就是它的子树的大小减一
1 //minamoto 2 #include<iostream> 3 #include<cstdio> 4 using namespace std; 5 #define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++) 6 char buf[1<<21],*p1=buf,*p2=buf; 7 int read(){ 8 int res,f=1;char ch; 9 while((ch=getc())>'9'||ch<'0')(ch=='-')&&(f=-1); 10 for(res=ch-'0';(ch=getc())>='0'&&ch<='9';res=res*10+ch-'0'); 11 return res; 12 }char sr[1<<21],z[20];int C=-1,Z=0; 13 inline void Ot(){fwrite(sr,1,C+1,stdout),C=-1;} 14 void print(int x){ 15 if(C>1<<20)Ot();if(x<0)sr[++C]='-',x=-x; 16 while(z[++Z]=x%10+48,x/=10); 17 while(sr[++C]=z[Z],--Z);sr[++C]='\n'; 18 } 19 const int N=65555; 20 int head[N<<2],Next[N<<2],ver[N<<2],tot,head1[N<<2],ver1[N<<2],Next1[N<<2],tot1,head2[N<<2],ver2[N<<2],Next2[N<<2],tot2; 21 inline void add(int u,int v){ver[++tot]=v,Next[tot]=head[u],head[u]=tot;} 22 inline void add1(int u,int v){ver1[++tot1]=v,Next1[tot1]=head1[u],head1[u]=tot1;} 23 inline void add2(int u,int v){ver2[++tot2]=v,Next2[tot2]=head2[u],head2[u]=tot2;} 24 int n,du[N],q[N],dep[N],fa[N][35],ans[N]; 25 void topo(){ 26 int h=1,t=0; 27 for(int i=1;i<=n;++i)if(!du[i])q[++t]=i; 28 while(h<=t){ 29 int u=q[h++]; 30 for(int i=head[u];i;i=Next[i]){ 31 int v=ver[i]; 32 if(--du[v]==0)q[++t]=v; 33 } 34 } 35 } 36 int LCA(int u,int v){ 37 if(dep[u]<dep[v])swap(u,v);int k=dep[u]-dep[v]; 38 for(int i=0;i<18;++i) 39 if((k>>i)&1)u=fa[u][i]; 40 if(u==v)return u; 41 for(int i=17;i>=0;--i) 42 if(fa[u][i]!=fa[v][i])u=fa[u][i],v=fa[v][i]; 43 return fa[u][0]; 44 } 45 void dfs(int u){ 46 for(int i=head2[u];i;i=Next2[i]) 47 dfs(ver2[i]),ans[u]+=ans[ver2[i]]; 48 ++ans[u]; 49 } 50 int main(){ 51 // freopen("testdata.in","r",stdin); 52 n=read(); 53 for(int i=1,u;i<=n;++i){ 54 while(true){ 55 if((u=read())==0)break; 56 ++du[i],add1(i,u),add(u,i); 57 } 58 } 59 topo(); 60 for(int i=1;i<=n;++i){ 61 int u=ver1[head1[q[i]]]; 62 for(int j=head1[q[i]];j;j=Next1[j])u=LCA(u,ver1[j]); 63 add2(u,q[i]),dep[q[i]]=dep[u]+1,fa[q[i]][0]=u; 64 for(int j=1;j<18;++j)fa[q[i]][j]=fa[fa[q[i]][j-1]][j-1]; 65 } 66 dfs(0); 67 for(int i=1;i<=n;++i)print(ans[i]-1); 68 return Ot(),0; 69 }
转载于:https://www.cnblogs.com/bztMinamoto/p/9888165.html
P2597 [ZJOI2012]灾难(倍增LCA+拓扑排序)相关推荐
- P2597 [ZJOI2012]灾难(top+lca)
阿米巴是小强的好朋友. 阿米巴和小强在草原上捉蚂蚱.小强突然想,果蚂蚱被他们捉灭绝了,那么吃蚂蚱的小鸟就会饿死,而捕食小鸟的猛禽也会跟着灭绝,从而引发一系列的生态灾难. 学过生物的阿米巴告诉小强,草原 ...
- [bzoj2815] [洛谷P2597] [ZJOI2012] 灾难
Description 阿米巴是小强的好朋友. 阿米巴和小强在草原上捉蚂蚱.小强突然想,如果蚂蚱被他们捉灭绝了,那么吃蚂蚱的小鸟就会饿死,而捕食小鸟的猛禽也会跟着灭绝,从而引发一系列的生态灾难. 学过 ...
- P2597 [ZJOI2012]灾难(支配树)
题目描述 阿米巴是小强的好朋友. 阿米巴和小强在草原上捉蚂蚱.小强突然想,如果蚂蚱被他们捉灭绝了,那么吃蚂蚱的小鸟就会饿死,而捕食小鸟的猛禽也会跟着灭绝,从而引发一系列的生态灾难. 学过生物的阿米巴告 ...
- BZOJ2815:[ZJOI2012]灾难(拓扑排序,LCA)
Description 阿米巴是小强的好朋友. 阿米巴和小强在草原上捉蚂蚱.小强突然想,如果蚂蚱被他们捉灭绝了,那么吃蚂蚱的小鸟就会饿死,而捕食小鸟的猛禽也会跟着灭绝,从而引发一系列的生态灾难. 学过 ...
- jzoj4240 [五校联考5day2]游行 拓扑排序+倍增lca+线段树优化建图
Description 恶梦是学校里面的学生会主席.他今天非常的兴奋,因为学校一年一度的学生节开始啦!! 在这次节日上总共有N个节目,并且总共也有N个舞台供大家表演.其中第i个节目的表演时间为第i个单 ...
- BZOJ2815: [ZJOI2012]灾难
BZOJ2815: [ZJOI2012]灾难 题目描述 阿米巴是小强的好朋友. 阿米巴和小强在草原上捉蚂蚱.小强突然想,如果蚂蚱被他们捉灭绝了,那么吃蚂蚱的小鸟就会饿死,而捕食小鸟的猛禽也会跟着灭绝, ...
- BZOJ 2144 跳跳棋(神仙建模题,倍增 LCA,二分)【BZOJ修复工程】
整理的算法模板合集: ACM模板 点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划 题目链接 https://hydro.ac/d/bzoj/p/2144 是 hydro 的 BZOJ ...
- 【基础练习】【拓扑排序】codevs3294 车站分级题解
题目来源:NOIP2013 普及第四题 题目描写叙述 Description 一条单向的铁路线上,依次有编号为1, 2, -, n的n个火车站.每一个火车站都有一个级别,最低为1级.现有若干趟车次在这 ...
- hihocoder 1343 : Stable Members【拓扑排序】
hihocoder #1343:题目 解释: 一个学习小组,一共有N个学员,一个主管.每个学员都有自己的导师(一个或者多个),导师可以是其他学员也可以是主管. 每周学员都要把自己的学习报告和收到的报告 ...
最新文章
- python转义符个数,python(五)——运算符、字符串、转义字符
- P7011-[CERC2013]Escape【堆,启发式合并】
- android按钮变输入框动画,AnimShopButton 仿饿了么加入购物车旋转控件 - 自带闪转腾挪动画 的按钮。自带海量可定制 UI 属性。在 Re @codeKK Android开源站...
- 解决django运行中连接mysql数据库超时报错pymysql.err.InterfaceError
- SIFT特征检测算子和sift = cv2.xfeatures2d.SIFT_create出错的解决办法
- 产品原型工具 AXURE9 Mac汉化
- 电源硬件设计----开关电源布局设计要领
- 文件系统性能测试工具 FIO工具
- Python - 使用Pylint检查分析代码
- Photoshop的安装教程
- 联发科mt8516价格_揭秘联发科MT8516单颗芯片破千万背后的故事
- 智能点餐mysql框架图_SpringBoot微信点餐系统--P3数据库设计
- ios开发的p12和provision
- DBLP数据集用weka数据挖掘 xml转csv格式文件
- 内外墙乳胶漆颜色怎么搭配,内外墙品牌涂料如何选购?
- 抖音seo源码,抖音矩阵,抖音seo系统,抖音搜索排名
- 数学速算法_小学初中高中 数学奥数教材及习题讲解(共42册PDF)
- 第十四课:树莓派网络视频监控(mjpg-streamer)
- 日更100天(5)每天进步一点点
- NGUI之UITable
热门文章
- SAP CRM Business Partner API的buffer工作原理
- 频域分析方法的理解(信号、频谱、能量谱、功率谱、倒频谱、小波分析)(终极教程)
- java并发排序_Java并发(三):重排序
- android 彩蛋 miui,MIUI12最新更新,安卓11彩蛋终于出现,流畅度提升
- java redis缓存实例_spring项目整合ehcache和redis缓存实例
- arraylist转int数组_LeetCode 1389. 按既定顺序创建目标数组
- python的scripts里没有pip_解决python scripts和pip缺失问题
- mysql授权无密码用户_MySQL下新建用户,授权,删除用户,修改密码
- java 管理员命令模式_java中命令模式详解和使用方法
- oracle查询orapw文件,oracle学习笔记《一》