洛谷P2279 [HNOI2003]消防局的设立

题目描述

2020年，人类在火星上建立了一个庞大的基地群，总共有n个基地。起初为了节约材料，人类只修建了n-1条道路来连接这些基地，并且每两个基地都能够通过道路到达，所以所有的基地形成了一个巨大的树状结构。如果基地A到基地B至少要经过d条道路的话，我们称基地A到基地B的距离为d。

由于火星上非常干燥，经常引发火灾，人类决定在火星上修建若干个消防局。消防局只能修建在基地里，每个消防局有能力扑灭与它距离不超过2的基地的火灾。

你的任务是计算至少要修建多少个消防局才能够确保火星上所有的基地在发生火灾时，消防队有能力及时扑灭火灾。

输入输出格式

输入格式：

输入文件名为input.txt。

输入文件的第一行为n （n<=1000），表示火星上基地的数目。接下来的n-1行每行有一个正整数，其中文件第i行的正整数为a[i]，表示从编号为i的基地到编号为a[i]的基地之间有一条道路，为了更加简洁的描述树状结构的基地群，有a[i]<i。

输出格式：

输出文件名为output.txt

输出文件仅有一个正整数，表示至少要设立多少个消防局才有能力及时扑灭任何基地发生的火灾。

输入输出样例

输入样例#1：

6
1
2
3
4
5

输出样例#1：

2

写树规无力，翻题解的时候发现了神奇的贪心法。

先建好树，从深度最大的地方开始贪心，能在祖先节点建消防局就在祖先节点建，否则在父节点建，否则在本身节点建。

可以用调整法证明正确性。

 /*by SilverN*/
 #include<algorithm>
 #include<iostream>
 #include<cstring>
 #include<cstdio>
 #include<cmath>
 #include<vector>
 using namespace std;
 const int mxn=;
 int read(){
     int x=,f=;char ch=getchar();
     while(ch<'' || ch>''){if(ch=='-')f=-;ch=getchar();}
     while(ch>='' && ch<=''){x=x*+ch-'';ch=getchar();}
     return x*f;
 }
 struct edge{
     int v,nxt;
 }e[mxn];
 int hd[mxn],mct=;
 void add_edge(int u,int v){
     e[++mct].v=v;e[mct].nxt=hd[u];hd[u]=mct;
     return;
 }
 int n,ans;
 int dep[mxn],fa[mxn],mxd;
 bool vis[mxn];
 vector<int>d[mxn];
 void Build(int u,int f){
     dep[u]=dep[f]+;
     mxd=max(mxd,dep[u]);
     d[dep[u]].push_back(u);
     for(int i=hd[u];i;i=e[i].nxt){
         int v=e[i].v;
         if(v==f)continue;
         fa[v]=u;
         Build(v,u);
     }
     return;
 }
 void DFS(int u,int f,int d){
     vis[u]=;
     if(d>)return;
     for(int i=hd[u];i;i=e[i].nxt){
         int v=e[i].v;
         if(v==f)continue;
         DFS(v,u,d+);
         vis[v]=;
     }
     return;
 }
 void solve(){
     int i,j;
     for(i=mxd;i;i--){
         for(j=;j<d[i].size();j++){
             int v=d[i][j];
             if(vis[v])continue;
             ++ans;
             if(!fa[v])DFS(v,,);
             else if(!fa[fa[v]])DFS(fa[v],,);
             else DFS(fa[fa[v]],,);
         }
     }
 }
 int main(){
     n=read();
     int i,j,v;
     for(i=;i<=n;i++){
         v=read();
         add_edge(i,v);
         add_edge(v,i);
     }
     Build(,);
     solve();
     printf("%d\n",ans);
     return ;
 }