学习了一下圆方树。

圆方树是一种可以处理仙人掌的数据结构,具体见这里:http://immortalco.blog.uoj.ac/blog/1955

简单来讲它是这么做的:用tarjan找环,然后对每个环建立一个新点,然后将环上的边删去,并环上的每个点都连到新点上。这样我们就可以把一个环缩成一个菊花图,重复这么做,一棵仙人掌就变成一棵树啦!这棵树就叫做圆方树,其中原点叫圆点,新点叫方点。

圆方树和原仙人掌很相似,而且它又是一棵树,于是我们就可以在上面dp啦!不过要注意的是对于方点的处理,不能直接更新,要作一个单调队列,因为环上有两种走法么。

期间调了很久,而且我还不会生成数据,orz cbh。

tarjan求DCC都不会写了,我好弱啊~

#include <bits/stdc++.h>
#define N 110000
using namespace std;
int n, m, nn;
vector <int> V[N], W[N];
int low[N], dfn[N], dfsnum;
int f[N];
int ans;
int shed[N], top;
deque <int> Q;
void tarjan(int t, int f)
{
dfn[t] = low[t] = ++ dfsnum;
shed[++ top] = t;
for (int p = ; p < V[t].size(); ++ p)
if (V[t][p] != f) if (!dfn[V[t][p]])
{
tarjan(V[t][p], t);
low[t] = min(low[t], low[V[t][p]]);
if (low[V[t][p]] >= dfn[t])
{
if (shed[top] != V[t][p])
{
nn ++;
int cur;
do
{
//W[shed[top]].push_back(nn);
W[nn].push_back(cur = shed[top]);
shed[top --] = ;
}
while (cur != V[t][p]);
W[t].push_back(nn);
//W[nn].push_back(t);
}
else
{
W[t].push_back(shed[top]);
//W[shed[top]].push_back(t);
shed[top --] = ;
}
}
}
else low[t] = min(low[t], dfn[V[t][p]]);
}
void dfs(int t)
{
f[t] = ;
int s = W[t].size();
if (t > n)
{
for (int p = ; p < s; ++ p)
dfs(W[t][p]);
for (int q = ; q < (s + ) / ; ++ q)
{
while (!Q.empty() && f[W[t][Q.back()]] + Q.back() < f[W[t][q]] + q) Q.pop_back();
Q.push_back(q);
}
for (int p = , q = (s + ) / ; p < s; ++ p, q = (q == s? : q + ) )
{
if (q != s)
{
while (!Q.empty() && f[W[t][Q.back()]] + (Q.back() < p? Q.back() + s + - p: Q.back() - p) < f[W[t][q]] + (q < p? q + s + - p: q - p)) Q.pop_back();
Q.push_back(q);
}
if (Q.front() == p) Q.pop_front();
if (!Q.empty()) ans = max(ans, f[W[t][p]] + f[W[t][Q.front()]] + (Q.front() > p? Q.front() - p: Q.front() + s + - p));
}
for (int p = ; p < (s + ) / ; ++ p) f[t] = max(f[t], f[W[t][p]] + p);
for (int p = (s + ) / ; p < s; ++ p) f[t] = max(f[t], f[W[t][p]] + s - - p);
Q.clear();
}
else
{
int mx1 = , mx2 = ;
for (int p = ; p < s; ++ p)
{
dfs(W[t][p]);
if (f[W[t][p]] + > mx1) mx2 = mx1, mx1 = f[W[t][p]] + ;
else if (f[W[t][p]] + > mx2) mx2 = f[W[t][p]] + ;
}
ans = max(ans, mx1 + mx2);
f[t] = mx1;
}
}
int main()
{
scanf("%d%d", &n, &m);
for (int i = ; i <= m; ++ i)
{
int k, a;
scanf("%d%d", &k, &a);
for (int i = ; i <= k; ++ i)
{
int b;
scanf("%d", &b);
V[a].push_back(b);
V[b].push_back(a);
a = b;
}
}
nn = n;
tarjan(, );
//puts("haha");
dfs();
printf("%d\n", ans);
}

