BZOJ2286:[SDOI2011]消耗战(树形DP,虚树)

BZOJ2286:[SDOI2011]消耗战(树形DP,虚树)

Description

在一场战争中,战场由n个岛屿和n-1个桥梁组成,保证每两个岛屿间有且仅有一条路径可达。现在,我军已经侦查到敌军的总部在编号为1的岛屿,而且他们已经没有足够多的能源维系战斗,我军胜利在望。已知在其他k个岛屿上有丰富能源,为了防止敌军获取能源,我军的任务是炸毁一些桥梁,使得敌军不能到达任何能源丰富的岛屿。由于不同桥梁的材质和结构不同,所以炸毁不同的桥梁有不同的代价,我军希望在满足目标的同时使得总代价最小。
侦查部门还发现,敌军有一台神秘机器。即使我军切断所有能源之后,他们也可以用那台机器。机器产生的效果不仅仅会修复所有我军炸毁的桥梁,而且会重新随机资源分布(但可以保证的是,资源不会分布到1号岛屿上)。不过侦查部门还发现了这台机器只能够使用m次,所以我们只需要把每次任务完成即可。

Input

第一行一个整数n,代表岛屿数量。

接下来n-1行,每行三个整数u,v,w,代表u号岛屿和v号岛屿由一条代价为c的桥梁直接相连,保证1<=u,v<=n且1<=c<=100000。

第n+1行,一个整数m,代表敌方机器能使用的次数。

接下来m行,每行一个整数ki,代表第i次后,有ki个岛屿资源丰富,接下来k个整数h1,h2,…hk,表示资源丰富岛屿的编号。

Output

输出有m行,分别代表每次任务的最小代价。

Sample Input

10

1 5 13

1 9 6

2 1 19

2 4 8

2 3 91

5 6 8

7 5 4

7 8 31

10 7 9

3

2 10 6

4 5 7 8 3

3 9 4 6

Sample Output

12

32

22

HINT

对于100%的数据,2<=n<=250000,m>=1,sigma(ki)<=500000,1<=ki<=n-1

Solution

虚树真是个好东西啊QAQ

推荐博客

建议对虚树的理解看第一个,构建看第二个QAQ

题解去看第二个吧我也懒得写了XD

顺带提一句因为本题特殊,所以建虚树的时候要写成57行那样,否则就把57行删掉改成58行就行了。

至于本题为什么要像57行那么写呢……

假设$x$点是$y$的祖先,如果$x$到根不连通,那么$y$到根一定不连通,所以$y$点也就没有加进去的必要了。而且加进去的话像我这样$DP$也就不对了啊啊QAQQQQ

Code

  1 #include<iostream>
  2 #include<cstring>
  3 #include<cstdio>
  4 #include<algorithm>
  5 #include<vector>
  6 #define N (250009)
  7 #define LL long long
  8 using namespace std;
  9 
 10 struct Edge{
   int to,next,len;}edge[N<<1];
 11 int n,m,k,u,v,l,dfs_num;
 12 int a[N],Depth[N],f[N][19],DFN[N];
 13 LL Min[N];
 14 int head[N],num_edge;
 15 
 16 void add(int u,int v,int l)
 17 {
 18     edge[++num_edge].to=v;
 19     edge[num_edge].next=head[u];
 20     edge[num_edge].len=l;
 21     head[u]=num_edge;
 22 }
 23 
 24 void DFS(int x,int fa)
 25 {
 26     f[x][0]=fa;
 27     for (int i=1; i<=18; ++i) f[x][i]=f[f[x][i-1]][i-1];
 28     DFN[x]=++dfs_num; Depth[x]=Depth[fa]+1;
 29     for (int i=head[x]; i; i=edge[i].next)
 30         if (edge[i].to!=fa)
 31         {
 32             Min[edge[i].to]=min(Min[x],(LL)edge[i].len);
 33             DFS(edge[i].to,x);
 34         }
 35 }
 36 
 37 int LCA(int x,int y)
 38 {
 39     if (Depth[x]<Depth[y]) swap(x,y);
 40     for (int i=18; i>=0; --i)
 41         if (Depth[f[x][i]]>=Depth[y]) x=f[x][i];
 42     if (x==y) return x;
 43     for (int i=18; i>=0; --i)
 44         if (f[x][i]!=f[y][i]) x=f[x][i], y=f[y][i];
 45     return f[x][0]; 
 46 }
 47 
 48 vector<int>E[N];
 49 void ADD(int x,int y) {E[x].push_back(y);}
 50 int stack[N],top;
 51 bool cmp(int x,int y) {
   return DFN[x]<DFN[y];}
 52 
 53 void Insert(int x)
 54 {
 55     if (top==1) {stack[++top]=x; return;}
 56     int lca=LCA(x,stack[top]);
 57     if (lca==stack[top]) return;
 58 //    if (lca==stack[top]) {stack[++top]=x; return;}
 59     while (top>1 && DFN[stack[top-1]]>=DFN[lca])
 60         ADD(stack[top-1],stack[top]), top--;
 61     if (lca!=stack[top]) ADD(lca,stack[top]), stack[top]=lca;
 62     stack[++top]=x;
 63 }
 64 
 65 void Build()
 66 {
 67     stack[top=1]=1;
 68     for (int i=1; i<=k; ++i) Insert(a[i]);
 69     while (top>=2) ADD(stack[top-1],stack[top]), top--;
 70 }
 71 
 72 LL DP(int x)
 73 {
 74     int sz=E[x].size();
 75     if (!sz) return Min[x];
 76     LL ans=0;
 77     for (int i=0; i<sz; ++i)
 78         ans+=DP(E[x][i]);
 79     E[x].clear();
 80     return min(ans,Min[x]);
 81 }
 82 
 83 int main()
 84 {
 85     Min[1]=1e18;
 86     scanf("%d",&n);
 87     for (int i=1; i<=n-1; ++i)
 88     {
 89         scanf("%d%d%d",&u,&v,&l);
 90         add(u,v,l); add(v,u,l);
 91     }
 92     DFS(1,0);
 93     scanf("%d",&m);
 94     for (int i=1; i<=m; ++i)
 95     {
 96         scanf("%d",&k);
 97         for (int j=1; j<=k; ++j) scanf("%d",&a[j]);
 98         sort(a+1,a+k+1,cmp);
 99         Build();
100         printf("%lld\n",DP(1));
101     }
102 }

