虚树dp

虚树的主要思想:

  • 不遍历没用的的节点以及没用的子树,从而使复杂度降低到\(\sum\limits k\)(k为询问的节点的总数)。

所以怎么办:

  • 只把询问节点和其LCA放入询问的数组中。

1、建虚树

q.clear();
int m;
scanf("%d",&m);
for(int i=1;i<=m;++i){
int x;
scanf("%d",&x);
v[x]=1;
q.push_back(x);
}
sort(q.begin(),q.end(),cmp);
for(int i=0;i<m-1;++i){
q.push_back(LCA(q[i],q[i+1]));
}
q.push_back(1);
sort(q.begin(),q.end());
q.erase(unique(q.begin(),q.end()),q.end());
sort(q.begin(),q.end(),cmp);

q是一个vector,我们开始先对所有节点按欧拉序(即深度优先搜索是访问的顺序)排序,然后对每两个相邻的节点将LCA放入q中(可知这样一定会将所有有效节点放入q中)。然后一波去重,再按欧拉序排序即可。

如果你还不会LCA的话请到这里

2.遍历虚树

然后我们得到了一个遍历表,向深度优先搜索一样搜一遍即可。

这里要注意每个节点x只有当他下一个节点y是他的子节点时(即\(dfn[x]+size[x]>=dfn[y]\)时,其中\(dfn\)为欧拉序,\(size\)为子树大小)才访问下一个节点,并用下一个点的信息更新当前节点。

对于两点间的最短树边,我们可以用倍增来寻找(当然也可以用st表\(O(1)\)求,但这题并不要求)。

long long getmin(int x,int lca){
int ret=1e18;
for(int i=t-1;~i;--i){
if(dep[x]-(1<<i)>=dep[lca]){
ret=min(ret,c[x][i]);
x=fa[x][i];
}
}
return ret;
}

3.树型DP

对于每个节点,如果他一定要被割掉,则当前点的最小花费为当前点到父亲节点的最小树边,否则为所有子节点的最小花费和当前点到父亲节点的最小树边的最小值。

void dfs1(){
int x=q[it];
long long ret=0;
while(1){
if(it+1==q.size())break;
if(dfn[q[it+1]]<=dfn[x]+sz[x]-1){
int y=q[++it];
if(v[y]==1){
dfs1();
dp[y]=getmin(y,x);
}
else dp[y]=1e18,dfs1(),dp[y]=min(dp[y],getmin(y,x));
ret+=dp[y];
}else break;
}
if(ret)dp[x]=min(dp[x],ret);
}

代码中的it为当前访问到的节点在q中的编号

然后就可以写出代码了,需要注意一些初始化的细节:

#include<bits/stdc++.h>
using namespace std;
const int N=600010,t=20;
int n;
int tot,bian[N<<1],nxt[N<<1],zhi[N<<1],head[N];
void add(int x,int y,int z){
tot++,bian[tot]=y,zhi[tot]=z,nxt[tot]=head[x],head[x]=tot;
}
int dfn[N],cnt;
int fa[N][t],c[N][t],dep[N],sz[N];
vector<int>q;
bool cmp(int x,int y){
return dfn[x]<dfn[y];
}
void dfs(int x,int f){
sz[x]=1;
dfn[x]=++cnt;
dep[x]=dep[f]+1;
fa[x][0]=f;
for(int i=1;i<t;++i){
fa[x][i]=fa[fa[x][i-1]][i-1];
c[x][i]=min(c[x][i-1],c[fa[x][i-1]][i-1]);
}
for(int i=head[x];i;i=nxt[i]){
int y=bian[i];
if(y==f)continue;
c[y][0]=zhi[i];
dfs(y,x);
sz[x]+=sz[y];
}
}
int LCA(int x,int y){
if(dep[x]<dep[y])swap(x,y);
for(int i=t-1;~i;--i){
if(dep[x]-(1<<i)>=dep[y]){
x=fa[x][i];
}
}
if(x==y)return x;
for(int i=t-1;~i;--i){
if(fa[x][i]!=fa[y][i]){
x=fa[x][i],y=fa[y][i];
}
}
return fa[x][0];
}
long long getmin(int x,int lca){
int ret=1e18;
for(int i=t-1;~i;--i){
if(dep[x]-(1<<i)>=dep[lca]){
ret=min(ret,c[x][i]);
x=fa[x][i];
}
}
return ret;
}
#define IT vector<int>::iterator
long long dp[N];
int it,v[N];
void dfs1(){
int x=q[it];
long long ret=0;
while(1){
if(it+1==q.size())break;
if(dfn[q[it+1]]<=dfn[x]+sz[x]-1){
int y=q[++it];
if(v[y]==1){
dfs1();
dp[y]=getmin(y,x);
}
else dp[y]=1e18,dfs1(),dp[y]=min(dp[y],getmin(y,x));
ret+=dp[y];
}else break;
}
if(ret)dp[x]=min(dp[x],ret);
}
int main(){
dfn[0]=1e9;
cin>>n;
for(int i=1;i<n;++i){
int x,y,z;
scanf("%d%d%d",&x,&y,&z);
add(x,y,z);
add(y,x,z);
}
dfs(1,0);
int T;
cin>>T;
while(T--){
q.clear();
int m;
scanf("%d",&m);
for(int i=1;i<=m;++i){
int x;
scanf("%d",&x);
v[x]=1;
q.push_back(x);
}
sort(q.begin(),q.end(),cmp);
for(int i=0;i<m-1;++i){
q.push_back(LCA(q[i],q[i+1]));
}
q.push_back(1);
sort(q.begin(),q.end());
q.erase(unique(q.begin(),q.end()),q.end());
sort(q.begin(),q.end(),cmp);
m=q.size();it=0;
dp[q[0]]=1e18;
dfs1();
for(int i=0;i<m;++i){
v[q[i]]=0;
}
printf("%lld\n",dp[*q.begin()]);
}
}

