【BZOJ】3835: [Poi2014]Supercomputer

题意

\(n(1 \le 1000000)\)个点的有根树,\(1\)号点为根,\(q(1 \le 1000000)\)次询问,每次给一个\(k\),每一次可以选择\(k\)个未访问的点,且父亲是访问过的,要求最少次数访问完所有的点。

分析

神题不会做。

题解

得到一个式子\(ans=max(i+ \left \lceil \frac{s[i]}{k} \right \rceil), 0 \le i \le maxh\),其中\(maxh\)是最大深度,\(s[i]\)是深度大于\(i\)的点的数量。证明如下:
定义关键层\(i\)表示拿了\(i\)次后、前\(i\)层已经拿完,以后每一次都可以拿\(k\)个(最后一次除外)。我们需要证明:1、存在关键层。2、关键层的解是最小解。
首先我们按照下面两个规则查询:

  1. 如果当前层\(i\)的结点不够\(k\),则查询完,而且如果之前有\(j < i\)层的点没查询,则查询完。
  2. 如果足够了\(k\),则在保证最终能在第\(maxh\)次遍历到第\(maxh\)层的情况下,随便选\(k\)个。

那么显然在最后一个执行1操作而且那层\(i\)不存在\(j < i\)层的点没查询的\(i\)就是关键层。由于这样的1操作至少有一个(即第1层肯定是执行的是这样的1操作),所以关键层肯定存在,而且只有一个。
至于关键层的解是否是最小解,感觉很显然,然而不会严格证明。
至于取max,不会严格证明。

最后原题可以转化为\(ans=\left \lceil max(i+\frac{s[i]}{k}) \right \rceil\),所以按照\(i+\frac{s[i]}{k}\)来斜率优化就行辣。

#include <bits/stdc++.h>
using namespace std;
inline int getint() {
	int x=0;
	char c=getchar();
	for(; c<'0'||c>'9'; c=getchar());
	for(; c>='0'&&c<='9'; x=x*10+c-'0', c=getchar());
	return x;
}
typedef long long ll;
const int N=1000005;
int a[N], d[N], q[N], c[N], ihead[N], cnt;
struct E {
	int next, to;
}e[N];
void add(int x, int y) {
	e[++cnt]=(E){ihead[x], y}; ihead[x]=cnt;
}
void dfs(int x) {
	for(int i=ihead[x]; i; i=e[i].next) {
		d[e[i].to]=d[x]+1;
		dfs(e[i].to);
	}
}
inline bool ok1(int i, int j, int k) {
	return (ll)(c[j]-c[i])*(k-i)>(ll)(c[k]-c[i])*(j-i);
}
inline bool ok2(int b, int j, int k) {
	return (ll)b*(j-k)>c[k]-c[j];
}
int main() {
	int n=getint(), Q=getint();
	for(int i=0; i<Q; ++i) {
		a[i]=getint();
	}
	d[0]=1;
	int mx=1;
	for(int i=1; i<n; ++i) {
		add(getint()-1, i);
	}
	dfs(0);
	for(int i=0; i<n; ++i) {
		++c[d[i]-1];
		mx=max(mx, d[i]);
	}
	for(int i=mx; i; --i) {
		c[i]+=c[i+1];
	}
	int *fr=q+1, *ta=q;
	for(int i=1; i<=mx; ++i) {
		for(; fr<ta && ok1(*(ta-1), i, *ta); --ta);
		*++ta=i;
	}
	for(int i=1; i<=n; ++i) {
		for(; fr<ta && ok2(i, *(fr+1), *fr); ++fr);
		d[i]=*fr+(c[*fr]+i-1)/i;
	}
	for(int i=0; i<Q; ++i) {
		printf("%d%c", a[i]>n?mx:d[a[i]], " \n"[i==Q-1]);
	}
	return 0;
}
posted @ 2015-11-22 16:07  iwtwiioi  阅读(493)  评论(0编辑  收藏  举报