jzoj 3104. 【NOIP2012提高组】疫情控制

这个博客讲得好（反正我是看懂了）
这题题目全是图片，就不贴了。。。
这题的思想很明显：

就是二分+贪心判断

但是细节很多，很坑的地方。
我们对于一个二分出来的答案，算出以下几个东东：
1.如果是能到达根节点的，那么求出它是从根节点的哪个儿子跳来的以及它还剩余的步数，存进can数组里。
2.如果不能，就将它能跳到的最高的那个点打一个标记（贪心思想）
而后，我们找出哪些根节点的儿子还没有被控制住，并存入na数组里。
然后对于这两个数组贪心乱搞即可。
我们当然知道，剩余步数小的控制剩余步数较小的是最优的，所以我们要排序。
然后对于当先点can[i]，如果它从x跳到根节点，但是x还没有被标记的话，我们就明白了：
x到根节点的距离一定比can[i]所剩余的步数要大，以小换大是个不错的选择，所以一定要判！！！

	sort(na+1,na+nat+1,cmp);
	sort(can+1,can+tot+1,cmp);
	int kk=1;
	for (int i=1;i<=tot;i++)
	{
		while (!bz[can[i].x] && i<=tot) bz[can[i].x]=1,i++;//这东东当时没判结果WA了。。。
		if (i>tot) break;//这是附带上面的赠送品。。。
		while (bz[na[kk].x]) kk++;
		if (can[i].len<na[kk].len) continue;
		bz[na[kk].x]=1;kk++;//控制了的节点一定要标记一下，没标记就WA了。。。
		if (kk>nat) break;
	}

就是上面的几个错误，搞低了我的正确率！！！

code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 100010
#define ll long long
using namespace std;
struct node{int v,fr,w;}e[N<<1];
struct ha{int x,len;}can[N<<1],na[N<<1];
int n,m,a[N],tail[N],fa[N][16],dep[N],cnt=0,tot,nat;
ll dis[N][16];
bool bz[N];

inline int read()
{
	int x=0; char c=getchar();
	while (c<'0' || c>'9') c=getchar();
	while (c>='0' && c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
	return x;
}

void add(int u,int v,int w) {e[++cnt]=(node){v,tail[u],w}; tail[u]=cnt;}

void dfs(int x)
{
	dep[x]=dep[fa[x][0]]+1;
	for (int i=0;fa[fa[x][i]][i];i++)
	{
		fa[x][i+1]=fa[fa[x][i]][i];
		dis[x][i+1]=dis[x][i]+dis[fa[x][i]][i];
	}
	for (int p=tail[x],v;p;p=e[p].fr)
	{
		v=e[p].v;
		if (v==fa[x][0]) continue;
		fa[v][0]=x,dis[v][0]=e[p].w;
		dfs(v);
	}
}

void ad(int x,int len) {can[++tot]=(ha){x,len};}

bool dg(int x)
{
	if (bz[x]) return 1;
	if (!e[tail[x]].fr) return 0;
	for (int p=tail[x],v;p;p=e[p].fr)
	{
		v=e[p].v;
		if (v==fa[x][0]) continue;
		if (!dg(v)) return 0;
	}
	bz[x]=1;
}

int cmp(ha x,ha y) {return x.len<y.len;}

int cmp1(ha x,ha y) {return x.x<y.x;}

bool check(int mid)
{
	tot=nat=0;
	memset(bz,0,sizeof(bz));
	for (int i=1,x,hav;i<=m;i++)
	{
		x=a[i],hav=0;
		for (int j=15;j>=0;j--)
			if (hav+dis[x][j]<=mid && dep[fa[x][j]]>1)
				hav+=dis[x][j],x=fa[x][j];
		if (fa[x][0]==1 && hav+dis[x][0]<=mid) can[++tot]=(ha){x,mid-hav-dis[x][0]};
		else bz[x]=1;
	}
	for (int p=tail[1],v;p;p=e[p].fr)
		if (!dg(v=e[p].v)) na[++nat]=(ha){v,dis[v][0]};
	sort(na+1,na+nat+1,cmp);
	sort(can+1,can+tot+1,cmp);
	int kk=1;
	for (int i=1;i<=tot;i++)
	{
		while (!bz[can[i].x] && i<=tot) bz[can[i].x]=1,i++;
		if (i>tot) break;
		while (bz[na[kk].x]) kk++;
		if (can[i].len<na[kk].len) continue;
		bz[na[kk].x]=1;kk++;
		if (kk>nat) break;
	}
	while (bz[na[kk].x]) kk++;
	if (kk>nat) return 1;
	else return 0;
}

int main()
{
	freopen("blockade.in","r",stdin);
	freopen("blockade.out","w",stdout);
	n=read();
	for (int i=1,u,v,w;i<n;i++)
		u=read(),v=read(),w=read(),add(u,v,w),add(v,u,w);
	dfs(1);
	m=read();
	for (int i=1;i<=m;i++) a[i]=read();
	int l=0,r=200000,mid;
	while (l<=r)
	{
		mid=l+r>>1;
		if (check(mid)) r=mid-1;
		else l=mid+1;
	}
	printf("%d
",l);
	return 0;
}