BZOJ1791: [Ioi2008]Island 岛屿
Description
Input
Output
Sample Input
3 8
7 2
4 2
1 4
1 9
3 4
2 3
Sample Output
HINT
题解Here!
- 位于某个子树中且不在环上。即:不经过环上的点。
- 部分经过环,且去掉处于环上的边后位于两颗不同的子树内。即:至少经过环上两个点。
所以我们得把环找出来。
同时求出:从每个环上节点出发在不经过环的前提下的最长链。
这里当然借鉴一下网上大佬得拓扑排序找环辣!
于是这个玩意直接丢给拓扑排序就好。
用求树的直径的方法更新答案处理第一种情况。
对于第二种情况,就等价于从环上找出两点$i,j$使得$dp[i]+dis(i,j)+dp[j]$最大。
$dis(i,j)$为两点在环上的最长距离。
$dp[i]$表示以$i$为根的子树内以一个端点为根的最长链。
记得考虑顺时针和逆时针两种走法。
#include<iostream>
#include<algorithm>
#include<cstdio>
#define MAXN 1000010
using namespace std;
int n,c=1,T=0;
int head[MAXN],degree[MAXN],colour[MAXN],que[MAXN<<1];
long long ans=0,dis[MAXN],dp[MAXN],f[MAXN<<1],g[MAXN<<1];
bool vis[MAXN];
struct Edge{
int next,to,w;
}a[MAXN<<1];
inline int read(){
int date=0,w=1;char c=0;
while(c<'0'||c>'9'){if(c=='-')w=-1;c=getchar();}
while(c>='0'&&c<='9'){date=date*10+c-'0';c=getchar();}
return date*w;
}
inline void add(int u,int v,int w){
a[c].to=v;a[c].w=w;a[c].next=head[u];head[u]=c++;
a[c].to=u;a[c].w=w;a[c].next=head[v];head[v]=c++;
degree[u]++;degree[v]++;
}
void bfs(int rt){
int l=1,r=1,u,v;
que[1]=rt;
colour[rt]=T;
while(l<=r){
u=que[l];
for(int i=head[u];i;i=a[i].next){
v=a[i].to;
if(!colour[v]){
que[++r]=v;
colour[v]=T;
}
}
l++;
}
}
void topsort(){
int l=1,r=0,u,v;
for(int i=1;i<=n;i++)if(degree[i]==1)que[++r]=i;
while(l<=r){
u=que[l];
for(int i=head[u];i;i=a[i].next){
v=a[i].to;
if(degree[v]>1){
dis[colour[u]]=max(dis[colour[u]],dp[u]+dp[v]+a[i].w);
dp[v]=max(dp[v],dp[u]+a[i].w);
degree[v]--;
if(degree[v]==1)que[++r]=v;
}
}
l++;
}
}
void solve(int rt,int colour){
int u,v,m=0;
u=rt;
for(int i=1;i;){
f[++m]=dp[u];
degree[u]=1;
for(i=head[u];i;i=a[i].next){
v=a[i].to;
if(degree[v]>1){
g[m+1]=g[m]+a[i].w;
u=v;
break;
}
}
}
if(m==2){
int len=0;
for(int i=head[u];i;i=a[i].next)if(a[i].to==rt)len=max(len,a[i].w);
dis[colour]=max(dis[colour],dp[rt]+dp[u]+len);
}
else{
int l=1,r=1;
que[1]=1;
for(int i=head[u];i;i=a[i].next)if(a[i].to==rt){
g[m+1]=g[m]+a[i].w;
break;
}
for(int i=1;i<m;i++){
f[m+i]=f[i];
g[m+i]=g[m+1]+g[i];
}
for(int i=2;i<2*m;i++){
while(l<=r&&i-que[l]>=m)l++;
dis[colour]=max(dis[colour],f[i]+f[que[l]]+g[i]-g[que[l]]);
while(l<=r&&f[que[r]]-g[que[r]]<=f[i]-g[i])r--;
que[++r]=i;
}
}
}
void work(){
for(int i=1;i<=n;i++)if(degree[i]>1&&!vis[colour[i]]){
vis[colour[i]]=true;
solve(i,colour[i]);
ans+=dis[colour[i]];
}
printf("%lld\n",ans);
}
void init(){
int x,w;
n=read();
for(int i=1;i<=n;i++){
x=read();w=read();
add(i,x,w);
}
for(int i=1;i<=n;i++)if(!colour[i]){T++;bfs(i);}
topsort();
}
int main(){
init();
work();
return 0;
}
BZOJ1791: [Ioi2008]Island 岛屿的更多相关文章
- [bzoj1791][ioi2008]Island 岛屿(基环树、树的直径)
[bzoj1791][ioi2008]Island 岛屿(基环树.树的直径) bzoj luogu 题意可能会很绕 一句话:基环树的直径. 求直径: 对于环上每一个点记录其向它的子树最长路径为$dp_ ...
- bzoj1791: [Ioi2008]Island 岛屿 单调队列优化dp
1791: [Ioi2008]Island 岛屿 Time Limit: 20 Sec Memory Limit: 162 MBSubmit: 1826 Solved: 405[Submit][S ...
- bzoj千题计划114:bzoj1791: [Ioi2008]Island 岛屿
http://www.lydsy.com/JudgeOnline/problem.php?id=1791 就是求所有基环树的直径之和 加手工栈 #include<cstdio> #incl ...
- BZOJ1791[Ioi2008]Island 岛屿 ——基环森林直径和+单调队列优化DP+树形DP
题目描述 你将要游览一个有N个岛屿的公园.从每一个岛i出发,只建造一座桥.桥的长度以Li表示.公园内总共有N座桥.尽管每座桥由一个岛连到另一个岛,但每座桥均可以双向行走.同时,每一对这样的岛屿,都有一 ...
- [BZOJ1791][IOI2008]Island岛屿(环套树DP)
同NOI2013快餐店(NOI出原题?),下面代码由于BZOJ栈空间过小会RE. 大致是对每个连通块找到环,在所有内向树做一遍DP,再在环上做两遍前缀和优化的DP. #include<cstdi ...
- bzoj1791[IOI2008]Island岛屿(基环树+DP)
题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=1791 题目大意:给你一棵n条边的基环树森林,要你求出所有基环树/树的直径之和.n< ...
- BZOJ1791 [Ioi2008]Island 岛屿[基环树+单调队列优化DP]
基环树直径裸题. 首先基环树直径只可能有两种形式:每棵基环树中的环上挂着的树的直径,或者是挂在环上的两个树的最大深度根之间的距离之和. 所以,先对每个连通块跑一遍,把环上的点找出来,然后对环上每个点跑 ...
- 【BZOJ 1791】 [Ioi2008]Island 岛屿
Description 你将要游览一个有N个岛屿的公园.从每一个岛i出发,只建造一座桥.桥的长度以Li表示.公园内总共有N座桥.尽管每座桥由一个岛连到另一个岛,但每座桥均可以双向行走.同时,每一对这样 ...
- bzoj 1791: [Ioi2008]Island 岛屿
#include<iostream> #include<cstdio> #define M 1000009 using namespace std; *M],cnt,n,hea ...
随机推荐
- 6. datasource - mysql【从零开始学Spring Boot】
在任何一个平台都逃离不了数据库的操作,那么在spring boot中怎么接入数据库呢? 很简单,我们需要在application.properties进行配置一下,application.proper ...
- 2016.3.16__HTML5新特性__第八天
HTML 5 + CSS 3 假设您认为这篇文章还不错,能够去H5专题介绍中查看很多其它相关文章. 今日代码非常冗杂,所以非常多内容直接摘自网上,假设造成您的不适.请留言告知. 非常感谢. 输入标签, ...
- Unity里面的自动寻路(二)
接着我的 上一篇自动寻路文章,这一次我们就来学习一下与自动寻路有关的组件吧.Unity中与自动寻路相关的组件主要有两个:NavMeshAgent ( 又称导航网格代理 ),Off Mesh Link ...
- C++Singleton的DCLP(双重锁)实现以及性能测评
本文系原创,转载请注明:http://www.cnblogs.com/inevermore/p/4014577.html 根据维基百科,对单例模式的描述是: 确保一个类只有一个实例,并提供对该 ...
- Python——os.path.dirname(__file__) 与 os.path.join(str,str)
Python os.path.dirname(__file__) Python os.path.join(str,str) (1).当"print os.path.dirname(__f ...
- STL学习笔记(算法概述)
算法头文件 要运用C++标准程序库的算法,首先必须包含头文件<algorithm> 使用STL算法时,经常需要用到仿函数以及函数配接器.它们定义域<functional>头文件 ...
- .aspx 页面引用命名空间
一.单个页面引用: <%@ Import Namespace="" %> 二.所有页面引用,Web.config配置如下: <system.web> < ...
- 公网通过代理访问阿里云vpc redis
前提条件 如果您需要从本地 PC 端访问 Redis 实例进行数据操作,可以通过在 ECS 上配置端口映射或者端口转发实现.但必须符合以下前提条件: 若 Redis 实例属于专有网络(VPC),ECS ...
- Android使用Fragment打造万能页面切换框架
首先我们来回顾一下传统用Activity进行的页面切换.activity之间切换.首先须要新建intent对象,给该对象设置一些必须的參数,然后调用startActivity方法进行页面跳转. 假设须 ...
- C++中没有定义类的引用。
在有时候由于类太大.须要在类在后面定义: 比如: class Y{ void f(X); }; class X{ //一些成员数据和函数 }; //error 由于c++要求不论什么一个变量在引用之前 ...