http://www.cnblogs.com/JS-Shining/archive/2013/01/12/2857429.html 题面

题解上写了用什么dominator tree,吓晕了,看了看,好像看不懂,于是看了看hzwer的代码,发现和dominator tree无半点关系。

-----------------------------------------------------------------------------------------------------------------------------------

首先我们假设有一个太阳,就是所有没有食物的物种的食物。

然后我们想一下 如果一个物种要灭绝,那么他的食物得先灭绝。他的食物怎么会灭绝?那么他的食物的食物也得灭绝,最后这些食物的根源肯定会汇集到一个点上,也就是所有食物得lca(你想啊一棵树上几个点不断往上爬肯定会汇集到一个点上(其实我也不太懂))。首先我们得拓扑排序一下,因为拓扑排序让一个物种能处于他的食物以及食物的食物之后(之后是指顺序上,也就是保证先处理好了自己食物及食物得的食物等,再处理自己)。然后扫描拓扑排序后的物种,对于每个物种,找他食物的lca,如果原图中没有食物,我们搞的太阳就变成了它的食物,然后把这个物种放在lca下面。这棵树可以帮我们统计答案,因为树的节点的儿子都是这个节点灭绝后会导致灭绝的物种。为什么放在lca下可以呢?因为我们要找lca的目的在于找到哪个物种能让自己灭绝,又因为lca能让这个物种灭绝,所以放在下面是对的。如何统计答案呢?因为每个点下挂的节点都是自己的灭绝能让他也灭绝,自己的儿子也是这样,所以以每个节点的子树大小-1就是答案。-1是把自己减去。

(我的代码是不是跟hzwer的很像?因为我看了他的,深受影响)还说一个可能的问题,为什么两个图能共用一个e数组?因为他们的head不一样

#include<bits/stdc++.h>
using namespace std;
#define N 100010
struct edge
{
int nxt,to;
}e[N<<];
int n,cnt=;
vector<int> q,s;
struct G
{
int dep[N],head[N],size[N],in[N];
int fa[N][];
void init()
{
memset(fa,-,sizeof(fa));
}
void link(int u,int v)
{
e[++cnt].nxt=head[u];
head[u]=cnt;
e[cnt].to=v;
in[v]++;
}
int lca(int u,int v)
{
if(u==-) return v;
if(v==-) return u;
if(dep[u]<dep[v]) swap(u,v);
int temp=dep[u]-dep[v];
for(int i=;i>=;--i) if(temp&(<<i))
u=fa[u][i];
if(u==v) return u;
for(int i=;i>=;--i) if(fa[u][i]!=fa[v][i])
{
u=fa[u][i]; v=fa[v][i];
}
return fa[u][];
}
void dfs(int u)
{
size[u]=;
for(int i=head[u];i;i=e[i].nxt)
{
dfs(e[i].to);
size[u]+=size[e[i].to];
}
}
void tp()
{
for(int i=;i<=n;++i) if(!in[i]) q.push_back(i);
while(!q.empty())
{
int u=q.back(); q.pop_back(); s.push_back(u);
for(int i=head[u];i;i=e[i].nxt)
{
in[e[i].to]--;
if(!in[e[i].to]) q.push_back(e[i].to);
}
}
}
void update(int u)
{
for(int i=;i<=;++i) if(fa[u][i-])
fa[u][i]=fa[fa[u][i-]][i-];
}
}G1,G2;
void build_tree()
{
for(int i=s.size()-;i>=;--i)
{
int u=s[i],lca=-;
for(int j=G1.head[u];j;j=e[j].nxt)
lca=G2.lca(e[j].to,lca);
if(lca==-) lca=;
G2.link(lca,u);
G2.fa[u][]=lca;
G2.dep[u]=G2.dep[lca]+;
G2.update(u);
}
}
int main()
{
scanf("%d",&n);
for(int i=;i<=n;++i)
{
int x; scanf("%d",&x);
while(x)
{
G1.link(i,x);
scanf("%d",&x);
}
}
G1.tp();
build_tree();
G2.dfs();
for(int i=;i<=n;++i)
{
if(G2.size[i]) printf("%d\n",G2.size[i]-);
else puts("");
}
return ;
}

