洛谷—— P1041 传染病控制

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

题目背景

近来,一种新的传染病肆虐全球。蓬莱国也发现了零星感染者,为防止该病在蓬莱国大范围流行,该国政府决定不惜一切代价控制传染病的蔓延。不幸的是,由于人们尚未完全认识这种传染病,难以准确判别病毒携带者,更没有研制出疫苗以保护易感人群。于是,蓬莱国的疾病控制中心决定采取切断传播途径的方法控制疾病传播。经过 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:



预处理子树的节点数,爆搜每棵子树,每次跟新出当前最优解
 1 #include <cstdio>
 2 #include <queue>
 3 
 4 inline void read(int &x)
 5 {
 6     x=0; register char ch=getchar();
 7     for(; ch>'9'||ch<'0'; ) ch=getchar();
 8     for(; ch>='0'&&ch<='9'; ch=getchar()) x=x*10+ch-'0';
 9 }
10 
11 const int N(326);
12 int head[N],sumedge;
13 struct Edge {
14     int v,next;
15     Edge(int v=0,int next=0):v(v),next(next){}
16 }edge[N<<1];
17 inline void ins(int u,int v)
18 {
19     edge[++sumedge]=Edge(v,head[u]);
20     head[u]=sumedge;
21     edge[++sumedge]=Edge(u,head[v]);
22     head[v]=sumedge;
23 }
24 
25 int size[N],cd[N],dad[N];
26 void DFS(int u)
27 {
28     size[u]=1; 
29     for(int v,i=head[u]; i; i=edge[i].next)
30     {
31         v=edge[i].v;
32         if(dad[u]==v) continue;
33         dad[v]=u; cd[u]++;
34         DFS(v); size[u]+=size[v];
35     }
36 }
37 bool vis[N];
38 std::queue<int>que;
39 void BFS(int s)
40 {
41     int ans=0;
42     que.push(s); vis[s]=1;
43     for(int u,v; !que.empty(); )
44     {
45         u=que.front(); que.pop();
46         ans++;
47         int t1=-1,t2=-1,t;
48         for(int i=head[u]; i; i=edge[i].next)
49         {
50             v=edge[i].v;
51             if(vis[v]) continue;
52             if(cd[v]>t2&&size[v]>=t1)
53                 t2=cd[v],t2=size[v],t=v;
54         }
55         for(int i=head[u]; i; i=edge[i].next)
56         {
57             v=edge[i].v;
58             if(vis[v]||v==t) continue;
59             que.push(v);vis[v]=1;
60         }
61     }
62     printf("%d
",ans);
63 }
64 
65 int AC()
66 {
67     int n,p; 
68     read(n),read(p);
69     for(int u,v; p--; )
70         read(u),read(v),ins(u,v);
71     DFS(1);
72     BFS(1);
73     return 0;
74 }
75 
76 int Aptal=AC();
77 int main(){;}
40分贪心,
 1 #include <cstdio>
 2 #include <vector>
 3 
 4 inline void read(int &x)
 5 {
 6     x=0; register char ch=getchar();
 7     for(; ch>'9'||ch<'0'; ) ch=getchar();
 8     for(; ch>='0'&&ch<='9'; ch=getchar()) x=x*10+ch-'0';
 9 }
10 
11 const int N(326);
12 int n,p,ans;
13 int head[N],sumedge;
14 struct Edge {
15     int v,next;
16     Edge(int v=0,int next=0):v(v),next(next){}
17 }edge[N<<1];
18 inline void ins(int u,int v)
19 {
20     edge[++sumedge]=Edge(v,head[u]);
21     head[u]=sumedge;
22     edge[++sumedge]=Edge(u,head[v]);
23     head[v]=sumedge;
24 }
25 
26 #define max(a,b) (a>b?a:b)
27 std::vector<int>vec[N];
28 int depth,size[N],dep[N],dad[N];
29 void Pre(int u)
30 {
31     size[u]=1;dep[u]=dep[dad[u]]+1;
32     depth=max(depth,dep[u]);
33     vec[dep[u]].push_back(u);
34     for(int v,i=head[u]; i; i=edge[i].next)
35     {
36         v=edge[i].v;
37         if(dad[u]==v) continue;
38         dad[v]=u; Pre(v); size[u]+=size[v];
39     }
40 }
41 bool vis[N];
42 void pushdown(int u,int x)
43 {
44     vis[u]=x;
45     for(int v,i=head[u]; i; i=edge[i].next)
46     {
47         v=edge[i].v;
48         if(dad[u]!=v) pushdown(v,x);
49     }
50 }
51 void DFS(int u,int sum)
52 {
53     ans=ans<(n-sum)?ans:(n-sum);
54     if(u>depth) return ;
55 //    if(1==size[u]) return ;
56     for(int i=0; i<vec[u].size(); ++i)
57     {
58         if(vis[vec[u][i]]) continue;
59         pushdown(vec[u][i],1);
60         DFS(u+1,sum+size[vec[u][i]]);
61         pushdown(vec[u][i],0);
62     }
63 }
64 
65 int AC()
66 {
67     read(n),read(p); ans=n;
68     for(int u,v; p--; )
69         read(u),read(v),ins(u,v);
70     Pre(1); DFS(2,0);
71     printf("%d
",ans);
72     return 0;
73 }
74 
75 int Aptal=AC();
76 int main(){;}
AC
——每当你想要放弃的时候,就想想是为了什么才一路坚持到现在。
原文地址:https://www.cnblogs.com/Shy-key/p/7503672.html