更好的阅读体验：https://www.cnblogs.com/To-Carpe-Diem/p/19635050

大意

给定一个序列，每个数有两种颜色，要找到一个位置，要找到一个位置使得 $\min \{ f(p), g(p) \}$ 尽可能大.

思路

我们在温度 $T$ 下，能参加的冰人必须是 $y \le T$ 的才能参赛，火人必须是 $y \ge T$ 才能参赛，我们定义能参赛的冰人的能量和为 $f(T)$，同理，火人为 $g(T)$。

那么我们的答案一定是 $2 \times \min \{ f(T), g(T)\}$，对于这个图，画出来是下面这样的：

那么，取 $\min$ 之后的图像是这样的：

就是绿色的部分，那么我们这样的话，答案一定在交点附近，那么我们在求这个交点的过程，可以尝试二分 + 树状数组，由于一个是递增，一个是递减，分别维护前缀和和后缀和即可，那么这样二分到交点 $p$ 的话，只需要看 $p$ 和 $p + 1$ 位置，相当于是冰人和火人的位置，火人的位置比较特殊，因为是可能后面还有段连续的区间，需要继续从 $p + 1$ 向右跳到最后一个合法的区间。

而上述的做法，不难发现时间复杂度是 $Q \log^2 Q$ 的，无法通过此题，只能有 $60pts$，于是我们考虑优化。

有个比较经典的思路（不过我也是刚学会的），在树状数组上倍增，这样的话，我们只需要单 $\log$ 的复杂度去跳交点，以及查找 $p + 1$ 的后面最后一个满足情况的点即可。

代码

#include<bits/stdc++.h>
using namespace std;

#define ll long long
const int MAXN = 2 * 1e6 + 5;
ll t1[MAXN], t2[MAXN];
int lsh[MAXN], tot;
ll sum = 0;

struct node{
	int op, t, x, y;
}q[MAXN];

int lowbit(int x){
	return x & -x;
}

void add(ll *bit, int x, ll k){
//	if(x <= 0) return;
	while(x <= tot + 1){
		bit[x] += k;
		x += lowbit(x);
	}
}

int Q;

int main(){
	cin.tie(0) -> ios::sync_with_stdio(0);
	
	cin >> Q;

	for(int i = 1;i <= Q;i ++){
		cin >> q[i].op;
		if(q[i].op == 1){
			cin >> q[i].t >> q[i].x >> q[i].y;
			lsh[++ tot] = q[i].x;
		}
		else{
			cin >> q[i].t;
		}
	}

	sort(lsh + 1, lsh + tot + 1);

	tot = unique(lsh + 1, lsh + tot + 1) - (lsh + 1);

	for(int i = 1;i <= Q;i ++){
		if(q[i].op == 1){
			q[i].x = lower_bound(lsh + 1, lsh + tot + 1, q[i].x) - lsh;
		}
	}

	for(int i = 1;i <= Q;i ++){
		if(q[i].op == 1){
			if(q[i].t == 0) add(t1, q[i].x, q[i].y);
			else add(t2, q[i].x + 1, q[i].y), sum += q[i].y;
		}
		else{
			if(q[q[i].t].t == 0) add(t1, q[q[i].t].x, -q[q[i].t].y);
			else add(t2, q[q[i].t].x + 1, -q[q[i].t].y), sum -= q[q[i].t].y;
		}
		ll p1 = 0, s1 = 0, s2 = 0;
		for(int j = 20;j >= 0;j --){
			int nxt = p1 + (1 << j);
			if(nxt <= tot && s1 + t1[nxt] < sum - (s2 + t2[nxt])){
				p1 = nxt, s1 += t1[nxt], s2 += t2[nxt];
			}
		}
		ll res1 = s1;
		ll res2 = 0, p2 = 0;
		if(p1 < tot){
			ll ss1 = 0, ss2 = 0;
			for(int j = p1 + 1;j;j -= lowbit(j)){
				ss1 += t1[j], ss2 += t2[j];
			}
			res2 = sum - ss2;
			if(res2 >= res1){
				p2 = 0, ss2 = 0;
				for(int j = 20;j >= 0;j --){
					int nxt = p2 + (1 << j);
					if(nxt <= tot && sum - (ss2 + t2[nxt]) >= res2){
						p2 = nxt;
						ss2 += t2[nxt];
					}
				}
			}
		}
		if(res1 == res2 && res1 == 0){
			cout << "Peace\n";
		}
		else if(res1 <= res2){
			cout << lsh[p2] << ' ' << res2 * 2 << '\n';
		}
		else{
			cout << lsh[p1] << ' ' << res1 * 2 << '\n';
		}
	}

	return 0;
}

这集神了。

先考虑一条链怎么做，即求 $\sum_{i=0}^k a_{x+i}\text{ xor } i$。因为先异或后求和，可以考虑拆位来算，不难发现，$\forall j\in [x,x+2^p)$，$a_j$ 的第 $p$ 位都异或上了 $0$，如果我们将第 $p$ 位去掉并只考虑小于 $p$ 位的贡献，不难发现区间分成了 $[x,x+2^{p-1}),[x+2^{p-1},x+2^p)$ 两个子区间，且他们要解决的问题结构完全相同，这启发我们用倍增去解决问题。