bzoj 2815的更多相关文章

  1. bzoj 2815 [ZJOI2012]灾难(构造,树形DP)

    [题意] 求把每个点删除后,不可达点的数目. [思路] 构造一棵“灭绝树”,要求这棵树满足如果删除根节点后则该子树内的所有结点都不可达.则答案为子树大小-1. 如何构造这棵“灭绝树”? 将原图拓扑排序 ...

  2. bzoj 2815 灭绝树

    对于一个食物网(一个DAG),一个物种死亡后,某些物种就必然死亡,求出必然死亡的是那些物种. 灭绝树的另一种含义是:“灭绝树跟节点到节点u的路径上的节点由那些原图中从根节点到节点u的所有路径中都经过了 ...

  3. 灾难 bzoj 2815

    灾难(1s 128MB)catas [样例输入] 5 0 1 0 1 0 2 3 0 2 0 [样例输出] 4 1 0 0 0 题解: 主要算法:拓扑排序:最近公共祖先(Lca): 先跑出拓扑序 我们 ...

  4. BZOJ 2815: [ZJOI2012]灾难

    呃,题面没了,大概就是给出一些生物之间的捕食关系,求灭绝树每个点的灾难值. 拓扑排序之后倒着加入点,动态维护fa[][]数组,倍增法求LCA,当然大佬愿意写动态树也是极好的…… #include &l ...

  5. bzoj 2815 灾难

    首先假设我们定义x灭绝后y会灭绝,那么离y最近的x就为y的父亲节点,那么如果我们可以求出每个节点的父亲节点,我们就得到了一棵树,然后每个节点的灾难值就是子树的大小-1. 我们将出度数为0的节点的父亲节 ...

  6. BZOJ 2815: [ZJOI2012]灾难 拓扑排序+倍增LCA

    这种问题的转化方式挺巧妙的. Code: #include <bits/stdc++.h> #define N 100000 #define M 1000000 #define setIO ...

  7. BZOJ2815: [ZJOI2012]灾难

    传送门 学LCA的时候根本没意识到LCA可以有这么多玩法. 这玩意据说是个高级数据结构(支配树)的弱化版,蒟蒻没学过呀.所以出题人提出一个概念叫灾难树. 我理解的灾难树的意思实际上是属于DAG的一个子 ...

  8. Luogu 2597 [ZJOI2012]灾难

    BZOJ 2815. 解法还是挺巧妙的. 放上写得很详细很好懂的题解链接  戳这里. 一个物种$x$如果要灭绝,那么沿着它的入边反向走走走,一定可以走到一个点$y$,如果这个点$y$的物种灭绝了,那么 ...

  9. 2815: [ZJOI2012]灾难 - BZOJ

    题目描述 Description 阿米巴是小强的好朋友.    阿米巴和小强在草原上捉蚂蚱.小强突然想,如果蚂蚱被他们捉灭绝了,那么吃蚂蚱的小鸟就会饿死,而捕食小鸟的猛禽也会跟着灭绝,从而引发一系列的 ...

随机推荐

  1. 【板+背包】多重背包 HDU Coins

    http://acm.hdu.edu.cn/showproblem.php?pid=2844 [题意] 给定n种价值为Ci,个数为Wi的硬币,问在1~V中的这些数中哪些数能由这些硬币组成? [思路] ...

  2. JPA框架下使用纯粹的原生SQL

    最近遇到一个需求,查询数据库中对应表的字段是动态的,项目使用的框架使用JPA+Spring Boot,JPA自带原生SQL支持的传入参数是强类型的,无法用于查询语句的字段更改,因为插入字符串的话带有单 ...

  3. No route info of this topic

    使用rocketmq时报错 com.alibaba.rocketmq.client.exception.MQClientException: No route info of this topic, ...

  4. BUPT2017 springtraining(16) #1 题解

    https://vjudge.net/contest/162590 A: 不难发现,当L=R时输出L,当L<R时输出2. B: 贪心得配对.1和n配 2和n-1配,对与对直接只要花1个代价就可以 ...

  5. React学习及实例开发(三)——用react-router跳转页面

    本文基于React v16.4.1 初学react,有理解不对的地方,欢迎批评指正^_^ 一.定义路由 1.安装react-router npm install react-router@ --sav ...

  6. 七天从零基础学习android(3)--实现过程

    首先这是我对自己编写程序的认识,要实现一个程序,根据之前编写C++的经验,要对所编写的软件有一个模糊的了解. 一个记账本软件,要实现的过程是,添加收支,显示本日,本月或本年的收支状态.然而基于是完全没 ...

  7. weblogic负载分发

    博客分类: weblogic 负载均衡的实现方式有很多种,这里只介绍三种相对来说成本较低的方案(维护成本以及费用成本)weblogic自带的proxy.apache.nginx 1.weblogic自 ...

  8. Python访问MySQL数据库并实现其增删改查功能

    概述:对于访问MySQL数据库的操作,我想大家也都有一些了解.不过,因为最近在学习Python,以下就用Python来实现它.其中包括创建数据库和数据表.插入记录.删除记录.修改记录数据.查询数据.删 ...

  9. 【scrapy】创建第一个项目

    1)创建项目命令: scrapy startproject tutorial 该命令将在当前目录下创建tutorial文件夹 2)定义Item Items are containers that wi ...

  10. ACM之数论数字根

    先来看一道杭电的数字根问题 此题的大大意是输入一个数.假设它不是一位的数字的话,那么我们就将它的每一位都相加,相加后假设还是两位或者很多其它的话那么我们继续取出它的每一位数字进行相加.知道等到单个数字 ...