bzoj3878: [Ahoi2014&Jsoi2014]奇怪的计算器 线段树

bzoj3878: [Ahoi2014&Jsoi2014]奇怪的计算器

Description

【故事背景】
JYY有个奇怪的计算器,有一天这个计算器坏了,JYY希望你能帮助他写
一个程序来模拟这个计算器的运算。
【问题描述】
JYY的计算器可以执行N条预设好的指令。每次JYY向计算器输入一个正
整数X,计算器就会以X作为初始值,接着依次执行预设的N条指令,最后把
最终得出的结果返回给JYY。
每一条指令可以是以下四种指令之一:(这里a表示一个正整数。)
1、+a:表示将当前的结果加上a;
2、-a:表示将当前的结果减去a;
3、*a:表示将当前的结果乘以a;
4、@a:表示将当前的结果加上a*X(X是一开始JYY输入的数)。
计算器用于记录运算结果的变量的存储范围是有限的,所以每次运算结束之
后会有计算结果溢出的问题。
JYY的计算器中,存储每计算结果的变量只能存储L到R之间的正整数,
如果一次指令执行过后,计算结果超过了R,那么计算器就会自动把结果变成R,然后再以R作为当前结果继续进行之后的计算。同理,如果运算结果小于L,计算器也会把结果变成L,再接着计算。
比如,假设计算器可以存储1到6之间的值,如果当前的计算结果是2,那
么在执行+5操作之后,存储结果的变量中的值将会是6。虽然2+5的实际结
果是7,但是由于7超过了存储范围的上界,所以结果就被自动更正成了上界的大小,也就是6。
JYY一共想在计算器上输入Q个值,他想知道这Q个值输入计算器之后,
分别会得到什么结果呢?

Input

输入文件的第一行包含三个正整数,N,L和R;
第接下来N行,每行一个指令,每个指令如题述,由一个字符和一个正整
数组成,字符和正整数中间有一个空格隔开;
第N+2行包含一个整数Q,表示JYY希望输入的数的数量;
第接下来Q行每行一个正整数,第k个正整数Xk表示JYY在第k次输入的
整数。

Output

输出Q行每行一个正整数,第k行的整数表示输入Xk后,依次经过N个指
令进行计算所得到的结果。

Sample Input

5 1 6
+ 5
- 3
* 2
- 7
@ 2
3
2
1
5

Sample Output

5
3
6

HINT

【样例说明】
当JYY输入2时,计算器会进行5次运算,每一次运算之后得到的结果分
别是6(实际计算结果为7但是超过了上界),3,6,1(实际结果为-1但是低于了下界)和5(由于一开始输入的是2,所以这一次计算为1+2×2)。
1<=N,Q<=10^5,1<=L<=Xk<=R<=10^9,1<=a<=10^9

分析

