2022杭电3 K - Taxi （曼哈顿距离+二分）

input:

1
3 4
1 5 7
5 1 6
2 3 9
1 5
2 2
4 3
10 10

output:

6
4
5
9

题意：多组样例，对于每组样例，先给出一个n和m，n代表点的个数，m代表询问的个数，接下来n行，每行3个数(xi,yi,wi)，分别代表第i个点的坐标和权值，对于每组询问，首先给出一个坐标，让我们求出这个点到n个点中的值的最大值，这个点到第i个点的值定义为两点曼哈顿距离和i点权值的较小值。
解题思路：

如果没有w的限制，那么是经典问题。根据|x|=max(x,−x)，有max{|x′−xi|+|y′−yi|}=max{max(x′−xi,−x′+xi)+max(y′−yi,−y′+yi)}=max{x′−xi+y′−yi,−x′+xi+y′−yi,x′−xi−y′+yi,−x′+xi−y′+yi)}=max{(x′+y′)+(−xi−yi),(x′−y′)+(−xi+yi),(−x′+y′)+(xi−yi),(−x′−y′)+(xi+yi)}分别记录xi+yi、xi−yi、−xi+yi、−xi−yi的最大值即可在O(1)时间内求出所有点到(x′,y′)的曼哈顿距离的最大值。现在考虑加入w的限制。将所有城镇按照w从小到大排序，并记录排序后每个后缀的xi+yi、xi−yi、−xi+yi、−xi−yi的最大值，用于O(1)求给定点(x′,y′)到该后缀中所有点的距离最大值。选取按w排序后的第k个城镇，O(1)求出给定点(x′,y′)到第k..n个城镇的距离最大值d，有两种情况：

•wk<d，那么第k..n个城镇对答案的贡献至少为wk。用wk更新答案后，由于第1..k个城镇的w值均不超过wk，因此它们不可能接着更新答案，考虑范围缩小至[k+1,n]。

•wk≥d，那么第k..n个城镇对答案的贡献为d。用d更新答案后，考虑范围缩小至[1,k−1]。容易发现每次考虑的范围都是一个区间，如果每次取k为区间的中点，那么迭代O(logn)次即可得到最优解。时间复杂度O((n+q)logn)。

上代码：

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
using namespace std;
const int N=1E5+10;
int x[N],y[N],w[N],maxv[N][4];
int n,q;
struct node{
	int x;
	int y ;
	int w;
}town[N];
bool cmp(struct node a,struct node b)
{
	return a.w<b.w;
}
int main()
{
	std::ios::sync_with_stdio(false);
	cin.tie(0); 
	cout.tie(0);
	int t;
	cin>>t;
	while(t--)
	{
		memset(maxv,-0x3f3f3f3f,sizeof maxv);
		cin>>n>>q;
		for(int i=1;i<=n;i++)
			cin>>town[i].x>>town[i].y>>town[i].w;
		sort(town+1,town+n+1,cmp);
		for(int i=n;i>=1;i--)
		{
			maxv[i][0]=max(maxv[i+1][0],town[i].x+town[i].y);
			maxv[i][1]=max(maxv[i+1][1],town[i].x-town[i].y);
			maxv[i][2]=max(maxv[i+1][2],-town[i].x+town[i].y);
			maxv[i][3]=max(maxv[i+1][3],-town[i].x-town[i].y);
		}
		for(int i=1;i<=q;i++)
		{
			int x,y;
			cin>>x>>y;
			int l=1,r=n;
			int ans=-0x3f3f3f3f;
			while(l<r)
			{
				int mid=l+r>>1;
				int d=max(max(-x-y+maxv[mid][0],-x+y+maxv[mid][1]),max(x-y+maxv[mid][2],x+y+maxv[mid][3]));
				if(d>town[mid].w)
				{
					ans=max(ans,town[mid].w);
					l=mid+1;
				}
				else
				{
					ans=max(ans,d);
					r=mid-1;
				}
			}	
			int d=max(max(-x-y+maxv[l][0],-x+y+maxv[l][1]),max(x-y+maxv[l][2],x+y+maxv[l][3]));
			d=min(d,town[l].w);
			ans=max(ans,d);
			cout<<ans<<endl;
		}
	}
	return 0;
 }