[ROI2019]课桌

前言

在 $2021/10/15$ 这一天的上午，钠悴 $\mathsf{D}\mathsf{D}\mathsf{G}$ 闪击 $\mathsf{O}\mathsf{n}\mathsf{e}\mathsf{I}\mathsf{n}\mathsf{D}\mathsf{a}\mathsf{r}\mathsf{k}$，在毫无准备的情况下，$\mathsf{O}\mathsf{n}\mathsf{e}\mathsf{I}\mathsf{n}\mathsf{D}\mathsf{a}\mathsf{r}\mathsf{k}$ 瞬间就被「随切」秒杀了！

随后 $\mathsf{D}\mathsf{D}\mathsf{G}$ 一路高歌猛进，在 $39\mathrm{s}$ 内拿下了号称 “地表最硬马桶” 的 $\mathsf{H}\mathsf{a}\mathsf{n}\mathsf{d}\mathsf{I}\mathsf{n}\mathsf{D}\mathsf{e}\mathsf{v}\mathsf{i}\mathsf{l}$ 。

接下来 $\mathsf{D}\mathsf{D}\mathsf{G}$ 又对 $\mathrm{p}\mathrm{r}\mathrm{i}\mathrm{n}\mathrm{c}\mathrm{i}\mathrm{p}\mathrm{a}\mathrm{l}$ 发动突然袭击，护校战争爆发！在全校同学的顽强抵抗下，取得了首都 $\mathrm{D}\mathrm{D}2023$ 保卫战的伟大胜利，成为了第二次反狗战争的转折点。

此时，$\mathsf{R}\mathsf{a}\mathsf{i}\mathsf{n}\mathsf{y}\mathsf{b}\mathsf{u}\mathsf{n}\mathsf{n}\mathsf{y}$ 成功登陆撤硕海滩，与 $\mathrm{p}\mathrm{r}\mathrm{i}\mathrm{n}\mathrm{c}\mathrm{i}\mathrm{p}\mathrm{a}\mathrm{l}$ 一同击退了 $\mathsf{D}\mathsf{D}\mathsf{G}$ 。第二次反狗战争结束了，但是对 $\mathsf{O}\mathsf{n}\mathsf{e}\mathsf{I}\mathsf{n}\mathsf{D}\mathsf{a}\mathsf{r}\mathsf{k}$ 带来的伤害是不可估量的，$\mathsf{O}\mathsf{n}\mathsf{e}\mathsf{I}\mathsf{n}\mathsf{D}\mathsf{a}\mathsf{r}\mathsf{k}$ 与全世界所有爱好和平的人类振臂高呼：危险的黑暗幽灵，不要再在人们头顶飘荡！

题目

传送门 to LOJ

思路

一上午全去推结论去了。结果旁边的所有碳基生物都说：“显而易见，一目了然。”

首先考虑把问题放在数轴上。——其实最初我在想是不是什么匹配，网络流啥的……然后又考虑 $\mathtt{d}\mathtt{p}$，发现不咋会……

那么此时 “不适度” 就是数轴上一个点移动到区间内走过的距离。相当于把 $2n$ 个点移动到别的位置，然后可以选择区间内两个点并移除，最终可以移除所有点。

第一个 $\mathrm{o}\mathrm{b}\mathrm{s}\mathrm{e}\mathrm{r}\mathrm{v}\mathrm{a}\mathrm{t}\mathrm{i}\mathrm{o}\mathrm{n}$ 是，点是没有区别的。所以两个点的移动路径是不会 “交叉” 的，也就是说，靠右的点所移动到的位置也更靠右。第二个 $\mathrm{o}\mathrm{b}\mathrm{s}\mathrm{e}\mathrm{r}\mathrm{v}\mathrm{a}\mathrm{t}\mathrm{i}\mathrm{o}\mathrm{n}$ 是，桌子对应的区间越大越好，也就是说桌子对应的区间不会有包含关系。

于是有了第一个性质：总是身高相邻的同学坐在同一张桌子旁。为什么呢？这就需要分析了。不妨设 $a⩽b⩽c⩽d$ 是四个同学的身高，四个人在数轴上走到的位置依次是 $x⩽y⩽z⩽w$ 。

• $a$ 和 $d$ 坐在 $i$ 号桌子旁，而 $b$ 和 $c$ 坐在 $j$ 号桌子旁。数学公式警告！

上面的条件等价于 $x,w\in [L_i,R_i]$ 而 $y,z\in [L_j,R_j]$ 。首先
$$y,z\in [x,w]⫅[L_i,R_i]$$
而 $y\in [x,w]$ 且 $y\in [L_j,R_j]$，说明 $[x,w]\cap [L_j,R_j]$ 非空。而 $[x,w]$ 不能包含 $[L_j,R_j]$，否则 $[L_i,R_i]⫆[x,w]⫆[L_j,R_j]$，与 observation 2 矛盾。所以
$$x\in [L_j,R_j]\vee w\in [L_j,R_j]$$
譬如 $x\in [L_j,R_j]$，那么 $x,y\in [L_j,R_j]$ 而 $z,w\in [L_i,R_i]$，即可以 $a,b$ 坐在桌子 $j$ 旁、$c,d$ 坐在桌子 $i$ 旁。

若是 $w\in [L_j,R_j]$，那么 $z,w\in [L_j,R_j]$ 而 $x,y\in [L_i,R_i]$，即可以 $a,b$ 坐在桌子 $i$ 旁、$c,d$ 坐在桌子 $j$ 旁。

