【LG2495】[SDOI2011]消耗战

题面

洛谷

题解

参考博客

题意

给你\(n\)个点的一棵树

\(m\)个询问,每个询问给出\(k\)个点

求将这\(k\)个点与\(1\)号点断掉的最小代价

其中\(n\leq250000\) \(m\geq1\) \(\Sigma k_i\leq500000\)

暴力

考虑直接暴力\(dp\)

设\(dp[i]\)表示处理完\(i\)的子树的最小代价

则\(dp[i]\)\(=\)\(min\Sigma dp[son_i],mn[i]\)其中\(mn[i]\)代表根节点到\(i\)节点路径上的最小边权

这样的话复杂度\(O(nm)\)无法通过此题

但观察数据范围

发现利用\(\Sigma k_i\)比较小的特点切入

虚树

思想

就是只将树上有用的点提取出来进行\(dp\),重构一棵树

这里指询问点和它们的\(lca\)

构建

考虑如何建一颗虚树。

首先我们要先对整棵树dfs一遍,求出他们的dfs序,然后对每个节点以dfs序为关键字从小到大排序

同时维护一个栈,表示从根到栈顶元素这条链。

设当前要加入的节点为\(p\),栈顶元素为\(x=s[top],lca\)为它们的最近公共祖先

因为我们按照\(dfs\)序遍历,所以\(p\)不可能是\(lca\)

那么现在会有两种情况

1.\(lca\)为\(x\),直接将\(p\)入栈

2.\(x\),\(p\)位于不同的子树中,则此时\(x\)所在子树已经遍历完了,我们需对其进行构造

设栈顶元素为\(x\),第二个为\(y\)

若\(dfn[y]>dfn[lca]\),可连边\(y->x\),将\(x\)出栈。

若\(dfn[y]=dfn[lca]\),\(y\)=\(lca\),连边\(lca->x\),子树构建完毕。

若\(dfn[y]<dfn[lca]\),即\(lca\)在\(x\)、\(y\)之间,连边\(lca->x\),\(x\)出栈,再将\(lca\)入栈,子树构建完毕。

然后重复这个过程就可以了,不理解的话可以自己手玩一下。。。

复杂度

大约是\(O(\Sigma k_i的)\),可能会带一些小常数

然后这题的代码贴在这里了:

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <climits>
#include <vector>
using namespace std; inline int gi() {
register int data = 0, w = 1;
register char ch = 0;
while (ch != '-' && (ch > '9' || ch < '0')) ch = getchar();
if (ch == '-') w = -1 , ch = getchar();
while (ch >= '0' && ch <= '9') data = data * 10 + (ch ^ 48), ch = getchar();
return w * data;
}
typedef long long ll;
#define int ll
const int MAX_N = 250005;
const int MAX_LOG_N = 19;
struct Graph { int to, cost, next; } e[MAX_N << 1]; int fir[MAX_N], e_cnt;
void clearGraph() { memset(fir, -1, sizeof(fir)); e_cnt = 0; }
void Add_Edge(int u, int v, int w) { e[e_cnt] = (Graph){v, w, fir[u]}; fir[u] = e_cnt++; }
int N, M, s[MAX_N], _top, dfn[MAX_N], mn[MAX_N];
namespace cpp1 {
int dep[MAX_N], top[MAX_N], fa[MAX_N], size[MAX_N], son[MAX_N], tim;
void dfs1(int x) {
dep[x] = dep[fa[x]] + 1; size[x] = 1;
for (int i = fir[x]; ~i; i = e[i].next) {
int v = e[i].to; if (v == fa[x]) continue;
fa[v] = x; mn[v] = min(mn[x], e[i].cost);
dfs1(v);
size[x] += size[v];
if (size[v] > size[son[x]]) son[x] = v;
}
}
void dfs2(int x, int tp) {
top[x] = tp; dfn[x] = ++tim;
if (son[x]) dfs2(son[x], tp);
for (int i = fir[x]; ~i; i = e[i].next) {
int v = e[i].to;
if (v == fa[x] || v == son[x]) continue;
dfs2(v, v);
}
}
int LCA(int x, int y) {
while (top[x] != top[y]) {
if (dep[top[x]] < dep[top[y]]) swap(x, y);
x = fa[top[x]];
}
return dep[x] > dep[y] ? y : x;
}
}
namespace cpp2 {
vector<int> G[MAX_N];
void add(int x, int y) { G[x].push_back(y); }
void ins(int x) {
if (_top == 1) { s[++_top] = x; return ; }
int lca = cpp1::LCA(x, s[_top]);
if (lca == s[_top]) return ;
while (_top > 1 && dfn[s[_top - 1]] >= dfn[lca]) add(s[_top - 1], s[_top]), _top--;
if (lca != s[_top]) add(lca, s[_top]), s[_top] = lca;
s[++_top] = x;
}
ll dfs(int x) {
if (G[x].size() == 0) return mn[x];
ll res = 0;
for (int i = 0; i < (int)G[x].size(); i++) res += dfs(G[x][i]);
G[x].clear();
return min(res, 1ll * mn[x]);
}
bool cmp(int a, int b) { return dfn[a] < dfn[b]; }
}
int a[MAX_N];
#undef int
int main () {
#define int ll
clearGraph();
N = gi();
for (int i = 1; i < N; i++) {
int u = gi(), v = gi(), w = gi();
Add_Edge(u, v, w);
Add_Edge(v, u, w);
}
fill(&mn[1], &mn[N + 1], LLONG_MAX / 2);
cpp1::dfs1(1); cpp1::dfs2(1, 1);
M = gi();
while (M--) {
int K = gi(); for (int i = 1; i <= K; i++) a[i] = gi();
sort(&a[1], &a[K + 1], cpp2::cmp);
s[_top = 1] = 1;
for (int i = 1; i <= K; i++) cpp2::ins(a[i]);
while (_top > 0) cpp2::add(s[_top - 1], s[_top]), _top--;
printf("%lld\n", cpp2::dfs(1));
}
return 0;
}

