BZOJ4627 权值线段树

4627: [BeiJing2016]回转寿司

Time Limit: 10 Sec Memory Limit: 256 MB
Submit: 1204 Solved: 475
[Submit][Status][Discuss]

Description

酷爱日料的小Z经常光顾学校东门外的回转寿司店。在这里，一盘盘寿司通过传送带依次呈现在小Z眼前。不同的寿

司带给小Z的味觉感受是不一样的，我们定义小Z对每盘寿司都有一个满意度，例如小Z酷爱三文鱼，他对一盘三文

鱼寿司的满意度为10；小Z觉得金枪鱼没有什么味道，他对一盘金枪鱼寿司的满意度只有5；小Z最近看了电影“美

人鱼”，被里面的八爪鱼恶心到了，所以他对一盘八爪鱼刺身的满意度是-100。特别地，小Z是个著名的吃货，他

吃回转寿司有一个习惯，我们称之为“狂吃不止”。具体地讲，当他吃掉传送带上的一盘寿司后，他会毫不犹豫地

吃掉它后面的寿司，直到他不想再吃寿司了为止。今天，小Z再次来到了这家回转寿司店，N盘寿司将依次经过他的

面前，其中，小Z对第i盘寿司的满意度为Ai。小Z可以选择从哪盘寿司开始吃，也可以选择吃到哪盘寿司为止，他

想知道共有多少种不同的选择，使得他的满意度之和不低于L，且不高于R。注意，虽然这是回转寿司，但是我们不

认为这是一个环上的问题，而是一条线上的问题。即，小Z能吃到的是输入序列的一个连续子序列；最后一盘转走

之后，第一盘并不会再出现一次。

Input

第一行包含三个整数N，L和R，分别表示寿司盘数，满意度的下限和上限。

第二行包含N个整数Ai，表示小Z对寿司的满意度。

N≤100000，|Ai|≤100000，0≤L, R≤10^9

Output

仅一行，包含一个整数，表示共有多少种选择可以使得小Z的满意度之和

不低于L且不高于R。

Sample Input

5 5 9
1 2 3 4 5

Sample Output

用权值线段树进行维护；

预处理出前缀和后，就是查询

sum[i]-R~sum[i]-L这个区间sum[j]的数量；

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstdlib>
#include<cstring>
#include<string>
#include<cmath>
#include<map>
#include<set>
#include<vector>
#include<queue>
#include<bitset>
#include<ctime>
#include<time.h>
#include<deque>
#include<stack>
#include<functional>
#include<sstream>
//#include<cctype>
//#pragma GCC optimize(2)
using namespace std;
#define maxn 200005
#define inf 10000000005ll
//#define INF 1e18
#define rdint(x) scanf("%d",&x)
#define rdllt(x) scanf("%lld",&x)
#define rdult(x) scanf("%lu",&x)
#define rdlf(x) scanf("%lf",&x)
#define rdstr(x) scanf("%s",x)
#define mclr(x,a) memset((x),a,sizeof(x))
typedef long long  ll;
typedef unsigned long long ull;
typedef unsigned int U;
#define ms(x) memset((x),0,sizeof(x))
const long long int mod = 1e9 + 7;
#define Mod 1000000000
#define sq(x) (x)*(x)
#define eps 1e-5
typedef pair<int, int> pii;
#define pi acos(-1.0)
//const int N = 1005;
#define REP(i,n) for(int i=0;i<(n);i++)
typedef pair<int, int> pii;inline int rd() {int x = 0;char c = getchar();bool f = false;while (!isdigit(c)) {if (c == '-') f = true;c = getchar();}while (isdigit(c)) {x = (x << 1) + (x << 3) + (c ^ 48);c = getchar();}return f ? -x : x;
}ll gcd(ll a, ll b) {return b == 0 ? a : gcd(b, a%b);
}
int sqr(int x) { return x * x; }/*ll ans;
ll exgcd(ll a, ll b, ll &x, ll &y) {if (!b) {x = 1; y = 0; return a;}ans = exgcd(b, a%b, x, y);ll t = x; x = y; y = t - a / b * y;return ans;
}
*/struct node {ll l, r, cnt;
}tr[maxn<<4];ll sum[maxn];
ll root = 1, cntrt = 1;
ll L, R;void Insert(ll &o, ll l, ll r, ll v) {if (!o)o = ++cntrt;tr[o].cnt++;if (l == r)return;ll mid = (l + r) >> 1;if (v <= mid)Insert(tr[o].l, l, mid, v);else Insert(tr[o].r, mid + 1, r, v);
}int query(ll l, ll r, ll LL, ll RR, ll rt) {if (LL <= l && r <= RR)return tr[rt].cnt;ll mid = (l + r) >> 1;int ans = 0;if (LL <= mid)ans += query(l, mid, LL, RR, tr[rt].l);if (mid < RR)ans += query(mid + 1, r, LL, RR, tr[rt].r);return ans;
}int main()
{// ios::sync_with_stdio(0);ll n; rdllt(n); rdllt(L); rdllt(R);for (int i = 1; i <= n; i++) {ll x; rdllt(x); sum[i] = sum[i - 1] + x;}Insert(root, -inf, inf, 0);ll ans = 0;for (int i = 1; i <= n; i++) {ans += 1ll * query(-inf, inf, sum[i] - R, sum[i] - L, root);Insert(root, -inf, inf, sum[i]);}printf("%lld\n", 1ll * ans);return 0;
}

转载于:https://www.cnblogs.com/zxyqzy/p/10381864.html