[IOI2008/BZOJ1791 岛屿](处理基环树的小技巧&基于bfs树形DP)
题目大意是在一个基环树森林里求每一棵基环树的直径①的和。
其实就是树的直径的基环树升级版。我们先把环找出来,然后从环上的每一个节点x出发,并且不经过环上其他节点,做一次树形DP,求出x的子树中到x最远的路径长d[x]和x的子树的直径dp[x]。
那么基环树的直径只有两种情况:
1.直径不经过环上的节点
对于每一个dp[x]取max就好了。
2.直径经过环上的节点
对于任意环上的两个节点x和y,都可以构成一条路径d[x]+d[y]+dist(x,y),其中dist(x,y)为x和y在环上的距离②。
可以直接n^2枚举找最大的吗?当然不能,N<=1e6耶。
这个时候就要用一点套路。首先把环复制一遍接在原先的环后,然后对dist数组做个前缀和s,那么一条路径就成了(x>y) d[x]+d[y]+s[x]-s[y]=(d[x]+s[x])+(d[y]-s[y]),也就是说对于每个x我们需要找到在x之前的一个y满足(d[y]-s[y])最大,并且x-y<=环长。单调队列优化的套路嘛!
①基环树的直径:即基环树中最长的简单路径,这里简单路径指的是不自交(不重复经过任何点或边)的路径。
②环上的距离:有顺时针和逆时针两种,因为求的是直径,所以取两者中较长的。
这是luogu的AC代码:(温馨提示:代码太丑并且不是重点,您大可直接跳过看后面)
#include<bits/stdc++.h>
using namespace std;
const int N=1e6+;
typedef long long LL;
struct edge{
int v,w,last;
}e[N<<];
int tot=,tail[N];
inline void add(int x,int y,int z){
e[++tot]=(edge){y,z,tail[x]};
tail[x]=tot;
}
inline int read(){
int x=,w=;char ch=;
while(!isdigit(ch)) w|=ch=='-',ch=getchar();
while(isdigit(ch)) x=(x<<)+(x<<)+(ch^),ch=getchar();
return w?-x:x;
}
bool incur[N];
int n,dn,cnt,dfn[N],tp,st[N][],c[N],s[N];
void getcur(int x,int pre){
dfn[x]=++dn;
for(int p=tail[x];p;p=e[p].last)if(p^pre^){
int &v=e[p].v,&w=e[p].w;
if(dfn[v]){
if(dfn[v]>=dfn[x]) continue;
c[cnt=]=v,s[cnt]=w,incur[v]=true;
for(int i=tp;st[i][]^v;--i) c[++cnt]=st[i][],s[cnt]=st[i][],incur[st[i][]]=true;
}
else st[++tp][]=v,st[tp][]=w,getcur(v,p);
}
--tp;
} int hd,rt,tl,q[N<<];
LL ans,t,d[N],f[N],dp[N],sum[N<<];
LL dfs(int x,int pre){
LL t=;d[x]=f[x]=;
for(int p=tail[x];p;p=e[p].last)if(p^pre^){
int &v=e[p].v,&w=e[p].w;
if(incur[v]||v==pre) continue;
t=max(t,dfs(v,x)),f[x]=max(f[x],d[x]+d[v]+w),d[x]=max(d[x],d[v]+w);
}
return max(f[x],t);
}
LL F(int id){
return dp[id]-sum[id];
}
int main(){
n=read();
for(int i=;i<=n;++i){
int x=read(),y=read();
add(x,i,y),add(i,x,y);
}
for(int i=;i<=n;++i)if(!dfn[i]){
memset(sum,,sizeof sum);
st[tp=][]=i,st[tp][]=;
cnt=;getcur(i,);
t=;
for(int j=;j<=cnt;++j) t=max(t,dfs(c[j],)),dp[j]=dp[j+cnt]=d[c[j]];
int hd=tl=;
for(int j=;j<=cnt*;++j){
sum[j]+=sum[j-]+s[(j-)>cnt?j--cnt:j-];
while((hd^tl)&&j-q[hd]+>cnt) ++hd;
if(hd^tl) t=max(t,F(q[hd])+dp[j]+sum[j]) ;
while((hd^tl)&&F(j)>F(q[tl-])) --tl;
q[tl++]=j;
}
ans+=t;
}
cout<<ans<<endl;
return ;
}
又是找环,又是树形DP,又是单调队列优化的,是不是很恶心?没错,虽然用的算法不算高级,但由于种种原因做此题的人各种RE&WA&MLE&TLE,因此在洛谷上都成黑题了。
并且我要告诉你更恶心的是,由于BZOJ栈的限制,上面这份dfs的代码在BZOJ上还会爆栈......
正是在这种山穷水尽的时候,我发现了一些黑科技。
首先我们要分析题目图的特性,这是一棵无向基环树(暂时不管森林),画出来应该是这个样子(没画边权,下同):
题目的输入方式有点特别,是第i+1行输入一个x和y表示i+1和x有一条长度为y的边。那么我们暂时把它当作有向边,令to[i+1]=x,w[i+1]=y,画出来就是这个样子:
它刚好是一棵内向树③!这是巧合吗?不是的。按照这种输入方式,每个点有且仅有一条出边,连出来就一定是一棵内向树。这样一来,我们就得到了每个节点x的父亲f[x]与它到父亲的边权w[x],当然对于环上的节点是到它相邻节点的,且一定是同一个方向(有点不太好表达,看代码就懂了)。于是我们就可以用按拓扑排序的顺序自底向上地DP了,也就是基于bfs的树形DP。代码:
③内向树:n个点,n条边,每个节点有且仅有一条出边的有向连通图就好像以”基环“为中心,有向内收缩的趋势,故称为”内向树“。
#include<bits/stdc++.h>
using namespace std;
#define rg register
const int N=1e6,M=1e6+;
typedef long long LL;
inline int read(){
int x=,w=;char ch=;
while(ch<''||ch>'') w|=ch=='-',ch=getchar();
while(ch>=''&&ch<='') x=(x<<)+(x<<)+(ch^),ch=getchar();
return w?-x:x;
}
bool v[M];
int n,cnt,dcnt,cur[M<<],to[M],w[M],in[M];
LL t,ans,s[M<<],d[M],f[M],dp[M];
struct Deque{//这个。。。手写的双端队列,不用管它
int l,r,q[N];
Deque(){l=r=;}
void clear(){l=r=;}
bool empty() {return !(l^r);}
void push_back(int v) {q[r++]=v;r%=N;}
void push_front(int v) {q[l=(l-+N)%N]=v;}
void pop_front() {++l;l%=N;}
void pop_back() {r=(r-+N)%N;}
int front() {return q[l];}
int back() {return q[r-];}
}q;
int main(){
n=read();
for(rg int i=;i<=n;++i) ++in[to[i]=read()],w[i]=read();
for(rg int i=;i<=n;++i) if(!in[i]) q.push_back(i);//入度为0的点入队
while(!q.empty()){
int x=q.front(),&fa=to[x],&l=w[x];q.pop_front();
if(!--in[fa]) q.push_back(fa);
//下面3行是DP,也就是更新父节点的信息
f[fa]=max(f[fa],d[fa]+d[x]+l); //f:x的子树中经过x节点的最长路径
d[fa]=max(d[fa],d[x]+l);//d: x的子树中离x最远的节点到x的距离
dp[fa]=max(dp[fa],max(dp[x],f[fa])); //dp:x的子树的直径,也就是对f[x]和dp[y](to[y]=x)取max
}
for(rg int i=;i<=n;++i)if(in[i]&&!v[i]){//基环树森林,不只一个环,用v标记
cnt=;t=;
for(rg int j=i;!v[j];j=to[j]){//这就是上面说的,这样一直沿着to走,一定能绕环一圈
cur[++cnt]=j,v[j]=true;//把环抠下来
t=max(t,dp[j]);//上面讨论的基环树直径的第一种情况
}
dcnt=cnt<<;
for(int i=cnt+;i<=dcnt;++i) cur[i]=cur[i-cnt];//复制一遍环
q.clear();//注意清空
for(rg int j=;j<=dcnt;++j){
s[j]=s[j-]+w[cur[j-]];//前缀和
while(!q.empty()&&j-q.front()>=cnt) q.pop_front();//使满足j-q.front()+1<=cnt
if(!q.empty())
t=max(t,d[cur[q.front()]]-s[q.front()]+d[cur[j]]+s[j]);
while(!q.empty()&&d[cur[j]]-s[j]>d[cur[q.back()]]-s[q.back()]) q.pop_back();
q.push_back(j);
}
ans+=t;//累加每一棵基环树的直径
}
printf("%lld\n",ans);
return ;
}
于是我们仅用这么点简洁的代码就A了一道黑题!这道题告诉我:拿到题先分析题目的特性,不要一来就按照传统方式存图,找环,dfs树形DP......
2019-06-26
[IOI2008/BZOJ1791 岛屿](处理基环树的小技巧&基于bfs树形DP)的更多相关文章
- bzoj1791[IOI2008]Island岛屿(基环树+DP)
题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=1791 题目大意:给你一棵n条边的基环树森林,要你求出所有基环树/树的直径之和.n< ...
- bzoj 1791: [Ioi2008]Island 岛屿【基环树+单调队列优化dp】
我太菜了居然调了一上午-- 这个题就是要求基环树森林的基环树直径和 大概步骤就是找环->dp找每个环点最远能到达距离作为点权->复制一倍环,单调队列dp 找环是可以拓扑的,但是利用性质有更 ...
- luogu 4381 [IOI2008]Island 单调队列 + 基环树直径 + tarjan
Description 你将要游览一个有N个岛屿的公园.从每一个岛i出发,只建造一座桥.桥的长度以Li表示.公园内总共有N座桥.尽管每座桥由一个岛连到另一个岛,但每座桥均可以双向行走.同时,每一对这样 ...
- 51nod 1673 树有几多愁(链表维护树形DP+状压DP)
题意 lyk有一棵树,它想给这棵树重标号. 重标号后,这棵树的所有叶子节点的值为它到根的路径上的编号最小的点的编号. 这棵树的烦恼值为所有叶子节点的值的乘积. lyk想让这棵树的烦恼值最大,你只需输出 ...
- 51nod 1812 树的双直径 题解【树形DP】【贪心】
老了-稍微麻烦一点的树形DP都想不到了. 题目描述 给定一棵树,边权是整数 \(c_i\) ,找出两条不相交的链(没有公共点),使得链长的乘积最大(链长定义为这条链上所有边的权值之和,如果这条链只有 ...
- 50分钟学会Laravel 50个小技巧(基于laravel5.2,仅供参考)
转载请注明:转载自 Yuansir-web菜鸟 | LAMP学习笔记 本文链接地址: 50分钟学会Laravel 50个小技巧 原文链接:< 50 Laravel Tricks in 50 Mi ...
- BZOJ4446:[SCOI2015]小凸玩密室(树形DP)
Description 小凸和小方相约玩密室逃脱,这个密室是一棵有n个节点的完全二叉树,每个节点有一个灯泡.点亮所有灯泡即可逃出密室. 每个灯泡有个权值Ai,每条边也有个权值bi.点亮第1个灯泡不需要 ...
- 模拟赛:树和森林(lct.cpp) (树形DP,换根DP好题)
题面 题解 先解决第一个子问题吧,它才是难点 Subtask_1 我们可以先用一个简单的树形DP处理出每棵树内部的dis和,记为dp0[i], 然后再用一个换根的树形DP处理出每棵树内点 i 到树内每 ...
- [jzoj 5776]【NOIP2008模拟】小x游世界树 (树形dp)
传送门 Description 小x得到了一个(不可靠的)小道消息,传说中的神岛阿瓦隆在格陵兰海的某处,据说那里埋藏着亚瑟王的宝藏,这引起了小x的好奇,但当他想前往阿瓦隆时发现那里只有圣诞节时才能到达 ...
随机推荐
- 创建一个py文件并运行
在 Linux 中,可以直接用vim 或者 vi 来编辑一个 python 文件 vim hello.py 进入编辑页面 #coding:utf-8 print("你好") (因为 ...
- 一个Accecc_Token生成和缓存和读取类,微信/小程序开发必须学
Access_Token是调用微信和小程序各种接口的临时凭证,有效期2小时(7200秒),很多接口都需要调用access_token接口生成一个access_token的,例如微信支付,微信分享,公众 ...
- RSA使用
RSA使用 今天在跟同事一起调试TCP通讯的时候,在RSA私钥解密这块,着实让我费了一番心思. 流程大致是这样的,终端登录的时候使用固定的des密码加密数据发送,平台接收后确认登录信息后,会返回一个字 ...
- 剑指offer-递归和循环-python
-斐波那契数列- 大家都知道斐波那契数列(1.1.2.3.5.8.13.21.34.……),现在要求输入一个整数n,请你输出斐波那契数列的第n项(从0开始,第0项为0). n<=39 斐波那契数 ...
- Oracle及SQLPLUS使用笔记
Oracle及SQLPLUS使用笔记 自己之前粗粗的学过MySQL,学校用的是Oracle,学生使用sqlplus,这是个命令行界面的数据库管理软件(为了学习嘛,不用图形化可以理解),这里记录一些使用 ...
- 22、nlpir 人工智能
练习介绍 [程序功能] 我们将完成一个和语义识别相关的爬虫程序,输入任意词汇.句子.文章或段落,会返回联想的词汇. [背景信息] 有一个非常牛的处理语言的网站nlpir,上面有非常多的处理语言的功能( ...
- CF528E Triangles 3000
cf luogu 既然要求三角形面积,不如考虑三角形的面积公式.因为是三条直线,所以可以考虑利用三个交点来算面积,如果这个三角形按照逆时针方向有\(ABC\)三点,那么他的面积为\(\frac{\ve ...
- PHP hex文件及bin文件读取
背景:做物联网时经常会有软件上传这种操作,上传的软件包文件常见的是hex和bin这两种. 一 hex文件读取 1 首先我们需要了解hex文件内容格式 (图及下面说明来自网络,侵权必删) :(冒号)每个 ...
- linux下重启tomcat命令
在Linux系统下,重启Tomcat使用命令操作的! 首先,进入Tomcat下的bin目录 cd /usr/local/tomcat/bin 使用Tomcat关闭命令 ./shutdown.sh 查看 ...
- 从零开始学MySQL(二)
鉴于上节篇幅以安装为主,因此对于调用mysql所需要使用的“命令”只是略微提及.随之而来就会带给读者诸多不解了,因为你会思考,这串长长的字符到底有什么特殊的含义呢?聪明的你可能早就抱着好奇心去“摆渡” ...