Codeforces 295E 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 线段树的更多相关文章
- Codeforces Beta Round #19D(Points)线段树
D. Points time limit per test 2 seconds memory limit per test 256 megabytes input standard input out ...
- codeforces 295E Yaroslav and Points (离线操作+离散化+区间合并)
参考链接:http://blog.csdn.net/dyx404514/article/details/8817717 写的很详细,这里就不再赘述,附上我的代码. #include <iostr ...
- codeforces Good bye 2016 E 线段树维护dp区间合并
codeforces Good bye 2016 E 线段树维护dp区间合并 题目大意:给你一个字符串,范围为‘0’~'9',定义一个ugly的串,即串中的子串不能有2016,但是一定要有2017,问 ...
- Codeforces 1140F Extending Set of Points 线段树 + 按秩合并并查集 (看题解)
Extending Set of Points 我们能发现, 如果把x轴y轴看成点, 那么答案就是在各个连通块里面的x轴的个数乘以y轴的个数之和. 然后就变成了一个并查集的问题, 但是这个题目里面有撤 ...
- CodeForces 19D Points (线段树+set)
D. Points time limit per test 2 seconds memory limit per test 256 megabytes input standard input out ...
- Codeforces 1140F Extending Set of Points (线段树分治+并查集)
这题有以下几个步骤 1.离线处理出每个点的作用范围 2.根据线段树得出作用范围 3.根据分治把每个范围内的点记录和处理 #include<bits/stdc++.h> using name ...
- codeforces 22E XOR on Segment 线段树
题目链接: http://codeforces.com/problemset/problem/242/E E. XOR on Segment time limit per test 4 seconds ...
- Codeforces 588E. A Simple Task (线段树+计数排序思想)
题目链接:http://codeforces.com/contest/558/problem/E 题意:有一串字符串,有两个操作:1操作是将l到r的字符串升序排序,0操作是降序排序. 题解:建立26棵 ...
- Codeforces Gym 100803G Flipping Parentheses 线段树+二分
Flipping Parentheses 题目连接: http://codeforces.com/gym/100803/attachments Description A string consist ...
随机推荐
- Android来电拦截及来电转移
1. 电话拦截这个功能大家可能都知道了,就是利用反射原理调用ITelephony的隐藏方法来实现.这个就不说了,在附件的代码里有.2.拦截后提示忙音/空号/已关机/已停机这个功能其实是要用到MMI指令 ...
- linux下export命令添加、删除环境变量(转载)
转自:http://blog.csdn.net/shenshendeai/article/details/49794699 export命令 功能说明:设置或显示环境变量. 语 法:export [- ...
- 【原创】大数据基础之Logstash(1)简介、安装、使用
Logstash 6.6.2 官方:https://www.elastic.co/products/logstash 一 简介 Centralize, Transform & Stash Yo ...
- Python基础-入门之路PYTHON-包 相对导入&绝对导入
什么是包 包也是一种模块,但本质上就是一个文件夹 对于使用者而言 使用方式和模块没有任何区别 本质上就是一个文件夹 不同之处在于 多了一个__init__.py 叫包的初始化文件 import导入模块 ...
- mysql-约束介绍
一.约束介绍 约束是一种限制,它通过对表的行或列的数据做出限制,来确保数据的完整性.一致性.约束条件与数据类型宽度一样都是可选参数. 常用约束: PRIMARY KEY (PK) 标识该字段为该表的主 ...
- 通过printf从目标板到调试器的输出
最近在SEGGER的博客上看到Johannes Lask写的一篇关于在调试时使用printf函数从目标MCU输出信息到调试器的文章,自我感觉很有启发,特此翻译此文并推荐给各位同仁.当然限于个人水平,有 ...
- STM32应用实例十一:基于SPI和AD7192的数据采集
在开发臭氧发生器的时,我们需要一个高分辨率的AD采集,于是选择了AD7192,选择这款ADC的原因比较简单.首先它是24位的符合我们的精度要求:其次它自带时钟,便于节省空间:第三他又4路单端或2路差分 ...
- django 中自定义过滤器
多参数过滤器
- Fatal error: Uncaught Error: Call to undefined function mysqli_connect()
sudo apt-get install php7.2-mysql //版本号肯能不一样 sudo dpkg-reconfigure php7.2-mysql sudo /etc/init.d/m ...
- Android Studio 打开activity_main.xml不能正常显示
操作系统:Windows 10 x64 IDE:Android Studio 3.2.1 解决方法:http://www.jcodecraeer.com/a/anzhuokaifa/Android_S ...