洛谷 P1041 [NOIP2003 T4] 传染病控制

题目背景

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

题目描述

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

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

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

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

你的程序要针对给定的树，找出合适的切断顺序。

输入输出格式

输入格式：

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

输出格式：

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

输入输出样例

输入样例#1：

7 6
1 2
1 3
2 4
2 5
3 6
3 7

输出样例#1：

3

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

dfs~

然而我第一次写的是贪心+dfs……90分……

实际上贪心每次都删去子节点最少的位置是错的，只能过9个点……反例就是菊花树，子节点很多但是只要一步就能终止……

dfs就是暴力得不能再暴力的那种，但是我怎么也该不对…… 还是更贪心吧……

90分贪心+dfs：

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;

int n,p,x,y,fi[301],w[10001],ne[10001],cnt,fa[301],tot;
int du[301],now[301],sec[301],po[301],tot1,c[10001],mx[301];
bool b[301];

void add(int u,int v)
{
	w[++cnt]=v;ne[cnt]=fi[u];fi[u]=cnt;
}

int findd(int u)
{
	int maxx=0;
	for(int i=fi[u];i;i=ne[i])
	  if(w[i]!=fa[u])
	  {
	  	fa[w[i]]=u;du[w[i]]=du[u]+1;now[w[i]]=sec[w[i]]=++cnt;
	  	sec[u]=findd(w[i]);
	  	if(sec[w[i]]-now[w[i]]+2>maxx) po[u]=w[i],mx[u]=maxx=sec[w[i]]-now[w[i]]+2;
	  }
	return sec[u];
}

int main()
{
	scanf("%d%d",&n,&p);
	for(int i=1;i<=p;i++)
	{
		scanf("%d%d",&x,&y);
		add(x,y);add(y,x);
	}
	cnt=0;now[1]=++cnt;tot=tot1=1;findd(1);b[1]=1;
	while(tot1)
	{
		int maxx=0,nowww;c[0]=0;
		for(int i=1;i<=n;i++)
		  if(b[i] && !du[i])
		  {
		  	if(mx[i]>maxx) maxx=mx[i],nowww=i;
		  	c[++c[0]]=i;
		  }
		for(int i=1;i<=c[0];i++)
		{
			tot1--;b[c[i]]=0;
			for(int j=fi[c[i]];j;j=ne[j])
			  if(w[j]!=po[nowww] && w[j]!=fa[c[i]])
			  {
			  	b[w[j]]=1;tot++;tot1++;du[w[j]]=0;
			  }
		}
	}
	printf("%d\n",tot);
	return 0;
}