P5459 [BJOI2016]回转寿司 (动态开点的权值线段树)

题目链接:点击这里

题目大意:
给定一个长度为 n n n 的序列 ,现在要从中选出一段连续子序列 [ l , r ] [l,r] [l,r] ,使得 L ≤ ∑ i = l r a i ≤ R L \le \sum^r_{i = l}a_i \le R Li=lraiR ,求符合条件的子序列数量

题目分析:
我们可以在 1 − n 1-n 1n 范围内对 r r r 进行枚举, l o g n logn logn 的复杂度求解此时符合条件的 l l l 的数目,求和即可
s u m x = ∑ i = 1 x a i sum_x=\sum^x_{i=1}a_i sumx=i=1xai ,故 ∑ i = l r a i = s u m r − s u m l − 1 \sum^r_{i=l}a_i=sum_r-sum_{l-1} i=lrai=sumrsuml1,故式子 L ≤ ∑ i = l r a i ≤ R L \le \sum^r_{i = l}a_i \le R Li=lraiR 可化为 L ≤ s u m r − s u m l − 1 ≤ R L \le sum_r-sum_{l-1} \le R Lsumrsuml1R 可继续化为 s u m r − R ≤ s u m l − 1 ≤ s u m r − L sum_r-R \le sum_{l-1} \le sum_r-L sumrRsuml1sumrL ,此时问题就转换成了在 [ s u m r − R , s u m r − L ] [sum_r-R,sum_r-L] [sumrR,sumrL] 中找符合条件的 s u m l − 1   ( l ∈ [ 1 , r ] ) sum_{l-1} \ (l \in [1,r]) suml1 (l[1,r]) 的数目
而此时还有一个问题就是这 L , R L,R L,R很大,并且没办法离散化(也可能是我太菜了),所以还要用动态开点建立线段树

具体细节见代码:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#define ll long long
#define inf 0x3f3f3f3f
#define Inf 0x3f3f3f3f3f3f3f3f
using namespace std;
int read()
{
	int res = 0,flag = 1;
	char ch = getchar();
	while(ch<'0' || ch>'9')
	{
		if(ch == '-') flag = -1;
		ch = getchar();
	}
	while(ch>='0' && ch<='9')
	{
		res = (res<<3)+(res<<1)+(ch^48);//res*10+ch-'0';
		ch = getchar();
	}
	return res*flag;
}
const int maxn = 1e5+5;
const int mod = 1e9+7;
const double pi = acos(-1);
const double eps = 1e-8;
struct segmenttree{
	int root,cnt;
	int lson[maxn<<8],rson[maxn<<8],sum[maxn<<8];
	void pushup(int now)
	{
		sum[now] = sum[lson[now]]+sum[rson[now]]; 
	}
	void update(int &now,ll pos,ll val,ll l = -Inf,ll r = Inf)
	{
		if(!now) now = ++cnt;
		if(l == r)
		{
			sum[now] += val;
			return ;
		}
		ll mid = (l+r)>>1;
		if(mid >= pos) update(lson[now],pos,val,l,mid);
		else update(rson[now],pos,val,mid+1,r);
		pushup(now);
	}
	ll query(int &now,ll ql,ll qr,ll l = -Inf,ll r = Inf)
	{
		if(!now) now = ++cnt;
		if(ql > r || qr < l) return 0;
		if(ql <= l && qr >= r) return sum[now];
		ll mid = (l+r)>>1;
		return query(lson[now],ql,qr,l,mid)+query(rson[now],ql,qr,mid+1,r);
	}
}tree;
int n,l,r;
ll ans,sum[maxn];
int main()
{
	n = read(),l = read(),r = read();
	for(int i = 1;i <= n;i++)
	{
		int x = read();
		sum[i] = sum[i-1]+x;
	}
	tree.update(tree.root,sum[0],1); //插入sum[0] 
	for(int i = 1;i <= n;i++) //枚举右端点 
	{
		ans += tree.query(tree.root,sum[i]-r,sum[i]-l);
		tree.update(tree.root,sum[i],1); //插入sum[i] 
	}
	printf("%lld\n",ans);
	return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值