【LG2495】[SDOI2011]消耗战的更多相关文章

  1. BZOJ 2286: [Sdoi2011]消耗战

    2286: [Sdoi2011消耗战 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 2082  Solved: 736[Submit][Status] ...

  2. bzoj 2286: [Sdoi2011]消耗战 虚树+树dp

    2286: [Sdoi2011]消耗战 Time Limit: 20 Sec  Memory Limit: 512 MB[Submit][Status][Discuss] Description 在一 ...

  3. bzoj千题计划254:bzoj2286: [Sdoi2011]消耗战

    http://www.lydsy.com/JudgeOnline/problem.php?id=2286 虚树上树形DP #include<cmath> #include<cstdi ...

  4. AC日记——[SDOI2011]消耗战 洛谷 P2495

    [SDOI2011]消耗战 思路: 建虚树走树形dp: 代码: #include <bits/stdc++.h> using namespace std; #define INF 1e17 ...

  5. [BZOJ2286][SDOI2011]消耗战(虚树DP)

    2286: [Sdoi2011]消耗战 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 4998  Solved: 1867[Submit][Statu ...

  6. BZOJ2286 [Sdoi2011]消耗战 【虚树 + 树形Dp】

    2286: [Sdoi2011]消耗战 Time Limit: 20 Sec  Memory Limit: 512 MB Submit: 4261  Solved: 1552 [Submit][Sta ...

  7. 【BZOJ2286】[Sdoi2011]消耗战 虚树

    [BZOJ2286][Sdoi2011]消耗战 Description 在一场战争中,战场由n个岛屿和n-1个桥梁组成,保证每两个岛屿间有且仅有一条路径可达.现在,我军已经侦查到敌军的总部在编号为1的 ...

  8. 2286: [Sdoi2011]消耗战

    2286: [Sdoi2011]消耗战 链接 分析 虚树练习题. 构建虚树,在虚树上DP. 跟着gxb学虚-tree... 代码 #include <cstdio> #include &l ...

  9. BZOJ2286 [Sdoi2011]消耗战 和 BZOJ3611 [Heoi2014]大工程

    2286: [Sdoi2011]消耗战 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 6371  Solved: 2496[Submit][Statu ...

随机推荐

  1. C++ 全局变量不明确与 using namespace std 冲突

    写了个汉诺塔,使用全局变量count来记录步数,结果Error:count不明确 #include <iostream> using namespace std; ; void hanoi ...

  2. Anaconda安装与常用命令及方法(深度学习入门1)

    Anaconda是一个软件发行版,它附带了 conda.Python 和 150 多个科学包及其依赖项. 安装Anaconda Anaconda分为Linux.Windows.Mac等版本,去 htt ...

  3. python-二维数组实现90度旋转

    本篇主要介绍了对一个N*N的数组,如果进行90度的旋转 首先,定义一个一维数组很简单,如下: a = [i for i in range(10)] print(a) -----结果----- 0, 1 ...

  4. HDU 2088 Box of Bricks(脑洞)

    传送门: http://acm.hdu.edu.cn/showproblem.php?pid=2088 Box of Bricks Time Limit: 1000/1000 MS (Java/Oth ...

  5. Many-to-many relationships in EF Core 2.0 – Part 1: The basics

    转载这个系列的文章,主要是因为EF Core 2.0在映射数据库的多对多关系时,并不像老的EntityFramework那样有原生的方法进行支持,希望微软在以后EF Core的版本中加入原生支持多对多 ...

  6. HTML的块状、内联、内联块状元素的特点

    元素分类及特点: 1.块级元素: 在html中<div>. <p>.<h1>.<form>.<ul> 和 <li>就是块级元素. ...

  7. linux设置容器(中间件)开机自启

    /etc/rc.d/rc.local   JAVA_HOME=/usr/java/jdk1.6.0_45 su - goldsign -c '/home/goldsign/Oracle/Middlew ...

  8. kali linux 安装 Mysql Can't read from messagefile 报错解决方案

    1.下载安装包 下载地点:https://dev.mysql.com/downloads/mysql/ 或者 wget http://dev.mysql.com/get/Downloads/MySQL ...

  9. Git 学习笔记–基本操作

    Git 与 SVN 不同,是分布式的版本控制系统,不需要主服务器即可工作,实际中为了方便各个工作者间同步工作,通常还是会设置主服务器. Git的设置及初始化: 设置全局用户信息: luojiahu@u ...

  10. 使用docker搭建laravel记叙

    第一步,先从dockerhub上pull一个docker镜 docker pull laraedit/laraedit 这个docker镜像已经安装了 nginx.laravel和mysql,所以不需 ...