Yaroslav and Points

明明区间合并一下就好的东西, 为什么我会写得这么麻烦的方法啊啊啊。

#include<bits/stdc++.h>
#define LL long long
#define fi first
#define se second
#define mk make_pair
#define PLL pair<LL, LL>
#define PLI pair<LL, int>
#define PII pair<int, int>
#define SZ(x) ((int)x.size())
#define ull unsigned long long using namespace std; const int N = 2e5 + ;
const int inf = 0x3f3f3f3f;
const LL INF = 0x3f3f3f3f3f3f3f3f;
const int mod = 1e9 + ;
const double eps = 1e-;
const double PI = acos(-); vector<int> oo;
int getPos(int x) {
return lower_bound(oo.begin(), oo.end(), x) - oo.begin() + ;
}
int getPosL(int x) {
return lower_bound(oo.begin(), oo.end(), x) - oo.begin() + ;
}
int getPosR(int x) {
return upper_bound(oo.begin(), oo.end(), x) - oo.begin();
} namespace SGT {
#define lson l, mid, rt << 1
#define rson mid + 1, r, rt << 1 | 1
LL sum[N << ], cnt[N << ];
LL lazya[N << ], lazyb[N << ], all[N << ];
inline void push(int rt) {
if(lazya[rt]) {
lazya[rt << ] += lazya[rt];
lazya[rt << | ] += lazya[rt];
all[rt << ] += cnt[rt << ] * lazya[rt];
all[rt << | ] += cnt[rt << | ] * lazya[rt];
lazya[rt] = ;
}
if(lazyb[rt]) {
lazyb[rt << ] += lazyb[rt];
lazyb[rt << | ] += lazyb[rt];
all[rt << ] += sum[rt << ] * lazyb[rt];
all[rt << | ] += sum[rt << | ] * lazyb[rt];
lazyb[rt] = ;
}
}
inline void pull(int rt) {
sum[rt] = sum[rt << ] + sum[rt << | ];
cnt[rt] = cnt[rt << ] + cnt[rt << | ];
all[rt] = all[rt << ] + all[rt << | ];
}
void modify1(int p, LL val, int l, int r, int rt) {
if(l == r) {
cnt[rt] += val;
sum[rt] += oo[l - ] * val;
if(!cnt[rt]) all[rt] = ;
return;
}
push(rt);
int mid = l + r >> ;
if(p <= mid) modify1(p, val, lson);
else modify1(p, val, rson);
pull(rt);
}
void modify2(int L, int R, LL val, int l, int r, int rt) {
if(L > R) return;
if(l >= L && r <= R) {
lazya[rt] += val;
all[rt] += val * cnt[rt];
return;
}
push(rt);
int mid = l + r >> ;
if(L <= mid) modify2(L, R, val, lson);
if(R > mid) modify2(L, R, val, rson);
pull(rt);
}
void modify3(int L, int R, LL val, int l, int r, int rt) {
if(L > R) return;
if(l >= L && r <= R) {
lazyb[rt] += val;
all[rt] += val * sum[rt];
return;
}
push(rt);
int mid = l + r >> ;
if(L <= mid) modify3(L, R, val, lson);
if(R > mid) modify3(L, R, val, rson);
pull(rt);
}
PLL query1(int L, int R, int l, int r, int rt) {
if(L > R) return mk(, );
if(l >= L && r <= R) return mk(sum[rt], cnt[rt]);
push(rt);
int mid = l + r >> ;
PLL ans = mk(, ), tmp = mk(, );
if(L <= mid) {
tmp = query1(L, R, lson);
ans.fi += tmp.fi; ans.se += tmp.se;
}
if(R > mid) {
tmp = query1(L, R, rson);
ans.fi += tmp.fi; ans.se += tmp.se;
}
return ans;
}
LL query2(int L, int R, int l, int r, int rt) {
if(L > R) return ;
if(l >= L && r <= R) return all[rt];
push(rt);
int mid = l + r >> ;
LL ans = ;
if(L <= mid) ans += query2(L, R, lson);
if(R > mid) ans += query2(L, R, rson);
return ans;
}
} int n, m, x[N], c[N];
pair<PII, int> qus[N]; int main() {
scanf("%d", &n);
for(int i = ; i <= n; i++) {
scanf("%d", &x[i]);
c[i] = x[i];
oo.push_back(x[i]);
}
scanf("%d", &m);
for(int i = ; i <= m; i++) {
scanf("%d%d%d", &qus[i].se, &qus[i].fi.fi, &qus[i].fi.se);
if(qus[i].se == ) {
c[qus[i].fi.fi] += qus[i].fi.se;
oo.push_back(c[qus[i].fi.fi]);
}
}
sort(oo.begin(), oo.end());
oo.erase(unique(oo.begin(), oo.end()), oo.end());
for(int i = ; i <= n; i++) SGT::modify1(getPos(x[i]), , , SZ(oo), );
for(int i = ; i <= n; i++) {
PLL add = SGT::query1(getPos(x[i]), SZ(oo), , SZ(oo), );
SGT::modify2(getPos(x[i]), getPos(x[i]), add.fi, , SZ(oo), );
SGT::modify3(getPos(x[i]), getPos(x[i]), -add.se, , SZ(oo), );
}
for(int i = ; i <= m; i++) {
if(qus[i].se == ) {
int p = qus[i].fi.fi, d = qus[i].fi.se;
if(d > ) {
SGT::modify1(getPos(x[p]), -, , SZ(oo), );
SGT::modify2(, getPos(x[p]) - , d, , SZ(oo), );
SGT::modify2(getPos(x[p]) + , getPos(x[p] + d) - , x[p] + d, , SZ(oo), );
SGT::modify3(getPos(x[p]) + , getPos(x[p] + d) - , -, , SZ(oo), );
x[p] += d;
SGT::modify1(getPos(x[p]), , , SZ(oo), );
PLL add = SGT::query1(getPos(x[p]), SZ(oo), , SZ(oo), );
SGT::modify2(getPos(x[p]), getPos(x[p]), add.fi, , SZ(oo), );
SGT::modify3(getPos(x[p]), getPos(x[p]), -add.se, , SZ(oo), );
} else if(d < ) {
d = -d;
SGT::modify1(getPos(x[p]), -, , SZ(oo), );
SGT::modify2(, getPos(x[p] - d) - , -d, , SZ(oo), );
SGT::modify2(getPos(x[p] - d) + , getPos(x[p]) - , -x[p], , SZ(oo), );
SGT::modify3(getPos(x[p] - d) + , getPos(x[p]) - , , , SZ(oo), );
x[p] -= d;
SGT::modify1(getPos(x[p]), , , SZ(oo), );
PLL add = SGT::query1(getPos(x[p]), SZ(oo), , SZ(oo), );
SGT::modify2(getPos(x[p]), getPos(x[p]), add.fi, , SZ(oo), );
SGT::modify3(getPos(x[p]), getPos(x[p]), -add.se, , SZ(oo), );
}
} else {
int L = qus[i].fi.fi, R = qus[i].fi.se;
L = getPosL(L); R = getPosR(R);
LL ans = SGT::query2(L, R, , SZ(oo), );
PLL tmpL = SGT::query1(L, R, , SZ(oo), );
PLL tmpR = SGT::query1(R + , SZ(oo), , SZ(oo), );
ans -= tmpL.se * tmpR.fi;
ans += tmpL.fi * tmpR.se;
printf("%lld\n", ans);
}
}
return ;
} /*
*/

