[HNOI2003]消防局的设立 树形dp // 贪心
https://www.luogu.org/problemnew/show/P2279
一开始就想到了贪心的方法,不过一直觉得不能证明。
贪心的考虑是在深度从深到浅遍历每个结点的过程中,对于每个没有覆盖的结点选择覆盖他的祖父结点。
仔细想想觉得这是正确的。
在实现的过程中有一个小技巧是o[i]记录i结点距离消防局最近的距离,如果o[i] > 2则需要在他的祖父结点建立一个消防站。用这种方法可以很方便的判断兄弟节点是否被覆盖。
一个细节是要给根节点1建立两个虚结点N + 1和N + 2作为他的父亲结点和祖父结点,不然在1结点寻找祖父节点的时候容易出现问题。
用这样的方法也可以在常数很小的时间内实现树上半径k的覆盖问题
#include <map>
#include <set>
#include <ctime>
#include <cmath>
#include <queue>
#include <stack>
#include <vector>
#include <string>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <sstream>
#include <iostream>
#include <algorithm>
#include <functional>
using namespace std;
inline int read(){int now=;register char c=getchar();for(;!isdigit(c);c=getchar());
for(;isdigit(c);now=now*+c-'',c=getchar());return now;}
#define For(i, x, y) for(int i=x;i<=y;i++)
#define _For(i, x, y) for(int i=x;i>=y;i--)
#define Mem(f, x) memset(f,x,sizeof(f))
#define Sca(x) scanf("%d", &x)
#define Sca2(x,y) scanf("%d%d",&x,&y)
#define Sca3(x,y,z) scanf("%d%d%d",&x,&y,&z)
#define Scl(x) scanf("%lld",&x);
#define Pri(x) printf("%d\n", x)
#define Prl(x) printf("%lld\n",x);
#define CLR(u) for(int i=0;i<=N;i++)u[i].clear();
#define LL long long
#define ULL unsigned long long
#define mp make_pair
#define PII pair<int,int>
#define PIL pair<int,long long>
#define PLL pair<long long,long long>
#define pb push_back
#define fi first
#define se second
typedef vector<int> VI;
const double eps = 1e-;
const int maxn =;
const int INF = 0x3f3f3f3f;
const int mod = 1e9 + ;
int N,M,K;
int fa[maxn];
int deep[maxn];
int a[maxn];
int o[maxn];
bool cmp(int a,int b){
return deep[a] > deep[b];
}
int main(){
Sca(N); o[] = INF;
fa[] = N + ;
fa[N + ] = N + ;
o[N + ] = o[N + ] = INF;
for(int i = ; i <= N; i ++) a[i] = i;
for(int i = ; i <= N ; i ++){
Sca(fa[i]);
deep[i] = deep[fa[i]] + ;
o[i] = INF;
}
sort(a + ,a + + N,cmp);
int ans = ;
for(int i = ; i <= N ; i ++){
int u = a[i];
int v = fa[u],w = fa[fa[u]];
o[u] = min(o[u],min(o[v] + ,o[w] + ));
if(o[u] > ){
o[w] = ;
o[fa[w]] = min(o[fa[u]],);
o[fa[fa[w]]] = min(o[fa[fa[w]]],);
ans++;
}
}
Pri(ans);
return ;
}
贪心
当然这道题出现在dp专题里,事实上也是一个很硬核的树dp,在没有想到或者证明贪心的时候,这样明显的树dp也不失为一种做法。
思路借鉴了luogu的题解,开始对于这样的树dp一直想着两次树dp,第一次记录根节点的answer,第二次利用根节点拓展到整棵树的answer,导致开始的思路仅限于dp[maxn][3]来记录这个点,这个点的儿子,这个点的孙子的消防站情况,使得第一次dfs就遇到了问题,可能可以做,但有些麻烦。
看了题解之后茅塞顿开,事实上树dp并不需要套路的总是两边dp,换一种思路同时维护两端的情况,就是除了原本的3个状态以外,重新记录dp[i][3]表示所有儿子都被覆盖的情况以及dp[i][4]所有孙子都被覆盖的情况。
仅用一遍dfs就可以解决这个问题。
#include <map>
#include <set>
#include <ctime>
#include <cmath>
#include <queue>
#include <stack>
#include <vector>
#include <string>
#include <cstdio>
#include<algorithm>
#include <cstdlib>
#include <cstring>
#include <sstream>
#include <iostream>
#include <functional>
using namespace std;
inline int read(){int now=;register char c=getchar();for(;!isdigit(c);c=getchar());
for(;isdigit(c);now=now*+c-'',c=getchar());return now;}
#define For(i, x, y) for(int i=x;i<=y;i++)
#define _For(i, x, y) for(int i=x;i>=y;i--)
#define Mem(f, x) memset(f,x,sizeof(f))
#define Sca(x) scanf("%d", &x)
#define Sca2(x,y) scanf("%d%d",&x,&y)
#define Sca3(x,y,z) scanf("%d%d%d",&x,&y,&z)
#define Scl(x) scanf("%lld",&x);
#define Pri(x) printf("%d\n", x)
#define Prl(x) printf("%lld\n",x);
#define CLR(u) for(int i=0;i<=N;i++)u[i].clear();
#define LL long long
#define ULL unsigned long long
#define mp make_pair
#define PII pair<int,int>
#define PIL pair<int,long long>
#define PLL pair<long long,long long>
#define pb push_back
#define fi first
#define se second
typedef vector<int> VI;
const double eps = 1e-;
const int maxn = ;
const int INF = 0x3f3f3f3f;
const int mod = 1e9 + ;
int N,M,K;
struct Edge{
int to,next;
}edge[maxn * ];
int head[maxn],tot;
void init(){
for(int i = ; i <= N ; i ++) head[i] = -;
tot = ;
}
void add(int u,int v){
edge[tot].to = v;
edge[tot].next = head[u];
head[u] = tot++;
}
int dp[maxn][]; // 0这是,1儿子是,2孙子是,3儿子全被覆盖,4孙子全被覆盖
int Min(int a,int b){
return a > b?b:a;
}
int Min(int a,int b,int c){
return Min(Min(a,b),c);
}
int Min(int a,int b,int c,int d){
return Min(Min(a,b,c),d);
}
int Min(int a,int b,int c,int d,int e){
return Min(Min(a,b,c),Min(d,e));
}
void dfs(int t){
if(head[t] == -){
dp[t][] = ;
dp[t][] = dp[t][] = INF;
dp[t][] = dp[t][] = ;
return;
}
dp[t][] = ;
int MAXSON = INF;
int MAXGS = INF;
for(int i = head[t]; ~i; i = edge[i].next){
int v = edge[i].to;
dfs(v);
dp[t][] += Min(dp[v][],dp[v][],dp[v][],dp[v][],dp[v][]);
dp[t][] += Min(dp[v][],dp[v][],dp[v][],dp[v][]);
MAXSON = Min(dp[v][] - Min(dp[v][],dp[v][],dp[v][],dp[v][]),MAXSON);
MAXGS = Min(dp[v][] - Min(dp[v][],dp[v][],dp[v][]),MAXGS);
dp[t][] += Min(dp[v][],dp[v][],dp[v][]);
dp[t][] += Min(dp[v][],dp[v][],dp[v][]);
dp[t][] += Min(dp[v][],dp[v][],dp[v][],dp[v][]);
}
dp[t][] += MAXSON;
dp[t][] += MAXGS;
}
int main(){
Sca(N); init();
for(int i = ; i <= N; i ++){
int x; Sca(x);
add(x,i);
}
dfs();
cout << Min(dp[][],dp[][],dp[][]);
return ;
}
[HNOI2003]消防局的设立 树形dp // 贪心的更多相关文章
- 【BZOJ1217】[HNOI2003]消防局的设立 树形DP
[BZOJ1217][HNOI2003]消防局的设立 Description 2020年,人类在火星上建立了一个庞大的基地群,总共有n个基地.起初为了节约材料,人类只修建了n-1条道路来连接这些基地, ...
- P2279 [HNOI2003]消防局的设立[树形dp]
题目描述 2020年,人类在火星上建立了一个庞大的基地群,总共有n个基地.起初为了节约材料,人类只修建了n-1条道路来连接这些基地,并且每两个基地都能够通过道路到达,所以所有的基地形成了一个巨大的树状 ...
- bzoj1217: [HNOI2003]消防局的设立 [树形dp]
Description 2020年,人类在火星上建立了一个庞大的基地群,总共有n个基地.起初为了节约材料,人类只修建了n-1条道路来连接这些基地,并且每两个基地都能够通过道路到达,所以所有的基地形成了 ...
- luogu 2279 [HNOI2003]消防局的设立 树形dp
就是细节多一些,思路都非常常规. Code: #include <bits/stdc++.h> #define N 1005 #define inf 1061109567 #define ...
- [HNOI2003] 消防局的设立 - 树形dp
仍然是点覆盖集问题,但覆盖半径变成了\(2\) 延续上一题的思路,只是式子更加复杂了 想体验一下min_element大法于是不想优化了 #include <bits/stdc++.h> ...
- [HNOI2003]消防局的设立 (贪心)
[HNOI2003]消防局的设立 题目描述 2020年,人类在火星上建立了一个庞大的基地群,总共有n个基地.起初为了节约材料,人类只修建了n-1条道路来连接这些基地,并且每两个基地都能够通过道路到达, ...
- BZOJ 1217: [HNOI2003]消防局的设立( 贪心 )
一个简单的贪心, 我们只要考虑2个消防局设立的距离为5时是最好的, 因为利用最充分. 就dfs一遍, 再对根处理一下就可以了. 这道题应该是SGU某道题的简化版...这道题距离只有2, 树型dp应该也 ...
- [luogu]P2279 [HNOI2003]消防局的设立[贪心]
[luogu]P2279 [HNOI2003]消防局的设立 题目描述 2020年,人类在火星上建立了一个庞大的基地群,总共有n个基地.起初为了节约材料,人类只修建了n-1条道路来连接这些基地,并且每两 ...
- BZOJ1217: [HNOI2003]消防局的设立
BZOJ1217: [HNOI2003]消防局的设立 Description 2020年,人类在火星上建立了一个庞大的基地群,总共有n个基地. 起初为了节约材料,人类只修建了n-1条道路来连接这些基地 ...
随机推荐
- k8s(一) kubeadm简单集群初始化
写给想入门kubernetes的同学们 # 系统版本 [root@master ~]# cat /etc/os-release NAME="CentOS Linux" VERSIO ...
- P1091 合唱队形 最长上升子序列
思路:最长上升子序列 正着做一遍 倒着做一遍 然后 取最大值 #include<bits/stdc++.h> using namespace std; const int maxn=105 ...
- Django+Xadmin打造在线教育系统(二)
基于xadmin的后台管理 先使用pip进行安装xadmin及其依赖包 pip install django-adminx 安装完成后卸载xadmin,保留依赖包即可 pip uninstall dj ...
- Codeforces543 B. Destroying Roads
传送门:>Here< 题意:给出一张无向图(边权为1),并给出两对起点和终点以及距离:s1,t1,l1; s2,t2,l2; 要求删除尽量多的边,使得dis(s1,t1)<=l1, ...
- 怎么用Verilog描述双向IO口
在很多情况下,需要使用双向IO.不过最好谨慎使用,在top层使用.网上很多描述的代码甚至是不可以综合并且有语法错误的,还是老实自己写个模块吧. 如果你需要一个口既做输入端口也做输出端口,那么你就需要去 ...
- Android assets的一个bug
摘要 Android assets目录下资源文件超过1MB的问题. 由于要显示一些奇奇怪怪的日文字符,我们在应用里放了一个字库文件,譬如叫做jp.ttf,放在assets目录下打包. 开发.调试一切正 ...
- [hgoi#2019/3/21]NOIP&NOI赛后总结
前言 今天做的是是2010年提高组和NOI的题目,做过几道原题,但是还是爆炸了,我真的太弱了. t1-乌龟棋 https://www.luogu.org/problemnew/show/P1541 这 ...
- PHP安装-phpMyAdmin+Discuz
PHP安装-phpMyAdmin+Discuz基于Apache和MySQL安装完成之后继续安装PHP.以构建LAMP动态网站平台.http:./configure --prefix=/usr/loca ...
- 买卖股票的最佳时机II
题目描述 给定一个数组,它的第 i 个元素是一支给定股票第 i 天的价格. 如果你最多只允许完成一笔交易(即买入和卖出一支股票),设计一个算法来计算你所能获取的最大利润. 注意你不能在买入股票前卖出股 ...
- 让Mac 可以使用mysql -u用户直接连接数据库
在执行完安装版本的mysql数据库后,会发现执行mysql还是会出现 command not found的错误:解决方案 方案1.设置软连接到/usr/local/bin下在命令行下输入如下 ln - ...