【题解】AT4520:[AGC032E] Modulo Pairing

原题传送门
首先想到的当然是二分答案+贪心
但是正解二分的并不是答案

结论:最终的正确答案一定存在一个分界点,满足
在这里插入图片描述
蓝线表示 x + y < M x+y<M x+y<M
红线表示 x + y > = M x+y>=M x+y>=M

给出证明

情况1:
在这里插入图片描述

任意两个数加起来都 < M <M <M的情况,Ⅰ最优
令从左到右分别为 a < b < c < d a<b<c<d a<b<c<d
计算出三种方案的价值
Ⅰ: m a x ( a + d , b + c ) max(a+d,b+c) max(a+d,b+c)
Ⅱ: m a x ( a + b , c + d ) = c + d > a + d , b + c , 即 Ⅱ > Ⅰ max(a+b,c+d)=c+d>a+d,b+c,即Ⅱ>Ⅰ max(a+b,c+d)=c+d>a+d,b+c,>
Ⅲ: m a x ( a + c , b + d ) = b + d > a + d , b + c , 即 Ⅲ > Ⅰ max(a+c,b+d)=b+d>a+d,b+c,即Ⅲ>Ⅰ max(a+c,b+d)=b+d>a+d,b+c,>
证毕

情况2:
在这里插入图片描述任意两个数加起来都 > = M >=M >=M的情况,Ⅰ最优
在证明之前先要发现,因为每个数 < M <M <M
所以对于任意 x + y > = M , x + y − M < m i n ( x , y ) x+y>=M,x+y-M<min(x,y) x+y>=M,x+yM<min(x,y)

Ⅰ: m a x ( a + d , b + c ) − M max(a+d,b+c)-M max(a+d,b+c)M
Ⅱ: m a x ( a + b , c + d ) − M max(a+b,c+d)-M max(a+b,c+d)M
Ⅲ: m a x ( a + c , b + d ) − M max(a+c,b+d)-M max(a+c,b+d)M
其实变成了情况1的证明

情况3:
在这里插入图片描述
Ⅰ: m a x ( a + d − M , b + c ) = b + c max(a+d-M,b+c)=b+c max(a+dM,b+c)=b+c
Ⅱ: m a x ( a + b , c + d − M ) = m a x ( a + b , e ) ( e 为 < c 的 一 个 数 ) max(a+b,c+d-M)=max(a+b,e)(e为<c的一个数) max(a+b,c+dM)=max(a+b,e)(e<c)
< < <

情况4:
在这里插入图片描述
Ⅰ: m a x ( b + d − M , a + c ) = a + c max(b+d-M,a+c)=a+c max(b+dM,a+c)=a+c
Ⅱ: m a x ( a + b , c + d − M ) = m a x ( a + b , e ) ( e 为 < c 的 一 个 数 ) max(a+b,c+d-M)=max(a+b,e)(e为<c的一个数) max(a+b,c+dM)=max(a+b,e)(e<c)
< < <

情况5:
在这里插入图片描述
Ⅰ: m a x ( b + c − M , a + d ) = a + d max(b+c-M,a+d)=a+d max(b+cM,a+d)=a+d
Ⅱ: m a x ( a + b , c + d − M ) = m a x ( a + b , e ) ( e 为 < c 的 一 个 数 ) max(a+b,c+d-M)=max(a+b,e)(e为<c的一个数) max(a+b,c+dM)=max(a+b,e)(e<c)
< < <

所有情况都罗列了
结论就证明好了

所以正确答案一定满足最开始说的那种情况,分界点左边全是蓝线,最大最小相连,次大次小相连……,右边全是红线,最大最小相连,次大次小相连……
然后分界点是越小越好,合法的分界点越小,左右两边的权值越小,同时的左边越可能满足蓝线
所以答案就一定是合法的最小的分界点

二分分界点,找到最小的就行了

Code:

#include <bits/stdc++.h>
#define maxn 200010
using namespace std;
int a[maxn], n, m;

inline int read(){
	int s = 0, w = 1;
	char c = getchar();
	for (; !isdigit(c); c = getchar()) if (c == '-') w = -1;
	for (; isdigit(c); c = getchar()) s = (s << 1) + (s << 3) + (c ^ 48);
	return s * w;
}

bool check(int mid){
	for (int i = mid + 1, j = n; i <= j; ++i, --j) if (a[i] + a[j] < m) return 0;
	return 1;
}

int main(){
	n = read() << 1, m = read();
	for (int i = 1; i <= n; ++i) a[i] = read();
	sort(a + 1, a + 1 + n);
	int l = 0, r = n >> 1, ans = 0;
	while (l <= r){
		int mid = (l + r) >> 1;
		if (check(mid << 1)) ans = mid << 1, r = mid - 1; else l = mid + 1;
	}
	int sum = 0;
	for (int i = 1, j = ans; i <= j; ++i, --j) sum = max(sum, a[i] + a[j]);
	for (int i = ans + 1, j = n; i <= j; ++i, --j) sum = max(sum, a[i] + a[j] - m);
	printf("%d\n", sum);
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值