luogu 4381 [IOI2008]Island 单调队列 + 基环树直径 + tarjan
Description
Input
Output
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define maxn 1000006
#define setIO(s) freopen(s".in","r",stdin)
int n,edges,scc,num,ncc;
int hd[maxn],to[maxn<<1],nex[maxn<<1];
long long val[maxn<<1];
int low[maxn],dfn[maxn],vis[maxn],st[maxn],siz[maxn],idx[maxn];
ll d1[maxn],d2[maxn];
long long sum[maxn<<1],dis[maxn<<1],dep[maxn<<1];
long long now=0;
bool vis_edge[maxn<<1], vis_edge2[maxn<<1],mark[maxn];
stack<int>S;
deque<int>Q;
struct Node
{
int cur;
long long val;
Node(int cur=0,ll val=0):cur(cur),val(val){}
}nd[maxn<<1];
void addedge(int u,int v,int c)
{
nex[++edges]=hd[u],hd[u]=edges,to[edges]=v,val[edges]=1ll*c;
}
void re()
{
for(int i=0;i<=num;++i) siz[i]=0;
scc=num=ncc=0;
}
int get()
{
return Q.size() ? Q.front() : 0;
}
void tarjan(int u)
{
low[u]=dfn[u]=++scc;
S.push(u);
vis[u]=1;
for(int i=hd[u];i;i=nex[i])
{
if(vis_edge[i]) continue;
vis_edge[i]=vis_edge[i^1]=1;
int v=to[i];
if(!vis[v]) tarjan(v), low[u]=min(low[u],low[v]);
else if(vis[v]==1) low[u]=min(low[u], dfn[v]);
}
if(low[u]==dfn[u])
{
++num;
st[num]=u;
siz[num]=0;
for(;;)
{
int x=S.top();S.pop();
++siz[num];
vis[x]=-1;
idx[x]=num;
if(x==u) break;
}
}
}
void get_arr(int u)
{
for(int i=hd[u];i;i=nex[i])
{
int v=to[i];
if(vis_edge2[i] || idx[v]!=idx[u]) continue;
vis_edge2[i]=vis_edge2[i^1]=1;
if(idx[v] == idx[u])
{
nd[++ncc]=Node(u, val[i]);
get_arr(v);
break;
}
}
}
ll tmp;
void dfs(int u)
{
mark[u]=1;
d1[u]=d2[u]=0;
for(int i=hd[u];i;i=nex[i])
{
int v=to[i];
if(mark[v]) continue;
dep[v]=dep[u]+1ll*val[i];
dfs(v);
if(d1[v] + val[i] >= d1[u])
{
d2[u]=d1[u], d1[u]=d1[v] + val[i];
}
else if(d1[v] + val[i] > d2[u])
{
d2[u]=d1[v] + val[i];
}
}
now=max(now, dep[u]);
tmp=max(tmp, d1[u]+d2[u]);
}
long long solve(int root)
{
tarjan(root);
for(int i=1;i<=num;++i)
{
if(siz[i] > 1)
{
get_arr(st[i]);
break;
}
}
int k=ncc;
for(int i=1;i<=k;++i) nd[++ncc]=nd[i];
for(int i=1;i<=k;++i) mark[nd[i].cur]=1;
long long ans=0;
for(int i=1;i<=ncc;++i)
{
tmp=0;
if(i>k) dis[i]=dis[i-k];
else dep[nd[i].cur]=now=0, dfs(nd[i].cur), dis[i]=now;
ans=max(ans, tmp);
sum[i+1]=sum[i]+nd[i].val;
while(!Q.empty() && i - Q.front() >= k) Q.pop_front();
ans=max(ans,sum[i]+dis[i]+dis[get()]-sum[get()]);
while(!Q.empty() && dis[i]-sum[i]>=dis[Q.back()]-sum[Q.back()]) Q.pop_back();
Q.push_back(i);
}
sum[1]=0;
while(!Q.empty()) Q.pop_front();
re();
return ans;
}
int main()
{
// setIO("input");
scanf("%d",&n);
edges=-1;
for(int i=1;i<=n;++i)
{
int x,y;
scanf("%d%d",&x,&y);
addedge(i,x,y), addedge(x,i,y);
}
long long ou=0;
for(int i=1;i<=n;++i)
{
if(!vis[i]) { ou+=solve(i); }
}
printf("%lld\n",ou);
return 0;
}
luogu 4381 [IOI2008]Island 单调队列 + 基环树直径 + tarjan的更多相关文章
- bzoj 1791: [Ioi2008]Island 岛屿【基环树+单调队列优化dp】
我太菜了居然调了一上午-- 这个题就是要求基环树森林的基环树直径和 大概步骤就是找环->dp找每个环点最远能到达距离作为点权->复制一倍环,单调队列dp 找环是可以拓扑的,但是利用性质有更 ...
- bzoj1791[IOI2008]Island岛屿(基环树+DP)
题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=1791 题目大意:给你一棵n条边的基环树森林,要你求出所有基环树/树的直径之和.n< ...
- 【题解】Luogu P4381 [IOI2008]Island
原题传送门 题意:求基环树森林的直径(所有基环树直径之和) 首先,我们要对环上所有点的子树求出它们的直径和最大深度.然后,我们只用考虑在环上至少经过一条边的路径.那么,这种路径在环上一定有起始点和终点 ...
- 『Island 基环树直径』
Island(IOI 2008) Description 你准备浏览一个公园,该公园由 N 个岛屿组成,当地管理部门从每个岛屿 i 出发向另外一个岛屿建了一座长度为 L_i 的桥,不过桥是可以双向行走 ...
- BZOJ1791或洛谷4381 [IOI2008]Island
一道基环树的直径 BZOJ原题链接 洛谷原题链接 又是一道实现贼麻烦的题.. 显然公园其实是基环树森林,求的最长距离其实就是求每一棵基环树的直径的总和. 对于每棵基环树,其直径要么经过环,要么是某个环 ...
- [洛谷P1886]滑动窗口 (单调队列)(线段树)
---恢复内容开始--- 这是很好的一道题 题目描述: 现在有一堆数字共N个数字(N<=10^6),以及一个大小为k的窗口. 现在这个从左边开始向右滑动,每次滑动一个单位,求出每次滑动后窗口中的 ...
- [Luogu P1886]滑动窗口--单调队列入门
题目描述 现在有一堆数字共N个数字(N<=10^6),以及一个大小为k的窗口.现在这个从左边开始向右滑动,每次滑动一个单位,求出每次滑动后窗口中的最大值和最小值. 例如: The array i ...
- POJ 2823 Sliding Window(单调队列 || 线段树)题解
题意:求每个长度为k的数组的最大值和最小值 思路: 1.用线段树创建维护最大值和最小值,遍历询问,简单复习了一下...有点手生 2.单调队列: 可以看一下详解 单调队列顾名思义就是一个单调递增或者递减 ...
- [luoguP1440] 求m区间内的最小值(单调队列 || 线段树)
传送门 这种水题没必要搞线段树了,单调队列就行啊. ——代码 #include <cstdio> ; , t = ; int a[MAXN], q[MAXN]; int main() { ...
随机推荐
- 2019暑假第三周(HDFS和HBase)
Hadoop的核心是HDFS和MapReduce. 1.分布式文件系统HDFS理论方面的认知学习. 2.HDFS编程实践. 3.分布式数据库HBase.
- Matlab——程序设计
M文件 我们之前所做的运算————> 算式不太长,或想以交谈式方式进行运算 如果算式很长或是需要一再执行的算式————> 采用M文件的方式 [将指令及算式写成巨集程式然后储存成一个特别的文 ...
- 2d平台怪物逻辑
2d来回巡逻 遇到坑会自动转向 可配置单次方向行走的时间,转向等待时间等 using System; using System.Collections; using System.Collection ...
- iptables规则
iptables命令是Linux上常用的防火墙软件,是netfilter项目的一部分 iptables文件设置路径:命令:vim /etc/sysconfig/iptables-config 0x02 ...
- 【Qt开发】状态栏设置
1.在Qt 里面,状态栏显示的信息有三种类型:临时信息.一般信息和永久信息. 其中,临时信息指临时显示的信息,比如QAction 的提示等,也可以设置自己的 临时信息,比如程序启动之后显示Read ...
- Golang中log与fmt区别
关于使用log与使用fmt的区别 最初的就是直接打印出来,之后一点点升级,比如加上输出的时间,加上goroutine之间的并发操作(打印信息并不能一定按照你规定好的顺序输出来 每次输出的顺序可能会不同 ...
- python+selenium切换窗口(获取句柄信息)
current_window_handle:获得当前窗口句柄: window_handles:返回所有窗口的句柄到当前会话: switch_to.window(suoyou[0]) ========= ...
- 【IO流】FileInputStream FileOutputStream BufferInputStream BufferOutputStream
FileInputStream IO流用来处理设备之间的数据传输. Java对数据的操作是通过流的方式. Java用于操作流的类都在IO包中. 流按流向分为两种:输入流,输出流. 流按操作类型分为两种 ...
- [pwnable.kr] - wtf
Q: I don't understand why my exploit is not working. I need your help. download : http://pwnable.kr/ ...
- 【转】MySQL my.cnf 参数 详解
[client] port = 3306 //客户端端口号为3306 socket = /tmp/mysql.sock //套接字文件(sockets),这种文件一般用在网络上的资料套 ...