设 $f_{i,j}$ 表示只考虑小于 $j$ 位的值对于 $\sum_{p=0}^{2^j-1} a_{i+p}\text{ xor } p$，设 $g_{i,j,0/1}$ 表示 $1\sim i$ 中第 $j$ 位为 $0/1$ 的位置个数，则有：
$$f_{i,j}=f_{i,j-1}+f_{i+2^{j-1},j-1}+(g_{i+2^{j-1}-1,j-1,1}-g_{i-1,j-1,1}+g_{i+2^j-1,j-1,0}-g_{i+2^{j-1},j-1,0})2^{j-1}$$

这个转移式子大概意思是：对于一个长度为 $2^j$ 的区间，先将这个东西小于 $j-1$ 位的贡献求出来，然后考虑第 $j-1$ 位的贡献，前一半区间只有第 $j-1$ 位为 $0$ 能贡献，后一半只有为 $1$ 才能贡献。

然后我们回到询问，从高到低位拆位算贡献，同时倍增跳，能跳就跳，不能跳就将当前位置到终点部分剩余位置在当前拆的这一位上的贡献算出来。

for(int j=30;j>=0;j--){
	if((k>>j)&1){
		nxt=u+(1<<j);ans+=f[j][u];
		c1=g[1][j][nxt-1]-g[1][j][u-1];
		c0=g[0][j][r]-g[0][j][nxt-1];
		ans+=(c1+c0)*(1ll<<j);u=nxt;
	}else{
		ans+=(g[1][j][r]-g[1][j][u-1])*(1ll<<j);
	}
}

具体的，我们将倍增数组 $f_{i,j}$ 的定义改为从 $i$ 到 $i$ 所在长链往下 $2^j$ 步。对于轻重儿子合并，可以直接将轻儿子的倍增数组加到重儿子的倍增数组上，然后直接对当前点的 $\log$ 个位置倍增。对于 $g$ 也这样更改定义，但是从前缀和变成了后缀和，这样可能会出现很多问题：如倍增数组空位无法合并，后缀和可能访问到不该访问的值，但总而言之是小问题，模仿上面的倍增过程跳下去即可。

时间复杂度为 $O(n\log V)$。

真（ru）的（yi）是（ke）签（ai）到（miao）。

看到 $n$ 大得吓人，于是我们不乱想，考虑贪心。具体来说就是我们发现每次从根结点往下走一定是走到一个叶子的地方，具体来说就是我们可以把树上的一堆边看作一条条链之后贪心的选。

具体来说就是对于一个结点记录以它为起点往下走到叶子能够得到的最优价值，所以我们可以遍历每个结点的时候可以更新其父结点能够得到的最优价值，因为当前父结点已经确定了除去走向当前结点的那条路径的最优价值，我们只需要看一下走向当前结点的是不是更优就可以了，注意要是我们成功更新了，就需要把原本的那条价值较高的链和父结点断开，因为我们到重复的点是不重复计算的，没有重复更新就断开当前结点和父亲的连边。

代码实现就是可以把当前结点确定的路径全部都设置成它原本的价值，然后遍历结点更新答案，每次确定了一条链的时候（假如我们知道了当前父节点连哪一条与儿子的边最大或者确定了某一个结点与它父亲的连边要断开就能确定）就直接把它丢进一个数组里，最后把这个数组排序取前 $k$ 大就做完了，因为我们最多访问 $k$ 个叶子。代码很好写，因为观察父节点的生成方式，$i$ 的父节点一定是小于 $i$ 的，那我们甚至连递归都不用了。

[COGS 4161] hope I can sort 题目解析

 引入

 注意到本题的01序列的排序过程可简化为错位数 k（前 c0 个位置中 1 的个数）的链。

 因此可以考虑动态规划。

 分析

#include<bits/stdc++.h>
using namespace std;

const int MOD=998244353,N=5005;
int a[N],dp[N],nd[N];

int qp(int a,int b){
	int r=1;
	while(b){
		if(b&1)r=1ll*r*a%MOD;
		a=1ll*a*a%MOD;
		b>>=1;
	}
	return r;
}

int main(){
    freopen("hopeicansort.in","r",stdin);
    freopen("hopeicansort.out","w",stdout); 
	int n,m;
	cin>>n>>m;
	int c0=0;
	for(int i=1;i<=n;++i)
    {
		cin>>a[i];
		if(!a[i]) ++c0;
	}
	int k0=0;
	for(int i=1;i<=c0;++i) 
        if(a[i]) 
            ++k0;
	int mx=k0<c0?k0:n-c0;
	if(mx<k0) mx=k0;
	int T=1ll*n*(n-1)/2%MOD,invT=qp(T,MOD-2);
	dp[k0]=1;
	for(int t=0;t<m;++t)
    {
		for(int k=0;k<=mx;++k) nd[k]=0;
		nd[0]=(dp[0]+1ll*dp[1]*invT)%MOD;
		for(int k=1;k<=mx;++k)
        {
			int d=1ll*k*k%MOD*invT%MOD;
			nd[k]=(1ll*dp[k]*(1-d+MOD)+1ll*dp[k+1]*(k+1)%MOD*(k+1)%MOD*invT)%MOD;
		}
		for(int k=0;k<=mx;++k) dp[k]=nd[k];
	}
	cout<<dp[0];
	return 0;
}