NOIP 2003 洛谷P1041 传染病控制

题目背景

    近来，一种新的传染病肆虐全球。蓬莱国也发现了零星感染者，为防止该病在蓬莱国大范围流行，该国政府决定不惜一切代价控制传染病的蔓延。不幸的是，由于人们尚未完全认识这种传染病，难以准确判别病毒携带者，更没有研制出疫苗以保护易感人群。于是，蓬莱国的疾病控制中心决定采取切断传播途径的方法控制疾病传播。经过 WHO（世界卫生组织）以及全球各国科研部门的努力，这种新兴传染病的传播途径和控制方法已经研究消楚，剩下的任务就是由你协助蓬莱国疾控中心制定一个有效的控制办法。

题目描述

研究表明，这种传染病的传播具有两种很特殊的性质；

第一是它的传播途径是树型的，一个人X只可能被某个特定的人Y感染，只要Y不得病，或者是XY之间的传播途径被切断，则X就不会得病。

第二是，这种疾病的传播有周期性，在一个疾病传播周期之内，传染病将只会感染一代患者，而不会再传播给下一代。

这些性质大大减轻了蓬莱国疾病防控的压力，并且他们已经得到了国内部分易感人群的潜在传播途径图（一棵树）。但是，麻烦还没有结束。由于蓬莱国疾控中心人手不够，同时也缺乏强大的技术，以致他们在一个疾病传播周期内，只能设法切断一条传播途径，而没有被控制的传播途径就会引起更多的易感人群被感染（也就是与当前已经被感染的人有传播途径相连，且连接途径没有被切断的人群）。当不可能有健康人被感染时，疾病就中止传播。所以，蓬莱国疾控中心要制定出一个切断传播途径的顺序，以使尽量少的人被感染。你的程序要针对给定的树，找出合适的切断顺序。

 

输入格式

    输入格式的第一行是两个整数n（1≤n≤300）和p。接下来p行，每一行有两个整数i和j，表示节点i和j间有边相连（意即，第i人和第j人之间有传播途径相连）。其中节点1是已经被感染的患者。

 

输出格式

    只有一行，输出总共被感染的人数。

样例输入

 7 6

 1 2

 1 3

 2 4

 2 5

 3 6

 3 7

 

样例输出

 3

 

题目链接

 

看数据范围可知，这道题搜索就能过去。但是显然我们不能直接爆搜，比如枚举哪些点删哪些边，那样肯定不行。

换一种思考方式，因为在一个感染周期内(其实就是一个时间单位了)，病毒只会传染一次。也就是说，病毒每次都是向下传染一层。而我们每次切边为了满足最优性，肯定要切掉更多的点，所以最优显然是在当前深度删边。

所以我们可以枚举每一层删除哪条边，记录下删除的点数然后更新答案。

我们可以先同时预处理深度和每个节点的子节点数，再用邻接表来存每一层的点。之后就逐层搜索，枚举每层删除的边，每次删边和还原边都要遍历一遍这条边下面的子树。由于我们已经处理好了每个节点的子节点数，所以可以O(1)算出当前情况已经删了多少点，然后每一层枚举结束后更行答案就好了。

 

#include<cstdio>
#include<iostream>
#include<queue>
#include<algorithm>
#include<cmath>
using namespace std;
inline int read()
{
    int sum=0;
    char ch=getchar();
    while(ch<'0'||ch>'9')
        ch=getchar();
    while(ch>='0'&&ch<='9')
    {
        sum=sum*10+ch-'0';
        ch=getchar();
    }
    return sum;
}
void write(int x)
{
    if(x>9) write(x/10);
    putchar(x%10+'0');
}
struct tree
{
    int next,node;
}h[605];
struct edge
{
    int next,node;
}Dpc[605];
int Head[305],son[305],Dhead[305],Depth[305];//Head&&Dhead用于邻接表
int n,m,top=0,tot=0,ans=192608170;//son记录子节点数，Depth记录深度
bool visit[305];//判断是否切掉。
inline void add(int u,int v)//邻接表存图
{
    h[++tot].next=Head[u];
    h[tot].node=v;
    Head[u]=tot;
}
inline void addd(int u,int v)//邻接表存每层的点数
{
    Dpc[++top].next=Dhead[u];
    Dpc[top].node=v;
    Dhead[u]=top;
}
int cts(int pos)//预处理子节点个数以及深度
{
    for(register int i=Head[pos];i;i=h[i].next)
    {
        int v=h[i].node;
        if(!Depth[v])
        {
            Depth[v]=Depth[pos]+1;
            addd(Depth[v],v);//存点
            son[pos]+=cts(v);
        }
    }
    return son[pos]+1;
}
void cut(int x,int flag)
{
    visit[x]=flag;//删边&还原边
    for(register int i=Head[x];i;i=h[i].next)
    {
        int v=h[i].node;
        if(Depth[v]==Depth[x]+1)
            cut(v,flag);
    }
}
void dfs(int deep,int cnt)//搜索每一层
{
    for(register int i=Dhead[deep];i;i=Dpc[i].next)
    {
        int v=Dpc[i].node;
        if(!visit[v])
        {
            cut(v,1);
            dfs(deep+1,cnt-son[v]-1);
            cut(v,0);
        }
    }
    ans=min(ans,cnt);//更新答案
}
int main()
{
    int x,y;
    n=read();
    m=read();
    for(register int i=1;i<=m;++i)
    {
        x=read();
        y=read();
        add(x,y);
        add(y,x);
    }
    visit[1]=1;
    Depth[1]=1;
    cts(1);
    dfs(2,n);//由于存的深度的关系(点1深度为1)所以我从第二层开始枚举
    cout<<ans<<endl;
    return 0;
}