《GMOJ-Junior-2470 【GDKOI2021 普及组 Day1】配对》题解

题目大意

给出一个区间$[l,r]$和$n$个非负整数$a_1,a_2,\cdots ,a_n$，问最少要删去多少个数，才能使剩下的数两两配对，且每一对数的和在$[l,r]$中。
对于$20\%$的数据，$n\le 10$；
对于$30\%$的数据，$n\le 20$；
对于$50\%$的数据，$n\le 100$；
对于$70\%$的数据，$n\le 10^3$；
对于$100\%$的数据，$n\le 10^6$，$a_i\le 10^9$。

分析

这是一道用到贪心算法的题。在解题之前，我们需要先寻找一些结论。
我们考虑把这些数从小到大排序，然后考虑给这些数分组。因为排完序的数组满足$a_1\le a_2\le \cdots \le a_n$，所以当我们发现有一个数对$(a_x,a_y)$满足$a_x+a_y>r$，那么$(a_x,a_{y+1}),(a_x,a_{y+2}),\cdots ,(a_x,a_n)$都满足$数对和>r$；当我们发现有一个数对$(a_x,a_y)$满足$a_x+a_y<l$，那么$(a_{x-1},a_y),(a_{x-2},a_y),\cdots ,(a_1,a_y)$都满足$数对和<l$。同时，我们可以发现当$p_1<p_2<p_3<p_4$时，若 { a p 1 + a p 3 , a p 2 + a p 4 } ∈ [ l , r ] \{ a_{p_1}+a_{p_3},a_{p_2}+a_{p_4} \} \in [l,r] {ap1​​+ap3​​,ap2​​+ap4​​}∈[l,r]，那么有$l\le a_{p_1}+a_{p_3}\le a_{p_2}+a_{p_3},a_{p_1}+a_{p_4}\le a_{p_2}+a_{p_4}\le r$，所以有 { a p 1 + a p 4 , a p 2 + a p 3 } ∈ [ l , r ] \{ a_{p_1}+a_{p_4},a_{p_2}+a_{p_3} \} \in [l,r] {ap1​​+ap4​​,ap2​​+ap3​​}∈[l,r]。
由此，我们可以在给$a_1,a_2,\cdots ,a_n$排序后，用两个指针$x,y$求出答案。初始时令$x=1,y=n$，然后我们考虑向右移动$x$，向左移动$y$。我们在处理的过程中会遇到如下几种情况：

1. 当我们发现$x>y$时，我们可以知道此时我们处理完了所有数。此时我们直接输出答案，退出程序即可。
2. 当我们当我们发现$x=y$时，我们可以知道此时我们处理剩了一个数。此时我们输出$答案+1$，然后退出程序。
3. 当我们发现$a_x+a_y\in [l,r]$时，我们发现我们找到了一个合法的数对。此时我们让$x$加$1$，让$y$减$1$。
4. 当我们发现$a_x+a_y>r$时，我们可以确定$a_y$一定是要被舍去的数。因为$a_x+a_y>r$时，$a_{y+1},a_{y+2},\cdots ,a_n$都无法与$a_x$配对，且不存在$x^{\prime }<x,y^{\prime }>y$，使得 { a x ′ + a y , a x + a y ′ } ∈ [ l , r ] \{ a_{x'}+a_{y},a_{x}+a_{y'} \} \in [l,r] {ax′​+ay​,ax​+ay′​}∈[l,r]。这时，我们把$a_y$删去，将答案加$1$，并使$y$减$1$。
5. 当我们发现$a_x+a_y<l$时，我们可以确定$a_x$一定是要被舍去的数。与情况$4$类似，因为$a_x+a_y<l$时，$a_{x-1},a_{x-2},\cdots ,a_1$都无法与$a_y$配对，且不存在$x^{\prime }>x,y^{\prime }<y$，使 { a x ′ + a y , a x + a y ′ } ∈ [ l , r ] \{ a_{x'}+a_{y},a_{x}+a_{y'} \} \in [l,r] {ax′​+ay​,ax​+ay′​}∈[l,r]。这时，我们把$a_x$删去，将答案加$1$，并使$x$加$1$。
   这样，我们就可以在$O(n\log n)$的时间复杂度内解决本题了。

代码

根据我们的思路，可以写出如下代码：

#include<cstdio>
char buf[16777216];
inline int Read() //快速读入
{
	static char* c = buf;
	while (*c < '0' || *c>'9')
	{
		++c;
	}
	int ans = *c ^ 48;
	while (*(++c) >= '0' && *c <= '9')
	{
		ans = ans * 10 + (*c ^ 48);
	}
	return ans;
}
int a[1000001];
inline void QuickSort(const int l, const int r) //快速排序
{
	const int standard = (a[l] + a[r]) >> 1;
	int i = l, j = r;
	while (i <= j)
	{
		while (a[i] < standard)
		{
			++i;
		}
		while (a[j] > standard)
		{
			--j;
		}
		if (i <= j)
		{
			const int t = a[i];
			a[i] = a[j];
			a[j] = t;
			++i;
			--j;
		}
	}
	if (i < r)
	{
		QuickSort(i, r);
	}
	if (j > l)
	{
		QuickSort(l, j);
	}
	return;
}
int main()
{
	static_cast<void>(freopen("match.in", "r", stdin)); //定义文件输入输出
	static_cast<void>(freopen("match.out", "w", stdout));
	fread(buf, 1, 16777216, stdin); //读入数据
	const int n = Read(), l = Read(), r = Read(); //读入n,l,r
	for (int i = 1; i <= n; ++i) //读入a
	{
		a[i] = Read();
	}
	QuickSort(1, n); //对a排序
	int ans = 0;
	for (int i = 1, j = n; i <= j;) //统计答案，其中i对应分析中的x，j对应分析中的y
	{
		if (i == j) //情况2
		{
			++ans;
			break;
		}
		while (i < j && (a[i] + a[j] >= l && a[i] + a[j] <= r)) //情况3
		{
			++i;
			--j;
		}
		while (i < j && a[i] + a[j] > r) //情况4
		{
			++ans;
			--j;
		}
		while (i < j && a[i] + a[j] < l) //情况5
		{
			++ans;
			++i;
		}
	}
	printf("%d", ans);
	return 0;
}

总结

这是一道需要一定思维能力的题，需要在分析后解决。