51nod找朋友 1463

题意:

1463 找朋友
1.5 秒 262,144.0 KB 80 分 5级题
给定:
两个长度为n的数列A 、B
一个有m个元素的集合K
询问Q次
每次询问[l,r],输出区间内满足|Bi-Bj|∈K 的最大Ai+Aj

数据约定:
n,Q<=100000
m <= 10
0<=A[i]<=1000000000
1<=B[i]<=n
1<=K[i]<=n
保证B[i]互不相等

输入
n Q m
A1 A2 …An
B1 B2 …Bn
K1 K2 …Km
l1 r1
l2 r2
.
.
lQ rQ
输出
Q行,每行一个整数表示相对应的答案。
如果找不到这样的两个数则输出0。
输入样例
4 2 2
1 2 3 4
3 2 1 4
1 3
1 4
2 3
输出样例
7
5

思路:

(1)首先,我们先搞清楚题意哈!!给了一个长度为n的A数组,给了一个长度为n的B数组,和一个长度为m的数组,有q次询问哦!!每次询问给出[l , r]的区间,在这个区间里,我们每次都要找到满足|Bi-Bj|∈K 的最大的Ai+Aj,问的是最大的Ai+Aj的值是多少?

(2)扫描线,什么是扫描线呢?想一想求多个矩形的周长,我们是把每一条平行于x轴的边都当作一条扫描线,将扫描线从小到大排序,然后就从下到上进行扫描

(3)在这道题里也是运用了扫描线的思想,我们把每个区间的r从小到大进行排序,然后从左到右进行扫描,求出每一段的每个点可以更新的最大值。。然后查询的时候就是查询区间[G[i].l,G[i].r]的最大值

(4)怎么利用线段树呢?线段树的每一个结点代表这个结点所代表的区间的最大值。

(5)最大值怎么找呢?我们每次找B[i] + k或B[i] - k在B数组中的位置pos,但是注意⚠️,pos < j哦!!,那么此时,A[i] + A[j] = A[i] + A[pos],然后更新pos位置的最大值就好了

代码实现:

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn = 2e5 + 5;

int n,q,m;
int a[maxn];
int b[maxn];
int id[maxn];
int k[50];
struct node{
	int l;
	int r;
	int id;
}G[maxn];
int sum[maxn<<2];
int maxx[maxn];
int ans[maxn];

bool cmp(node a,node b){
	return a.r < b.r;
}

void update(int pos,int v,int l,int r,int k){
	if(l == r){
		sum[k] = max(sum[k],v);
		return ;
	}
	int mid = (l + r) >> 1;
	if(pos <= mid) update(pos,v,l,mid,k << 1);
	else update(pos,v,mid + 1,r,k << 1|1);
	sum[k] = max(sum[k << 1],sum[k << 1|1]);
}

int query(int L,int R,int l,int r,int k){
	if(L <= l&&r <= R){
		return sum[k];
	}
	int mid = (l + r) >> 1;
	int res = 0;
	if(R <= mid) res = max(query(L,R,l,mid,k<<1),res);
	else if(L > mid) res = max(query(L,R,mid+1,r,k << 1|1),res);
	else{
		res = max(res,query(L,R,l,mid,k<<1));
		res = max(res,query(L,R,mid + 1,r,k<<1|1));
	}
	return res;
}

int main(){
	scanf("%d%d%d",&n,&q,&m);
	for(int i = 1;i <= n;i++){
		scanf("%d",&a[i]);
	}
	for(int i = 1;i <= n;i++){
		scanf("%d",&b[i]);
		id[b[i]] = i;
	}
	for(int i = 1;i <= m;i++){
		scanf("%d",&k[i]);
	}
	for(int i = 1;i <= q;i++){
		scanf("%d%d",&G[i].l,&G[i].r);
		G[i].id = i;
	}
	sort(G + 1,G + 1 + q,cmp);
	int l = 1;
	for(int i = 1;i <= q;i++){
		int r = G[i].r;
		for(int j = l;j <= r;j++){
			for(int s = 1;s <= m;s++){
				int num = b[j] + k[s];
				int pos = id[num];
				if(num <= n&&pos < j&&a[j] + a[pos] > maxx[pos]){//找pos < j的数
					maxx[pos] = a[j] + a[pos];
					update(pos,maxx[pos],1,n,1);
				}
				num = b[j] - k[s];
				pos = id[num];
				if(num >= 1&&pos < j&&a[j] + a[pos] > maxx[pos]){//找pos < j的数
					maxx[pos] = a[j] + a[pos];
					update(pos,maxx[pos],1,n,1);
				}
			}
		}
		ans[G[i].id] = query(G[i].l,G[i].r,1,n,1);
		l = G[i].r;
	}
	for(int i = 1;i <= q;i++){
		printf("%d\n",ans[i]);
	}
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值