Codeforces 295E Yaroslav and Points 线段树的更多相关文章

  1. Codeforces Beta Round #19D(Points)线段树

    D. Points time limit per test 2 seconds memory limit per test 256 megabytes input standard input out ...

  2. codeforces 295E Yaroslav and Points (离线操作+离散化+区间合并)

    参考链接:http://blog.csdn.net/dyx404514/article/details/8817717 写的很详细,这里就不再赘述,附上我的代码. #include <iostr ...

  3. codeforces Good bye 2016 E 线段树维护dp区间合并

    codeforces Good bye 2016 E 线段树维护dp区间合并 题目大意:给你一个字符串,范围为‘0’~'9',定义一个ugly的串,即串中的子串不能有2016,但是一定要有2017,问 ...

  4. Codeforces 1140F Extending Set of Points 线段树 + 按秩合并并查集 (看题解)

    Extending Set of Points 我们能发现, 如果把x轴y轴看成点, 那么答案就是在各个连通块里面的x轴的个数乘以y轴的个数之和. 然后就变成了一个并查集的问题, 但是这个题目里面有撤 ...

  5. CodeForces 19D Points (线段树+set)

    D. Points time limit per test 2 seconds memory limit per test 256 megabytes input standard input out ...

  6. Codeforces 1140F Extending Set of Points (线段树分治+并查集)

    这题有以下几个步骤 1.离线处理出每个点的作用范围 2.根据线段树得出作用范围 3.根据分治把每个范围内的点记录和处理 #include<bits/stdc++.h> using name ...

  7. codeforces 22E XOR on Segment 线段树

    题目链接: http://codeforces.com/problemset/problem/242/E E. XOR on Segment time limit per test 4 seconds ...

  8. Codeforces 588E. A Simple Task (线段树+计数排序思想)

    题目链接:http://codeforces.com/contest/558/problem/E 题意:有一串字符串,有两个操作:1操作是将l到r的字符串升序排序,0操作是降序排序. 题解:建立26棵 ...

  9. Codeforces Gym 100803G Flipping Parentheses 线段树+二分

    Flipping Parentheses 题目连接: http://codeforces.com/gym/100803/attachments Description A string consist ...

