P2495 [SDOI2011]消耗战-虚树+树形dp

https://www.luogu.com.cn/problem/P2495

虚树:当我们在解决树形dp的问题的时候,题目中会给出一些询问,询问涉及的关键节点不多,并保证总的点数规模的时候,我们就可以使用虚数,如果每次询问都对整个树进行dp的话,显然会超市。要注意每一次的初始化

题目:给出n点,有n-1条带权边,有q个询问,每次给出k个点,询问要使得这k个点和1号节点不连通,需要移除的边的总权值最小是多少。

思路:

令dp[n]表示从nn开始不能到达其子树中的关键点所需切断的最小边权和.

令me[u]表示切断11到uu的路径中的边权最小值.

设vv是uu的直接儿子.

如果vv是关键节点,那么dp[u]+=me[v],否则dp[u]+=min(me[v],dp[v])
(第2个转移方程的解释:要么直接切断1−v的路径,要么使得从v出发不能到达其子树的关键点.)

显然我们不能针对每个询问对整颗子树进行dpdp,时间复杂度过高,而我们发现那些非关键点我们没有必要在dp的时候考虑,所以使用虚树.

参考:大神题解

#include <iostream>
#include <cstdio>
#include <fstream>
#include <algorithm>
#include <cmath>
#include <deque>
#include <vector>
#include <queue>
#include <string>
#include <cstring>
#include <map>
#include <stack>
#include <set>
#include <cstdlib>
#define INF 0x3f3f3f3f3f3f3f3f
#define inf 0x3f3f3f3f
#define FILL(a,b) (memset(a,b,sizeof(a)))
#define re register
#define lson rt<<1
#define rson rt<<1|1
#define lowbit(a) ((a)&-(a))
#define ios std::ios::sync_with_stdio(false);std::cin.tie(0);std::cout.tie(0);
#define fi first
#define rep(i,n) for(int i=0;(i)<(n);i++)
#define rep1(i,n) for(int i=1;(i)<=(n);i++)
#define se second
#define scd(a) scanf("%d",&a)
#define scdd(a,b) scanf("%d%d",&a,&b)
#define scddd(a,b,c) scanf("%d%d%d",&a,&b,&c)
#define ac cout<<ans<<"\n"
#define F(x) ((x)/3+((x)%3==1?0:tb))
#define G(x) ((x)<tb?(x)*3+1:((x)-tb)*3+2)
using namespace std;
typedef long long  ll;
typedef unsigned long long  ull;
typedef pair<ll,ll> pii;
int dx[4]= {-1,1,0,0},dy[4]= {0,0,1,-1};
const ll mod=1e9+7;
const ll N =250007;
const ll M =250000;
const double eps = 1e-4;
//const double pi=acos(-1);
ll qk(ll a,ll b){ll ans=1;while(b){if(b&1) ans=(ans*a)%mod;a=(a*a)%mod;b/=2;}return ans%mod;}
int n;
int cnt,m;
int ldfn[N],rdfn[N];//当前节点的dfs序
int fa[N][20];//f[i][j]表示i节点向上2^j步的祖先
int dep[N];//节点深度
int stk[N];//栈,用于构建虚树
int d[N<<1];//用于存储虚树的节点编号。由于LCA可能有重复,因此需要开2倍
int vis[N];//标记当前节点是否为关键点
vector<pii> RG[N];
vector<int> VG[N];
ll me[N];
void predfs(int u,int f,int depth)
{ldfn[u]=++cnt;dep[u]=depth;fa[u][0]=f;for(int j=1;j<20;j++){fa[u][j]=fa[fa[u][j-1]][j-1];}for(pii k:RG[u]){int v=k.fi;if(v==f) continue;me[v]=k.second;if(u!=1) me[v]=min(me[u],me[v]);predfs(v,u,depth+1);}rdfn[u]=cnt;
}int lca(int u,int v)
{if(dep[u]<dep[v])swap(u,v);int delta=dep[u]-dep[v];for(int j=0;j<20&&delta;j++){if(delta&(1<<j))u=fa[u][j];}if(u==v) return u;for(int j=20-1;j>=0;j--){if(fa[u][j]!=fa[v][j]){u=fa[u][j];v=fa[v][j];}}return fa[u][0];
}bool cmp(int x,int y)
{return ldfn[x]<ldfn[y];
}void build()
{sort(d+1,d+m+1,cmp);int keynum=m;for(int i=1;i<keynum;i++) d[++m]=lca(d[i],d[i+1]);sort(d+1,d+m+1,cmp);m=unique(d+1,d+m+1)-d-1;int top=0;stk[++top]=d[1];for(int i=2;i<=m;i++){while(top&&rdfn[stk[top]]<ldfn[d[i]])top--;if(top)VG[stk[top]].push_back(d[i]);stk[++top]=d[i];}
}
ll DP(int u){ll cost=0;for(int v:VG[u]){cost+=min(me[v],DP(v));}VG[u].clear();//初始化if(vis[u]) return me[u];else return cost;
}
void sovle(){cin>>n;for(int i=1;i<n;i++){int u,v,c;cin>>u>>v>>c;RG[u].push_back({v,c});RG[v].push_back({u,c});}predfs(1,0,0);int q;cin>>q;while(q--){cin>>m;m++;d[1]=1;for(int i=2;i<=m;i++){cin>>d[i];vis[d[i]]=1;}build();cout<<DP(1)<<endl;for(int i=2;i<=m;i++) vis[d[i]]=0;//初始化部分很关键cnt=0;}
}
int main()
{
#ifdef LOCALfreopen("in.txt", "r", stdin);
#else// iosint t=1;// cin>>t;while(t--) sovle();
#endif // LOCALreturn 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/315307.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

DDD中的聚合和UML中的聚合以及组合的关系

UML&#xff1a;聚合关系&#xff1a;成员对象是整体的一部分&#xff0c;但是成员对象可以脱离整体对象独立存在。如汽车&#xff08;Car&#xff09;与引擎&#xff08;Engine&#xff09;、轮胎&#xff08;Wheel&#xff09;、车灯&#xff08;Light&#xff09;之间的关系…

Network 黑暗爆炸 - 3732 倍增lca || Kruskal重构树

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 两点间最长边最小值一定是最小生成树上两点间的最大值&#xff0c;这个比较容易证&#xff0c;就不多说了。 知道这个结论后&#xff0c; 我们直接跑一个KruskalKruskalKruskal建树&#x…

Manacher入门

写在前面 manachermanachermanacher比想象中好理解得多 至少它给了我学习字符串的信心 能干啥 manachermanachermanacher&#xff0c;中文马拉车&#xff08;您别说&#xff0c;这名字还挺形象&#xff09;&#xff0c;主要用于计算字符串每一个位置为对称中心的回文串长度&a…

HDU.6761.Minimum Index(Lyndon分解)

题目大意&#xff1a;给一个字符串&#xff0c;求字符串的所有前缀的最小后缀&#xff1b; 思路&#xff1a;主要还是要理解Lyndon串的算法的整一个过程 参考&#xff1a;搬运来自 #include <iostream> #include <cstdio> #include <fstream> #include <…

C#规范整理·异常与自定义异常

这里会列举在C#中处理CLR异常方面的规范&#xff0c;帮助大家构建和开发一个运行良好和可靠的应用系统。前言迄今为止&#xff0c;CLR异常机制让人关注最多的一点就是“效率”问题。其实&#xff0c;这里存在认识上的误区&#xff0c;因为正常控制流程下的代码运行并不会出现问…

Codeforces Round #619 (Div. 2) D. Time to Run 矩阵回路构造

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个n∗mn*mn∗m的矩阵&#xff0c;每两个相邻点之间有个双向边&#xff0c;问你能不能走满kkk步&#xff0c;每一步走的边不同&#xff0c;但是点可以相同&#xff0c;换句话说就是走的边不能再走了。输…

HDU - 6756 Finding a MEX-分块思想

https://vjudge.net/problem/HDU-6756 题目大意&#xff1a;给你一个无向图&#xff0c;每个点有权值a&#xff0c;将f&#xff08;u&#xff09;定义为对u的邻居的集合求mex&#xff1b; 有两个操作&#xff1a; 1&#xff1a;将u的权值修改为x 2&#xff1a;查询f&#x…

中高级数论 [欧拉函数线性筛,二次剩余]

欧拉函数线性筛 对于素数ppp, φ(p∗i){p−1i1p∗φ(i)p∣i(p−1)∗φ(i)p∤i\varphi (p*i) \begin{cases} p-1& i1\\ p*\varphi(i)& p \mid i\\ (p-1)*\varphi(i) & p \nmid i \end{cases}φ(p∗i)⎩⎪⎨⎪⎧​p−1p∗φ(i)(p−1)∗φ(i)​i1p∣ip∤i​ 证明&am…

C#中await/async闲说

自从C#5.0增加异步编程之后&#xff0c;异步编程越来越简单&#xff0c;async和await用的地方越来越多&#xff0c;越来越好用&#xff0c;只要用异步的地方都是一连串的异步&#xff0c;如果想要异步编程的时候&#xff0c;需要从底层开始编写&#xff0c;这样后边使用的时候就…

Codeforces Round #619 (Div. 2) E. Nanosoft 思维 + 二维前缀和

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 考虑到最大面积是由四种颜色构成的&#xff0c;且四种颜色可以从中心扩展出去&#xff0c;所以我们分别维护四种颜色的二维前缀和&#xff0c;O(1)O(1)O(1)计算矩阵内颜色的个数。现在我们…

CF 613D - Kindom and its Cities-虚树

https://codeforces.com/problemset/problem/613/D 题目大意&#xff1a;有n个点&#xff0c;每次询问要最少移除多少个点使得关键点分隔 思路&#xff1a;首先&#xff0c;如果关键点是相连的话&#xff0c;那么一定是不行的。 然后因为询问的总节点数不多&#xff0c;可以…

FFT:从入门到沉迷

终于学会了FFTFFTFFT并无法自拔 啥是FFT FFTFFTFFT指快速傅里叶变换 在OIOIOI中的应用是O(nlogn)O(nlogn)O(nlogn)计算函数卷积 人话&#xff1a;多项式乘法 多项式毒瘤模板题的万恶之源 系数表达式与点值表达式 系数表达式就是平常的表示方法 f(x)∑i0naixif(x)\sum_{i0…

HDU - 6769-In Search of Gold-二分+树形dp

https://vjudge.net/problem/HDU-6769 题目大意&#xff1a;给你n个点&#xff0c;有n-1条边&#xff0c;每条边有a&#xff0c;b两个权值&#xff0c;给你一个k&#xff0c;恰好有k条边的权值取a&#xff0c;其余取b的时候&#xff0c;树的直径的最小值。 思路&#xff1a;答…

Codeforces Round #701 (Div. 2) E. Move and Swap 思维 + dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 由于是按层来的&#xff0c;所以我们肯定先按照层来分组。 定义dp[i]dp[i]dp[i]为红棋在位置iii的时候的最大得分和。 先考虑不换的情况&#xff0c;我们对于每个点都从他的父节点转移过来…

用 docker-compose 启动 WebApi 和 SQL Server

本系列文章所要做出的演示架构基于 .NET Core Web Api、MSSQL、Skywalking 和 nginx &#xff0c;这些都会通过docker-compose一键创建/启动容器&#xff0c;然后用 Azure DevOps 发布上线。所以本系列文章重点并不是如何写好.NET Core&#xff0c;而是围绕着 .NET Core 的容器…

OI群论:从入门到自闭

怎么这么多阅读啊…… 这篇文章是群论&#xff08;其实只有Polya&#xff09;在信息学奥赛中计数的运用&#xff0c;并不是群论的讲义…… 群论可以说是数学中最高深的内容&#xff0c;一个oier去深入了解是不现实的。因此&#xff0c;我们只需要知道结论。 本文主要讲解群论…

2020牛客暑期多校训练营(第二场)Just Shuffle

https://ac.nowcoder.com/acm/contest/5667/J 题目大意&#xff1a;给你一个置换A&#xff0c;使得置换P^kA,让你求出置换P。 思路&#xff1a;我们根据置换A再置换z次&#xff0c;那么就等于置换p 置换z*k次,如果z*k%len0&#xff0c;那么将会回到单位序列&#xff0c;那我们…

Codeforces Round #619 (Div. 2) F. Super Jaber 多源bfs + 思维转换

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个矩阵&#xff0c;每个格子都有一个颜色kkk&#xff0c;每秒可以移动到相邻矩阵或者瞬移到同一颜色的任意矩阵。有qqq个询问&#xff0c;每次询问给出两个点&#xff0c;求从一个点到另一个点的最短时…

你可以保持沉默,但你所说的一切都将成为呈堂证供——浅谈Azure WORM保护

本文作者|Yuan Han本文来源|Reid爸的菜园子美国安然事件后&#xff0c;电子数据的合规性保存越来越受到重视&#xff1b;各国政府制定了一系列的法律&#xff0c;如美国《赛班斯法案》等&#xff0c;对于不同类型的电子数据保留期限做了严格规定&#xff1b;国内也没落后&#…

后缀自动机:从入门到放弃

写在前面 后缀自动机&#xff0c;简称SAMSAMSAM,是一种十分优秀的字符串匹(shu)配(ju)算(jie)法(gou) 字符串界的bossbossboss&#xff0c;几乎可以解决全部正常的字符串题目 至少我前前后后学了一年&#xff0c;听过444次课&#xff0c;几度怀疑自己不适合oioioi 请做好心…