bzoj1023: [SHOI2008]cactus仙人掌图的更多相关文章

  1. bzoj千题计划113:bzoj1023: [SHOI2008]cactus仙人掌图

    http://www.lydsy.com/JudgeOnline/problem.php?id=1023 dp[x] 表示以x为端点的最长链 子节点与x不在同一个环上,那就是两条最长半链长度 子节点与 ...

  2. BZOJ1023:[SHOI2008]cactus仙人掌图(圆方树,DP,单调队列)

    Description 如果某个无向连通图的任意一条边至多只出现在一条简单回路(simple cycle)里,我们就称这张图为仙人掌图(cactus). 所谓简单回路就是指在图上不重复经过任何一个顶点 ...

  3. BZOJ1023: [SHOI2008]cactus仙人掌图(仙人掌dp)

    Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 3467  Solved: 1438[Submit][Status][Discuss] Descripti ...

  4. BZOJ1023[SHOI2008]cactus仙人掌图 【仙人掌DP】

    题目 如果某个无向连通图的任意一条边至多只出现在一条简单回路(simple cycle)里,我们就称这张图为仙人掌 图(cactus).所谓简单回路就是指在图上不重复经过任何一个顶点的回路. 举例来说 ...

  5. BZOJ1023: [SHOI2008]cactus仙人掌图(仙人掌)

    Description 如果某个无向连通图的任意一条边至多只出现在一条简单回路(simple cycle)里,我们就称这张图为仙人掌图(cactus).所谓简单回路就是指在图上不重复经过任何一个顶点的 ...

  6. [bzoj1023][SHOI2008]cactus 仙人掌图 (动态规划)

    Description 如果某个无向连通图的任意一条边至多只出现在一条简单回路(simple cycle)里,我们就称这张图为仙人图(cactus).所谓简单回路就是指在图上不重复经过任何一个顶点的回 ...

  7. 2018.10.29 bzoj1023: [SHOI2008]cactus仙人掌图(仙人掌+单调队列优化dp)

    传送门 求仙人掌的直径. 感觉不是很难. 分点在环上面和不在环上分类讨论. 不在环上直接树形dpdpdp. 然后如果在环上讨论一波. 首先对环的祖先有贡献的只有环上dfsdfsdfs序最小的点. 对答 ...

  8. bzoj1023 [SHOI2008]cactus仙人掌图 & poj3567 Cactus Reloaded——求仙人掌直径

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=1023    http://poj.org/problem?id=3567 仙人掌!直接模仿 ...

  9. bzoj千题计划224:bzoj1023: [SHOI2008]cactus仙人掌图

    又写了一遍,发出来做个记录 #include<cstdio> #include<algorithm> #include<iostream> using namesp ...

随机推荐

  1. 嵌入式Linux驱动学习之路(十一)按键驱动-中断机制

    轮询方式: 和led驱动不同的是在配置IO引脚的时候,把LED的输出引脚换成输入,在read函数中向外发送io的状态.必须由应用程序不断的来查询当前IO口的状态来判断. 中断方式: /* file n ...

  2. [bzoj1597][usaco2008 mar]土地购买 (动态规划+斜率优化)

    Description 农夫John准备扩大他的农场,他正在考虑N (1 <= N <= 50,000) 块长方形的土地. 每块土地的长宽满足(1 <= 宽 <= 1,000, ...

  3. Javascript函数式编程要掌握的知识点讲解

    一:理解call和apply 及arguments.callee ECMAScript3给Function的原型定义了两个方法,他们是Function.prototype.call 和 Functio ...

  4. jquery的show/hide/toggle详解

    通过阅读源码我们发现show,hide,toggle调用了showHide和isHidden这2个方法,所以我们要搞明白原理必须先看一下这2个方法. jQuery.fn.extend({ ...... ...

  5. 快速部署tomcat项目的Shell脚本

    为了做集群测试,在每台机器上装了3个tomcat,每次发布项目的时候都要反复敲一些命令,重启tomcat之前先检查tomcat进程有没有停掉,没有还要手动kill该进程. 发布次数多了,操作就比较繁琐 ...

  6. vc++2013中使用MySQL connector/C++ 1.1.4静态链接报错

    包含头文件 #include <mysql_connection.h> #include <mysql_driver.h> #include <cppconn/state ...

  7. 2.1 python使用MongoDB 示例代码

    import pymongo client = pymongo.MongoClient('localhost', 27017) # MongoDB 客户端 walden = client['walde ...

  8. Java调用Linux命令

    // int tp = 1 返回执行结果  非1 返回命令执行后的输出 public static String runCommand(String cmd, int tp) { StringBuff ...

  9. 最为简易的yii 教程(一)

    了解目录的框架结构 framework主要有 base          框架核心组件 caching        缓存组件 db                数据库组件 gii          ...

  10. jeecg单步调试

    自己没本事写开发平台,用别人的又各种担心,想学着别人弄个单步调试,老是出现这个"source not found"-- 百度各种方法都搞不定,担心是自己安装错了?这个jeecg本身 ...