灾难 BZOJ 2815

灾难

【样例输入】

5

0

1 0

1 0

2 3 0

2 0

【样例输出】

4

1

0

0

0


题解:

先跑出拓扑序

我们按拓扑序建立一棵“灭绝树”

灭绝树含义是当一个点灭绝时,它的子树将会全部灭绝

所以答案就是点在灭绝树中的子树大小

一个点如果灭绝,那么需要所有指向它的点灭绝

由于拓扑序的关系,指向它的点已经加入过了"灭绝树”中

所以这个点要灭绝,就需要所有指向它的点全部灭绝,即这些点的最近公共祖先

那么直接我们将这个祖先与此点连边,更新Lca

最后求出子树大小,即统计答案

  1 #include<algorithm>
  2 #include<iostream>
  3 #include<cstring>
  4 #include<cstdlib>
  5 #include<cstdio>
  6 #include<cmath>
  7 using namespace std;
  8 inline int Get()
  9 {
 10     int x = 0;
 11     char c = getchar();
 12     while('0' > c || c > '9') c = getchar();
 13     while('0' <= c && c <= '9')
 14     {
 15         x = (x << 3) + (x << 1) + c - '0';
 16         c = getchar();
 17     }
 18     return x;
 19 }
 20 const int me = 1000233;
 21 int n;
 22 int head, tail;
 23 int in[me];
 24 int ue[me];
 25 int de[me];
 26 int si[me];
 27 int fat[me][21];
 28 int tot, nex[2][me], fir[2][me], to[2][me];
 29 inline void Ins(int x, int y, int z)
 30 {
 31     nex[z][++tot] = fir[z][x];
 32     fir[z][x] = tot;
 33     to[z][tot] = y;
 34 }
 35 inline void Topo()
 36 {
 37     head = 0, tail = 0;
 38     for(int i = 1; i <= n; ++i)
 39         if(!in[i])
 40             ue[++tail] = i;
 41     while(head < tail)
 42     {
 43         int u = ue[++head];
 44         for(int i = fir[0][u]; i; i = nex[0][i])
 45         {
 46             int v = to[0][i];
 47             --in[v];
 48             if(!in[v]) ue[++tail] = v;
 49         }
 50     }
 51 }
 52 inline int Lca(int x, int y)
 53 {
 54     if(x < 0) return y;
 55     if(de[x] < de[y]) swap(x, y);
 56     for(int i = 20; i >= 0; --i)
 57         if(de[fat[x][i]] >= de[y])
 58             x = fat[x][i];
 59     for(int i = 20; i >= 0; --i)
 60         if(fat[x][i] != fat[y][i])
 61         {
 62             x = fat[x][i];
 63             y = fat[y][i];
 64         }
 65     if(x == y) return x;
 66     return fat[x][0];
 67 }
 68 inline void Update(int u, int v)
 69 {
 70     fat[v][0] = u;
 71     de[v] = de[u] + 1;
 72     for(int i = 1; i <= 20; ++i)
 73         fat[v][i] = fat[fat[v][i - 1]][i - 1];
 74 }
 75 inline void Build()
 76 {
 77     while(tail)
 78     {
 79         int u = ue[tail];
 80         int lca = -1;
 81         for(int i = fir[0][u]; i; i = nex[0][i])
 82         {
 83             int v = to[0][i];
 84             lca = Lca(lca, v);
 85         }
 86         if(lca < 0) lca = 0;
 87         Ins(lca, u, 1);
 88         Update(lca, u);
 89         --tail;
 90     }
 91 }
 92 void Ergo(int u)
 93 {
 94     si[u] = 1;
 95     for(int i = fir[1][u]; i; i = nex[1][i])
 96     {
 97         int v = to[1][i];
 98         Ergo(v);
 99         si[u] += si[v];
100     }
101 }
102 int main()
103 {
104     n = Get();
105     for(int i = 1; i <= n; ++i)
106     {
107         int x = Get();
108         while(x)
109         {
110             ++in[x];
111             Ins(i, x, 0);
112             x = Get();
113         }
114     }
115     Topo();
116     Build();
117     Ergo(0);
118     for(int i = 1; i <= n; ++i)
119         printf("%d\n", si[i] - 1);
120 }

 

