洛谷1041 传染病控制

NOIP2003 

原题地址

https://www.luogu.org/problem/show?pid=1041

搜索

题目描述 Description

【问题背景】

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

【问题描述】

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

输入描述 Input Description

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

输出描述 Output Description

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

样例输入 Sample Input

7 6

1 2

1 3

2 4

2 5

3 6

3 7

样例输出 Sample Output

3

数据范围及提示 Data Size& Hint

n<=300

 

解题思路

乍一看好像是树P加贪心，然而事实上树P还差不多，贪心就扯远了。仔细想一想这题压根不符合贪心的性质，有人举反例说菊花形的就显然不符合。所以，排除了贪心，再看看这良心的数据范围，跑一跑DFS就可以了。

1. 首先建树（双向边），然后跑出每个点子树的size（包含本身），同时确定父子关系，开一个vector存每一层（同一深度）的点的编号。

2. 以深度deep和当前可感染人数cnt为关键词进行DFS：枚举深度为deep的所有点，如果满足条件（没有被阻断），截断父亲节点连向该点的边——将该点极其所有儿子节点打上标记（为啥是儿子节点而不是子树上的所有节点？我们知道dfs是一层一层来的，判断是否能进行dfs的条件是由上一层的选择情况决定的，所以只与上一层的操作有关，所以只给儿子节点打标记就好了），同时可感染人数cnt减去该点子树规模size[i]，进行下一层的搜索。注意这是一个回溯的过程，dfs之后要清空所打的标记。每一层DFS时要判断当前cnt是否能更新答案。

3. 输出答案，return 0。

参考代码

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<vector>
using namespace std;
int num=0,ans,n,m,v[350],size[350],fa[350];
bool vis[350];
vector<int>d[300];
struct mc
{
    int x,y,ne;
}e[10086];
 
void put(int x,int y)
{
    num++;
    e[num].x=x;
    e[num].y=y;
    e[num].ne=v[x];
    v[x]=num;
}
 
void make(int x,int deep)
{
    size[x]=1;
    if (v[x]==1&&x!=1)return;
    for (int i=v[x];i;i=e[i].ne)
    {
       int y=e[i].y;
       if (!vis[y])
       {
           fa[y]=x;
           d[deep+1].push_back(y);
           vis[y]=1;
           make(y,deep+1);
           size[x]+=size[y];
       }
    }
}
 
void clean(int x,int t)
{
    vis[x]=t;
    for (int i=v[x];i;i=e[i].ne)
    {
       int y=e[i].y;
       if (y!=fa[x])
       {
           vis[y]=t;
           clean(y,t);
       }
    }
}
 
void dfs(int deep,int cnt)
{
    ans=min(ans,cnt);
    for (inti=0;i<d[deep].size();i++)
    {
       int y=d[deep][i];
       if (!vis[y])
       {
           clean(y,1);
           dfs(deep+1,cnt-size[y]);
           clean(y,0);
       }
    }
}
 
int main()
{
    int x,y;
    cin>>n>>m;
    for (int i=1;i<=m;i++)
    {
       cin>>x>>y;
       put(x,y);
       put(y,x);
    }
    memset(vis,false,sizeof(vis));
    d[1].push_back(1);
    vis[1]=1;
    make(1,1);
    ans=n;
    memset(vis,false,sizeof(vis));
    vis[1]=1;
    dfs(2,n);
    cout<<ans;
    return 0;
}