询问?操作?当然是高级数据结构
乘法?加法?结果?当然是线段树
可是上下界是什么操作?
不会。
怎么办?
我还是会看题解的
当然是要挖掘题目本身的性质。
为什么加和减要分开?
突然发现a是正整数。
为什么a是正整数。
于是就顺利成章地想到了其实所有操作都是单调的(不考虑上下界的情况)(鬼能想得到哦
于是我们考虑吧所有数排序,那么每次越界的都一定是线段树左右两边,可以用类似平衡树的操作维护最大最小值来判断是否超界。
然而这题还有一个难点,就是大标记。
突然发现,我们需要支持的修改操作好多。
乘法,加法,带权乘法,还有区间值的同化。
首先发现带权乘法并不难,因为我们只需要维护区间最小最大值,而每个区间排序后都是单调的,所以最大最小值一定在左右两个端点。只需要用左右端点的值作为权就可以搞定了。
但是我们发现四种标记是稳不行的,因为区间值的同化和乘法和加法的优先级很迷。
这个时候有一个雕虫小技,区间值的同化相当于想把整个区间乘上零然后再加上同化的值,这样的话不仅可以省一个标记,还可以解决优先级的问题。
然后就是码农的日常了。

代码

/**************************************************************
    Problem: 3878
    User: 2014lvzelong
    Language: C++
    Result: Accepted
    Time:1972 ms
    Memory:20632 kb
****************************************************************/
#include<iostream>
#include<cstdlib>
#include<algorithm>
#include<cstdio>
#include<cstring>
using namespace std;
const int N = 110000, M = 110000, T = 440000;
int s[1100];
int read() {
    char ch = getchar(); int x = 0, f = 1;
    while(ch < '0' || ch > '9') {if(ch == '-') f = -1; ch = getchar();}
    while(ch >= '0' && ch <= '9') {x = x * 10 + ch - '0'; ch = getchar();}
    return x * f;
}
int readchar() {char ch = getchar(); while(!s[ch]) ch = getchar(); return s[ch];}
struct data {
    int v, p;
    data(int val = 0, int pos = 0) : v(val), p(pos) {}
    bool operator < (data a) const {return v < a.v;}
}a[N];
int opt[M], x[M], ans[N], n, m, low, hig;
long long mx[T], mn[T], C[T], K[T], B[T];
void Update(int p) {mn[p] = mn[p << 1]; mx[p] = mx[p << 1 | 1];}
void Build(int p, int L, int R) {
    C[p] = 1;
    if(L == R) {mx[p] = mn[p] = a[L].v; return;}
    int mid = L + R >> 1;
    Build(p << 1, L, mid); Build(p << 1 | 1, mid + 1, R);
    Update(p);
}
void Paint(int p, int L, int R, long long c, long long k, long long b) {
    mn[p] = mn[p] * c + a[L].v * k + b;
    mx[p] = mx[p] * c + a[R].v * k + b;
    C[p] *= c; K[p] *= c; B[p] *= c;
    K[p] += k; B[p] += b;
}
void Pushdown(int p, int L, int mid, int R) {
    Paint(p << 1, L, mid, C[p], K[p], B[p]);
    Paint(p << 1 | 1, mid + 1, R, C[p], K[p], B[p]);
    C[p] = 1; K[p] = B[p] = 0;
}
void ChangeL(int p, int L, int R) {
    if(L == R) {Paint(p, L, R, 0, 0, low); return;}
    int mid = L + R >> 1, lc = p << 1, rc = lc | 1;
    Pushdown(p, L, mid, R);
    if(mn[rc] < low) Paint(lc, L, mid, 0, 0, low), ChangeL(rc, mid + 1, R);
    else ChangeL(lc, L, mid);
    Update(p);
}
void ChangeR(int p, int L, int R) {
    if(L == R) {Paint(p, L, R, 0, 0, hig); return;}
    int mid = L + R >> 1, lc = p << 1, rc = lc | 1;
    Pushdown(p, L, mid, R);
    if(mx[lc] > hig) Paint(rc, mid + 1, R, 0, 0, hig), ChangeR(lc, L, mid);
    else ChangeR(rc, mid + 1, R);
    Update(p);
}

void Dfs(int p, int L, int R) {
    if(L == R) {ans[a[L].p] = mn[p]; return ;}
    int mid = L + R >> 1; Pushdown(p, L, mid, R);
    Dfs(p << 1, L, mid); Dfs(p << 1 | 1, mid + 1, R);
}

int main() {
    s['+'] = 1; s['-'] = 2; s['*'] = 3; s['@'] = 4;
    m = read(); low = read(); hig = read();
    for(int i = 1;i <= m; ++i) opt[i] = readchar(), x[i] = read();
    n = read();
    for(int i = 1;i <= n; ++i) a[i] = data(read(), i);
    sort(a + 1, a + n + 1);
    Build(1, 1, n);
    for(int i = 1;i <= m; ++i) {
        if(opt[i] == 1) Paint(1, 1, n, 1, 0, x[i]);
        if(opt[i] == 2) Paint(1, 1, n, 1, 0, -x[i]);
        if(opt[i] == 3) Paint(1, 1, n, x[i], 0, 0);
        if(opt[i] == 4) Paint(1, 1, n, 1, x[i], 0);
        if(mn[1] < low) ChangeL(1, 1, n); 
        if(mx[1] > hig) ChangeR(1, 1, n); 
    }
    Dfs(1, 1, n);
    for(int i = 1;i <= n; ++i) printf("%d\n", ans[i]);
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值