疫情控制 blockade
疫情控制 blockad
题目描述
H 国有 n 个城市,这 n 个城市用 n-1 条双向道路相互连通构成一棵树, 1 号城市是首都, 也是树中的根节点。
H 国的首都爆发了一种危害性极高的传染病。当局为了控制疫情,不让疫情扩散到边境 城市(叶子节点所表示的城市),决定动用军队在一些城市建立检查点,使得从首都到边境 城市的每一条路径上都至少有一个检查点,边境城市也可以建立检查点。但特别要注意的是, 首都是不能建立检查点的。 现在,在 H 国的一些城市中已经驻扎有军队,且一个城市可以驻扎多个军队。一支军队可以在有道路连接的城市间移动,并在除首都以外的任意一个城市建立检查点,且只能在 一个城市建立检查点。一支军队经过一条道路从一个城市移动到另一个城市所需要的时间等 于道路的长度(单位:小时)。
请问最少需要多少个小时才能控制疫情。注意:不同的军队可以同时移动。
输入
第一行一个整数 n,表示城市个数。
接下来的 n-1 行,每行 3 个整数,u、v、w,每两个整数之间用一个空格隔开,表示从 城市 u 到城市 v 有一条长为 w 的道路。数据保证输入的是一棵树,且根节点编号为 1。
接下来一行一个整数 m,表示军队个数。
接下来一行 m 个整数,每两个整数之间用一个空格隔开,分别表示这 m 个军队所驻扎 的城市的编号。
输出
共一行,包含一个整数,表示控制疫情所需要的最少时间。如果无法控制疫情则输出-1。
样例输入
4
1 2 1
1 3 2
3 4 3
2
2 2
样例输出
3
提示
【样例解释】
第一支军队在 2 号点设立检查点,第二支军队从 2 号点移动到 3 号点设立检查点,所需 时间为 3 个小时
【数据范围】
保证军队不会驻扎在首都。
对于 20%的数据,2≤ n≤ 10;
对于 40%的数据,2 ≤n≤50,0<w <105;
对于 60%的数据,2 ≤ n≤1000,0<w <106;
对于 80%的数据,2 ≤ n≤10,000;
对于 100%的数据,2≤m≤n≤50,000,0<w <109。
来源
solution
容易发现答案具有单调性。
二分答案,然后尝试构造最优方案。
我们应该让所有点尽可能往上跳,如果步数充裕还可以去到别的子树去驻扎。
可以用倍增的方法算出每个点会跳到哪里,如果能跳到根就记下还能走几步。
那么跳不到根的,他的贡献是确定的,所以可以先确定根的那些儿子是已被处理好了。
然后我们把剩余的点按剩余步数排序。记剩余步数为res
如果一个点跳不回自己来的儿子,并且这个儿子还没处理好,那就让他跳回去。
因为让另一个res>它的点去肯定不优。
接下来就对每一个还未被占领的儿子按占领代价(也就是根到它的长度)排序。
贪心匹配即可。
总结:先找性质,二分完构造最优解。细节需想清楚。
1 #include<cstdio> 2 #include<iostream> 3 #include<cstdlib> 4 #include<cstring> 5 #include<algorithm> 6 #include<cmath> 7 #define maxn 50004 8 #define ll long long 9 using namespace std; 10 int n,m,head[maxn],tot,f[maxn][22],be[maxn]; 11 int match[maxn],top; 12 ll l[maxn][22],len[maxn]; 13 bool flag[maxn]; 14 struct node{ 15 int v,nex,w; 16 }e[maxn*2]; 17 struct no{ 18 int st,ed; 19 ll res; 20 }s[maxn]; 21 void lj(int t1,int t2,int t3){ 22 e[++tot].v=t2;e[tot].w=t3;e[tot].nex=head[t1];head[t1]=tot; 23 } 24 void dfs(int k,int fa){ 25 f[k][0]=fa; 26 if(fa==1)be[k]=k; 27 else be[k]=be[fa]; 28 for(int i=head[k];i;i=e[i].nex){ 29 if(e[i].v!=fa){ 30 l[e[i].v][0]=e[i].w; 31 dfs(e[i].v,k); 32 } 33 } 34 } 35 void work(int k){ 36 bool F=0,fl=0; 37 for(int i=head[k];i;i=e[i].nex){ 38 if(e[i].v==f[k][0])continue; 39 work(e[i].v); 40 if(!fl)fl=1,F=1; 41 F&=flag[e[i].v]; 42 } 43 flag[k]|=F; 44 } 45 bool Mr(no a,no b){ 46 return a.res<b.res; 47 } 48 bool cmp(int a,int b){ 49 return len[a]<len[b]; 50 } 51 bool pd(ll mid){ 52 53 for(int i=1;i<=n;i++)flag[i]=0; 54 for(int i=1;i<=m;i++){ 55 int k=s[i].st;ll la=mid; 56 for(int x=20;x>=0;x--){ 57 if(f[k][x]&&l[k][x]<=la){ 58 la-=l[k][x],k=f[k][x]; 59 } 60 } 61 s[i].ed=k;s[i].res=la; 62 if(k>1)flag[k]=1; 63 } 64 work(1); 65 sort(s+1,s+m+1,Mr); 66 for(int i=1;i<=m;i++){ 67 if(s[i].ed!=1||flag[be[s[i].st]])continue; 68 int b=be[s[i].st]; 69 if(s[i].res<len[b])flag[b]=1,s[i].ed=b; 70 } 71 top=0; 72 for(int i=head[1];i;i=e[i].nex){ 73 if(flag[e[i].v])continue; 74 match[++top]=e[i].v; 75 } 76 sort(match+1,match+top+1,cmp); 77 int j=1; 78 for(int i=1;i<=top;i++){ 79 80 int L=len[match[i]]; 81 while((s[j].ed!=1||s[j].res<L)&&j<=m)j++; 82 j++; 83 if(j>m+1)return 0; 84 } 85 return 1; 86 } 87 int main(){ 88 cin>>n; 89 for(int i=1,t1,t2,t3;i<n;i++){ 90 scanf("%d%d%d",&t1,&t2,&t3); 91 lj(t1,t2,t3);lj(t2,t1,t3); 92 } 93 dfs(1,0); 94 for(int i=head[1];i;i=e[i].nex)len[e[i].v]=e[i].w; 95 for(int j=1;j<=20;j++) 96 for(int i=1;i<=n;i++){ 97 f[i][j]=f[f[i][j-1]][j-1]; 98 l[i][j]=l[i][j-1]+l[f[i][j-1]][j-1]; 99 } 100 cin>>m; 101 for(int i=1;i<=m;i++)scanf("%d",&s[i].st); 102 ll l=0,r=1e16; 103 while(l<r){ 104 ll mid=l+r>>1; 105 if(pd(mid))r=mid; 106 else l=mid+1; 107 } 108 if(l==1e16)puts("-1"); 109 else cout<<l<<endl; 110 return 0; 111 } 112