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 ...
随机推荐
- centOS docker运行Asp.net Core程序
[root@localhost chenhua]# docker run -it --rm -p : --name aspnetcore_sample microsoft/dotnet-samples ...
- C#一元二次方程
- Windows 批处理大全(附各种实例)
Windows 批处理大全(附各种实例) 2009年07月19日 21:31:00 阅读数:2552 批处理文件是无格式的文本文件,它包含一条或多条命令.它的文件扩展名为 .bat 或 .cmd.在命 ...
- Mybatis调用PostgreSQL存储过程实现数组入参传递
注:本文来源于 < Mybatis调用PostgreSQL存储过程实现数组入参传递 > 前言 项目中用到了Mybatis调用PostgreSQL存储过程(自定义函数)相关操作,由于Pos ...
- Confluence 6 创建站点的导出文件
希望为你的站点创建一个 XML 导出文件: 进入 > 基本配置(General Configuration) > 备份和恢复(Backup & Restore). 选择 归档到备 ...
- Confluence 6 自定义管理员联系信息
你可以自定义在 联系站点管理员(Contact Site Administrators)页面中显示的消息. 希望编辑这个管理员联系消息: 在屏幕的右上角单击 控制台按钮 ,然后选择 General C ...
- MySQL、MongoDB、Redis 数据库之间的区别
NoSQL 的全称是 Not Only SQL,也可以理解非关系型的数据库,是一种新型的革命式的数据库设计方式,不过它不是为了取代传统的关系型数据库而被设计的,它们分别代表了不同的数据库设计思路. M ...
- 使用react 在页面上引用静态图片,图片不显示
const url='../assets/logo.png'; <img src={url} alt=''/> 1.使用require <img src={require('../a ...
- 原码、补码、反码的概念和java数的存储方式
原码:用符号位和数值位表示一个带符号数,整数符号->0,负数符号->1,数值一般用二进制形式表示 [+10011]原=00010011 [-10011]原=10010011 反码:正 ...
- loadrunner出现报错operands of = have illegal types `pointer to char' and `int'
原始代码: void split(char * p,char * str){ /* 传入一个数组进行p和一个以什么进行分割的str,返回切片后的值 */ int i = 0, j = 0; char ...