转载于:https://www.cnblogs.com/refun/p/10064045.html

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/101285.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • osip状态机分析「建议收藏」

    osip状态机分析「建议收藏」OSIP的核心是系统状态机,在不同情况下,系统处于不同的状态,在某一状态下当系统发生某一个动作后(如接受或者发送了消息),状态机做相应的跳转。的状态机在不同的状态下,对某一动作的响应也是不一样的。由于SIP的复杂性,为了降低代码实现难度,也为了能够保证代码的模块化,Osip把整个SIP的实现分成了四个状态机,分别是:      (1)ICT(带invit事件的out处理)

    2022年6月15日
    23
  • excel图片链接显示为图片_怎样将图片拼接成长图

    excel图片链接显示为图片_怎样将图片拼接成长图由于公司需要统计根据门店业绩和装修档次定制扶持政策,所以要到处携带门店门头照片+业绩的报表,而数据库存储的是图片的链接。由于后台统计报表相对应的库暂时没有插入图片的功能。所以通过宏来完成将Excel中url替换为插入图片,又为了避免插入图片太多,导致Excel大小暴增,所以在选择了对应门店门头照片链接时才插入图片。具体步骤如下:1、进入对应sheet页的ViewCode一种…

    2022年9月22日
    0
  • linux golang环境安装_python环境搭建

    linux golang环境安装_python环境搭建Golang环境搭建

    2022年8月31日
    1
  • 一起学JAVA 反射学习(超详细)

    一起学JAVA 反射学习(超详细)1什么是反射?Reflection(反射)是Java程序开发语言的特征之一,它允许运行中的Java程序对自身进行检查,或者说“自审”,也有称作“自省”。反射非常强大,它甚至能直接操作程序的私有属性。我们前面学习都有一个概念,被private封装的资源只能类内部访问,外部是不行的,但这个规定被反射赤裸裸的打破了。反射就像一面镜子,它可以在运行时获取一个类的所有信息,可以获取到任何定义的信息(包括成员变量,成员方法,构造器等),并且可以操纵类的字段、方法、构造器等部分。2为什么需要反射?

    2022年6月11日
    32
  • IntelliJ IDEA常用快捷键汇总

    IntelliJ IDEA常用快捷键汇总在使用IntelliJIdea的时候,使用快捷键是必不可少的。掌握一些常用的快捷键能大大提高我们的开发效率。有些快捷键可以熟练的使用,但是还有另外一些快捷键虽然很好用,但是由于因为没有形成使用习惯或者没有理解快捷键的用法,甚至之前对一些快捷键根本没有概念,导致不会去使用。对于这些快捷键,如果能够用好,编辑代码的效率必能提高一个水平。所以在此梳理出来,加强自己的使用,形成习惯。(注:有些操作的快捷键

    2022年5月15日
    53
  • ioctl函数详解_函数concat的作用

    ioctl函数详解_函数concat的作用ioctl函数的作用特殊的read,write,当你用read,write不能完成某一功能时,就用ioctl我这里说的ioctl函数是在驱动程序里的,因为我不知道还有没有别的场合用到了ioctl,所以就规定了我们讨论的范围。为什么要写篇文章呢,是因为我前一阵子被ioctl给搞混了,这几天才弄明白它,于是在这里清理一下头脑。什么是ioctl。ioctl是设备驱动程序中对设备…

    2022年10月18日
    0

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号