转载于:https://www.cnblogs.com/lytccc/p/6253270.html

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/402346.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

centos关于”running yum-complete-transaction first...

2019独角兽企业重金招聘Python工程师标准>>> 今天在用yum安装软件出错几次户&#xff0c;总是有提示信息&#xff1a; There are unfinished transactions remaining. You might consider running yum-complete-transaction first to finish them. The program yum…

js身份证号、电话脱敏处理(用*替换中间数据)

数字类型 certificatecodecopy certificatecode.replace(/^(.{6})(?:\d)(.{4})$/, "\$1****\$2"); 所有类型 enginenocopy engineno.replace(/^(.{2})(?:\w)(.{1})$/, "\$1****\$2"); enginenocopy engineno.replace(/^(.{4})(?:\w)(.{4})$/, &…

Linux 系统应用编程——网络编程(I/O模型)

Unix下可用的5种I/O模型&#xff1a;阻塞I/O非阻塞I/OI/O复用(select和poll)信号驱动I/O(SIGIO)异步I/O(POSIX的aio_系列函数)一个输入操作通常包括两个不同的阶段&#xff1a;1&#xff09;等待数据准备好&#xff1b;2&#xff09;从内核向进程复制数据&#xff1b;对于一个套…

Windows 7 OpenGL配置

Windows 7 OpenGL配置&#xff0c;解决“无法启动此程序,因为计算机中丢失glut32.dll。”转载于:https://www.cnblogs.com/yangai/p/6253332.html

servlet乱码 解决方法 2种方法

public class ResponseDemo1 extends HttpServlet {public void doGet(HttpServletRequest req, HttpServletResponse resp) throws ServletException, IOException {test1(resp);}//方法1:public void test1(HttpServletResponse resp) throws IOException, UnsupportedEncod…

C#通过属性名字符串获取、设置对象属性值

目录 #通过反射获取对象属性值并设置属性值 0、定义一个类1、通过属性名&#xff08;字符串&#xff09;获取对象属性值2、通过属性名&#xff08;字符串&#xff09;设置对象属性值#获取对象的所有属性名称及类型#判断对象是否包含某个属性回到顶部 #通过反射获取对象属性值…

Linux 系统应用编程——网络编程(TCP/IP 数据包格式解析)

图中括号中的数字代表的是当前域所占的空间大小&#xff0c;单位是bit位。 黄色的是数据链路层的头部&#xff0c;一共14字节 绿色的部分是IP头部&#xff0c;一般是20字节 紫色部分是TCP头部&#xff0c;一般是20字节 最内部的是数据包内容 黄色部分&#xff1a;链路层 目的MA…

mongodb防火墙配置

http://ruby-china.org/topics/20128 https://docs.mongodb.com/manual/tutorial/configure-linux-iptables-firewall/转载于:https://www.cnblogs.com/diyunpeng/p/6256928.html

【python】动态调用函数名

环境&#xff1a; C:\Users\DELL\Desktop>python -V Python 3.9.10 源码&#xff1a; #!/bin/env python # encoding utf-8 import sys import socket# 获取本机ip地址 def get_host_ip():try:s socket.socket(socket.AF_INET, socket.SOCK_DGRAM)s.connect((8.8.8.8, …

AC Again hdoj 1582 搜索

AC Again Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 306 Accepted Submission(s): 84 Description “AC”这个单词大概是我们 ACMers 最希望听到的回答了。不过如果把它写下来&#xff0c;你还能认识吗&am…

Linux 系统应用编程——网络编程(利用TCP/IP 模型分析数据传输过程)

TCP/IP参考模型是一个非常基础&#xff0c;而且也非常重要的基础框架&#xff0c;要想入门数通这是个必须掌握的基本概念&#xff0c;本文档通过一个简单的示例&#xff0c;结合参考模型来分析一下数通的基本过程。网络环境非常简单&#xff0c;如下图所示&#xff0c;我们现在…

Python SQLAlchemy --3

本文為 Python SQLAlchemy ORM 一系列教學文&#xff1a; 刪除 學會如何查詢之後&#xff0c;就能夠進行後續的刪除、更新等操作。 同樣地&#xff0c;以幾個範例做為學習的捷徑。 123456789 user_1 User(user1, username1, password_1)user_2 User(user2, username2, passw…

C#DateTime.ToString格式化时间字符串

DateTime.ToString(string format) 中 format 参数&#xff08;区分大小写&#xff09;&#xff1a; 参数 含义 d 月中的某一天。一位数的日期没有前导零。 dd 月中的某一天。一位数的日期有一个前导零。 ddd 周中某天的缩写名称&#xff0c;在 AbbreviatedDayNam…

Wordpress中显示页面当前位置

Wordpress中显示页面当前位置 在浏览网站的时候&#xff0c;很多网站都提供了页面定位的功能&#xff0c;一般在文章的上方以这样的形式显示&#xff1a; 首页 >> 分类 >> 子分类 >> „ >> 阅读正文 那如何在wordpress中实现这样一个功能呢&#xff1…

CentOS8下 nginx+php7.4安装pdo_sqlsrv 连接 SQL SERVER2012数据库

参考自&#xff1a;https://blog.csdn.net/lwprain/article/details/109598560 关键步骤是一二三&#xff0c;结合 pecl install 正常安装(先安装dnf install unixODBC-devel)即可搞定 ------------------------------------------------------------------------------------…

ArcGIS Engine 编辑- ITask

转自原文ArcGIS Engine 编辑- ITask 下面的代码是我们定制的一个工作流-给等高线赋值 namespace EngineApplication { [Guid("5b0c0692-eaf7-4d64-9cee-c8c1afaf06f4")] [ClassInterface(ClassInterfaceType.None)] [ProgId("EditeTest.ContourTask")] pu…

Linux 系统应用编程——网络编程(TCP 协议三次握手过程)

TCP(Transmission Control Protocol) 传输控制协议 TCP是主机对主机层的传输控制协议&#xff0c;提供可靠的连接服务&#xff0c;采用三次握手确认建立一个连接: 位码即tcp标志位,有6种标示: SYN ( synchronous 建立联机 ) ACK ( acknowledgement 确认 ) PSH ( push 传送…

知识管理在客户关系管理中的应用领域

CRM中的知识获取 CRM前端业务系统&#xff08;包括销售管理、营销管理、客户服务三个模块&#xff09;为与客户互动、获取客户知识提供了各种途径&#xff0c;如销售自动化系统、电子商务网站、呼叫中心、网络社区、博客等等&#xff0c;另外&#xff0c;企业还直接从外部环境、…

基于ELK的简单数据分析

原文链接&#xff1a; http://www.open-open.com/lib/view/open1455673846058.html 环境 CentOS 6.5 64位JDK 1.8.0_20Elasticsearch 1.7.3LogStash 1.5.6Kibana 4.1.4介绍 ElasticSearch是有名的开源搜索引擎&#xff0c;现在很多公司使用ELK技术栈做日志分析&#xff0c;比如…

Linux 系统应用编程——网络编程(高级篇)

一、网络超时检测 在网络通信过程中&#xff0c;经常会出现不可预知的各种情况。例如网络线路突发故障、通信一方异常结束等。一旦出现上述情况&#xff0c;很可能长时间都不会收到数据&#xff0c;而且无法判断是没有数据还是数据无法到达。如果使用的是TCP协议&#xff0c;可…