• $a$ 和 $c$ 坐在 $i$ 号桌子旁，而 $b$ 和 $d$ 坐在 $j$ 号桌子旁。

这个比较简单，由 $y\in [x,z]⫅[L_i,R_i]$ 且 $z\in [y,w]⫅[L_j,R_j]$ 知可以 $a,b$ 坐 $i$ 号桌子、$c,d$ 坐在 $j$ 号桌子旁。

证毕，证毕，证毕。重要的事说三遍，不重要的事也可以说三遍。

这样可以解决 $m=1$ 的情况，但是又怎么拓展呢？于是考虑桌子的选用。事实上，依次选用左端点递增的桌子 就好了。它的证明也是类似的，而且更简单。

• $a,b$ 走到 $x,y\in [L_i,R_i]$，而 $c,d$ 走到 $z,w\in [L_j,R_j]$，满足 $a⩽b⩽c⩽d$ 然而 $L_i>L_j$ 。根据 observation 2 同时也满足 $R_i>R_j$ 。

由于 $z,w⩾y⩾L_i$ 且 $z,w⩽R_j<R_i$ 可知 $z,w\in [L_i,R_i]$ 。

由于 $x,y⩾L_i>L_j$ 且 $x,y⩽z⩽R_j$ 可知 $x,y\in [L_j,R_j]$ 。

所以直接把 $x,y$ 和 $z,w$ 所对应的桌子互换，并不增加代价。证毕。

这说明什么呢？假如你买好了桌子，那么 任意一组的第 $2i-1$ 名和第 $2i$ 名同学都是用的同一桌子。于是每个这样的 $i$ 是一个独立的问题，只需要求出最优的那张桌子。

而根据前面的描述，决策是单调的。二分即可。时间复杂度 $\mathcal{O}(n\log m+nm+k\log n+k\log k)$ 。为什么没有额外的 $\log$ 呢？因为我们已经用 $\mathcal{O}(k\log k)$ 将桌子排序过了，就可以直接用类似双指针的方式求值（而不是 $\mathrm{l}\mathrm{o}\mathrm{w}\mathrm{e}\mathrm{r}\_\mathrm{b}\mathrm{o}\mathrm{u}\mathrm{n}\mathrm{d}$ 现场查）。

代码

#include <cstdio>
#include <iostream>
#include <vector>
#include <cstring>
#include <algorithm>
using namespace std;
# define rep(i,a,b) for(int i=(a); i<=(b); ++i)
# define drep(i,a,b) for(int i=(a); i>=(b); --i)
typedef long long int_;
inline int readint(){
	int a = 0, c = getchar(), f = 1;
	for(; '0'>c||c>'9'; c=getchar())
		if(c == '-') f = -f;
	for(; '0'<=c&&c<='9'; c=getchar())
		a = (a<<3)+(a<<1)+(c^48);
	return a*f;
}
inline void writeint(int_ x){
	if(x > 9) writeint(x/10);
	putchar(int((x-x/10*10)^48));
}

const int MAXN = 200005;
struct Table{
	int l, r;
	bool operator < (const Table &t) const {
		if(l != t.l) return l < t.l;
		return r > t.r; // bigger
	}
};
Table table[MAXN];
vector<int> a[MAXN];
int tmp[MAXN<<1], n, m, k;

int_ ans;
void solve(int l,int r,int ql,int qr){
	if(l > r) return ; // empty
	int mid = (l+r)>>1;
	vector<int> &v = a[mid];
	sort(v.begin(),v.end());
	int lid = 0, rid = 0; int_ now_val = 0;
	while(lid != (m<<1) && v[lid] < table[ql].l)
		now_val += table[ql].l-v[lid], ++ lid;
	while(rid != (m<<1) && v[rid] < table[ql].r)
		now_val += table[ql].r-v[rid], ++ rid;
	rep(i,lid,(m<<1)-1) now_val += v[i]-table[ql].l;
	rep(i,rid,(m<<1)-1) now_val += v[i]-table[ql].r;
	int best = ql; int_ best_val = (now_val >>= 1);
	best_val -= m*int_(table[ql].r-table[ql].l);
	for(int i=ql+1; i<=qr; ++i){
		now_val += (lid-m)*int_(table[i].l-table[i-1].l);
		now_val += (rid-m)*int_(table[i].r-table[i-1].r);
		while(lid != (m<<1) && v[lid] < table[i].l)
			now_val += table[i].l-v[lid], ++ lid;
		while(rid != (m<<1) && v[rid] < table[i].r)
			now_val += table[i].r-v[rid], ++ rid;
		const int_ addend = m*int_(table[i].r-table[i].l);
		if(now_val-addend < best_val)
			best_val = now_val-addend, best = i;
	}
	ans += best_val; // current value
	solve(l,mid-1,ql,best);
	solve(mid+1,r,best,qr);
}

int main(){
	m = readint(), n = readint();
	k = readint();
	rep(i,1,k){
		table[i].l = readint();
		table[i].r = readint();
	}
	sort(table+1,table+k+1);
	for(int i=2,d=0; i<=k; ++i){
		table[i] = table[i+d];
		if(table[i].r <= table[i-1].r)
			++ d, -- k, -- i; // reload
	}
	rep(i,1,m){
		rep(j,1,n<<1) tmp[j] = readint();
		sort(tmp+1,tmp+(n<<1)+1);
		rep(j,1,n<<1) a[(j+1)>>1].push_back(tmp[j]);
	}
	solve(1,n,1,k);
	printf("%lld\n",ans);
	return 0;
}