[SDOI2011]消耗战(虚树+树形动规)的更多相关文章

  1. BZOJ.2286.[SDOI2011]消耗战(虚树 树形DP)

    题目链接 BZOJ 洛谷P2495 树形DP,对于每棵子树要么逐个删除其中要删除的边,要么直接断连向父节点的边. 如果当前点需要删除,那么直接断不需要再管子树. 复杂度O(m*n). 对于两个要删除的 ...

  2. BZOJ 2286: [Sdoi2011]消耗战 虚树 树形dp 动态规划 dfs序

    https://www.lydsy.com/JudgeOnline/problem.php?id=2286 wa了两次因为lca犯了zz错误 这道题如果不多次询问的话就是裸dp. 一棵树上多次询问,且 ...

  3. BZOJ2286: [Sdoi2011]消耗战(虚树/树形DP)

    Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 5246  Solved: 1978[Submit][Status][Discuss] Descript ...

  4. 【BZOJ-2286】消耗战 虚树 + 树形DP

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

  5. 【BZOJ2286】【SDOI2011】消耗战 [虚树][树形DP]

    消耗战 Time Limit: 20 Sec  Memory Limit: 512 MB[Submit][Status][Discuss] Description 在一场战争中,战场由n个岛屿和n-1 ...

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

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

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

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

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

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

  9. BZOJ 2286 消耗战 (虚树+树形DP)

    给出一个n节点的无向树,每条边都有一个边权,给出m个询问,每个询问询问ki个点,问切掉一些边后使得这些顶点无法与顶点1连接.最少的边权和是多少.(n<=250000,sigma(ki)<= ...

随机推荐

  1. 视频基础知识:浅谈视频会议中H.264编码标准的技术发展

    浅谈视频会议中H.264编码标准的技术发展 浅谈视频会议中H.264编码标准的技术发展 数字视频技术广泛应用于通信.计算机.广播电视等领域,带来了会议电视.可视电话及数字电视.媒体存储等一系列应用,促 ...

  2. make/makefile中的加号+,减号-和at号@的含义

    http://www.crifan.com/order_make__makefile_in_the_plus__minus_-_and_at_the_meaning_of_numbers/ 在看mak ...

  3. NOIP水题测试(2017082501)

    日常水题测试又来了! 以后答案都以单题形式公布. 下面看今天的水题: 时间限制:5小时 题目一:无法形容的水 题目二:比上一题还水 题目三:一元三次方程求解 题目四:单词接龙 题目五:统计单词个数 题 ...

  4. 指令发布中如何实现new新消息的提醒?

    设计思路:反馈后,最急需了解反馈结果的是申请人,故给每一条反馈信息添加一个查看状态的字段,如CK_STATUS,并为这个状态设计为char(1)类型,java bean中使用integer可以实现默认 ...

  5. JMeter测试工具.jmx文件详解

    摘要:了解.jmx文件格式类型,对jmeter二次开发与拓展有很大的帮助,当然也可以利用python对其进行一些处理(生成一些测试用例,对jmx文件进行 ”增删改查“). 一个完整用例的.jmx文件基 ...

  6. excel2007VBA绘图2

    '--------------------------------------------------------Module1------------------------------------ ...

  7. 未找到导入的项目“C:\Program Files\MSBuild\Microsoft\Silverlight\v5.0\Microsoft.Silverlight.CSharp.targets”。

    问题描述: 原先创建的Silverlight程序,后来系统重装了,再打开Silverlight程序时提示:C:\Users\yzg\Desktop\ocr\TJSilverlight\TJSilver ...

  8. bzoj2938(ac自动机)

    刚学了ac自动机,去hzwer上找了道练习题: 串是安全的就说明ac自动机不会找到匹配,考虑ac自动机的匹配过程: 我们把val等于1的点删掉和fail指针指向被删掉的点删掉: 如果剩下的图有环,就有 ...

  9. 学习c++需要知道

    1.c++投入和经济产出比目前较低,c++很复杂,学习难度高.周期长.而且开发效率低. 2.但是c++趋紧于完美语言,它囊括了多种编程范型.它运行效率高,对计算机有着极致的控制能力.适用于低层. 3. ...

  10. linux-CentOS初学terminal命令(3)rm、chmod、mkdir、who、w、id、systemctl、

    PS 1:windows不允许出现字母相同,但是大小写不同的文件名,因为在windows下会将它们认作是同名. 但是linux允许出现字母相同,大小写不同的文件名. ps 2:prompt 提示 1. ...