随机推荐

  1. python安装提示No module named setuptools,wget提示ERROR 403: SSL is required

    在下载安装一个python工具时提示报错No module named setuptools [root@kermit supervisor-3.3.0]$ sudo python setup.py ...

  2. 阿里云主机Nginx下配置NodeJS、Express和Forever

    https://cnodejs.org/topic/5059ce39fd37ea6b2f07e1a3 AngularJS中文社区即运行在阿里云主机上,本站使用Nginx引擎,为了AngularJS,我 ...

  3. ASP.NET提供三种主要形式的缓存

    ASP.NET提供三种主要形式的缓存:页面级输出缓存.用户控件级输出缓存(或称为片段缓存)和缓存API.

  4. LA 4108 (线段树)

    区间更新 + 统计更新长度 稍稍不注意就T了 #include<bits/stdc++.h> #define lson l, m, rt<<1 #define rson m+1 ...

  5. HDU 5297

    用x ^ (1 / n) 来求个数,容斥原理 , Num会向后移动, 迭代到不再变化,退出循环 #include<iostream> #include<cstdio> #inc ...

  6. Ex 2_16 给定一个无穷数组..._第二次作业

    先比较数组的A[0]元素,若不相等接下来比较A[1],A[2],A[4],A[8]…,若找到一个区间A[2n-1]<x<A[2n],再对这个区间进行折半查找操作.总的时间为O(logn). ...

  7. Python-Pool类

    目录: multiprocessing模块 Pool类 apply apply_async map close terminate join 进程实例 multiprocessing模块 如果你打算编 ...

  8. mysql 安装问题一:由于找不到MSVCR120.dll,无法继续执行代码.重新安装程序可能会解决此问题。

    这种错误是由于未安装  vcredist  引起的 下载  vcredist  地址:https://www.microsoft.com/zh-CN/download/details.aspx?id= ...

  9. day03 变量 运算符 基本数据类型 输出功能 格式化输出

    变量补充 变量的命名 1变量名的命名的大前提:应该能够反映出变量值所记录的状态 具体的1.变量名由字母数字下划线组成 2.不能以数字开头 3.不能使用关键字命名为变量名 两种写法 1.驼峰体(由字母组 ...

  10. json的转换操作

    toJSON 把JS对象{ 'x': 2, 'y': 3 }转为JSON对象格式的字符串   不能转化字符串 比如"{ 'x': 2, 'y': 3 }" 可以转格式不标准的jso ...