CF803G - Periodic RMQ Problem 动态开点线段树或离线

题意

有一个长度为n × k （<=1E9）的数组，有区间修改和区间查询最小值的操作。

思路

由于数组过大，直接做显然不行。

有两种做法，可以用动态开点版本的线段树，或者离线搞（还没搞）（搞好了）。

注意只有1E5次操作，所以真正被更新到的区间并不多，最差单次新开2×log（1E9）。

对于新开的区间的最小值，可以这样计算，如果区间表示的值大于n，那就是原来长度为n的区间的最小值，小于n的话，在ST表中查询即可。

#include <bits/stdc++.h>

using namespace std;

#define pb push_back

#define fi first

#define se second

#define debug(x) cerr<<#x << " := " << x << endl;

#define bug cerr<<"-----------------------"<<endl;

#define FOR(a, b, c) for(int a = b; a <= c; ++ a)

typedef long long ll;

typedef long double ld;

typedef pair<int, int> pii;

typedef pair<ll, ll> pll;

template<class T> void _R(T &x) { cin >> x; }

void _R(int &x) { scanf("%d", &x); }

void _R(ll &x) { scanf("%lld", &x); }

void _R(double &x) { scanf("%lf", &x); }

void _R(char &x) { scanf(" %c", &x); }

void _R(char *x) { scanf("%s", x); }

void R() {}

template<class T, class... U> void R(T &head, U &... tail) { _R(head); R(tail...); }

template<typename T>

inline T read(T&x){

    x=;int f=;char ch=getchar();

    while (ch<''||ch>'') f|=(ch=='-'),ch=getchar();

    while (ch>=''&&ch<='') x=x*+ch-'',ch=getchar();

    return x=f?-x:x;

}

const int inf = 0x3f3f3f3f;

const int mod = 1e9+;

/**********showtime************/

            const int maxn = 1e5+;

            int b[maxn];

            int st[maxn][], Log[maxn];

            void init_st(int n) {

                Log[] = -;

                for(int i=; i<=n; i++) {

                    Log[i] = Log[i>>] + ;

                    st[i][] = b[i];

                }

                for(int j=; (<<j) <= n; j++) {

                    for(int i=; i + (<<j) - <= n; i++) {

                        st[i][j] = min(st[i][j-], st[i+ (<<(j-))][j-]);

                    }

                }

            }

            int rmq_st(int L, int R) {

                int k = Log[R - L + ];

                return min(st[L][k], st[R-(<<k)+][k]);

            }

            int n,k;

            int getmin(int le, int ri) {

                if(ri - le +  >= n) return rmq_st(, n);

                int L = le % n; if(L == ) L = n;

                int R = ri % n; if(R == ) R = n;

                if(L <= R) return rmq_st(L, R);

                return min(rmq_st(L, n), rmq_st(, R));

            }

            struct Node{

                int le, ri;

                int lc,rc;

                int val, tag;

            } tree[maxn * ];

            int tot = ;

            int newNode(int le, int ri) {

                tot++;

                tree[tot].le = le;

                tree[tot].ri = ri;

                tree[tot].lc =  tree[tot].rc = ;

                tree[tot].val = getmin(le, ri);

                tree[tot].tag = ;

                return tot;

            }

            void pushdown(int rt) {

                tree[tree[rt].lc].val = tree[tree[rt].lc].tag = tree[rt].tag;

                tree[tree[rt].rc].val = tree[tree[rt].rc].tag = tree[rt].tag;

                tree[rt].tag = ;

            }

            void update(int L, int R, int b, int rt) {

                if(L<=tree[rt].le && tree[rt].ri <= R) {

                    tree[rt].val = tree[rt].tag = b;

                    return;

                }

                int mid = (tree[rt].le + tree[rt].ri) >> ;

                if(tree[rt].lc == ) tree[rt].lc = newNode(tree[rt].le, mid);

                if(tree[rt].rc == ) tree[rt].rc = newNode(mid+, tree[rt].ri);

                if(tree[rt].tag) pushdown(rt);

                if(mid >= L) update(L, R, b, tree[rt].lc);

                if(mid < R) update(L, R, b, tree[rt].rc);

                tree[rt].val = min(tree[tree[rt].lc].val, tree[tree[rt].rc].val);

            }

            int query(int L, int R, int rt) {

                if(L <= tree[rt].le && tree[rt].ri <= R) {

                    return tree[rt].val;

                }

                int mid = (tree[rt].le + tree[rt].ri) >> ;

                if(tree[rt].lc == ) tree[rt].lc = newNode(tree[rt].le, mid);

                if(tree[rt].rc == ) tree[rt].rc = newNode(mid+, tree[rt].ri);

                if(tree[rt].tag) pushdown(rt);

                int res = inf;

                if(mid >= L) res = min(res, query(L, R, tree[rt].lc));

                if(mid < R) res = min(res, query(L, R, tree[rt].rc));

                tree[rt].val = min(tree[tree[rt].lc].val, tree[tree[rt].rc].val);

                return res;

            }

int main(){

            scanf("%d%d", &n, &k);

            for(int i=; i<=n; i++) scanf("%d", &b[i]);

            init_st(n);

            int q;  scanf("%d", &q);

            newNode(, n * k);

            while(q--) {

                int op;

                scanf("%d", &op);

                if(op == ) {

                    int le,ri,x;

                    scanf("%d%d%d", &le, &ri, &x);

                    update(le, ri, x, );

                }

                else {

                    int le, ri;

                    scanf("%d%d", &le, &ri);

                    printf("%d\n", query(le, ri, ));

                }

            }

            return ;

}

离线的话，我们可以记录下所有被问到的点，然后我们可以压缩原来长度为1E9的数组：问到的点保持不变，而相邻两点之间的区间压缩成一个点，保存这段区间的最小值即可。

#include <iostream>

#include <vector>

#include <queue>

#include <algorithm>

using namespace std;

#define pb push_back

#define fi first

#define se second

#define debug(x) cerr<<#x << " := " << x << endl;

#define bug cerr<<"-----------------------"<<endl;

#define FOR(a, b, c) for(int a = b; a <= c; ++ a)

typedef long long ll;

typedef long double ld;

typedef pair<int, int> pii;

typedef pair<ll, ll> pll;

template<class T> void _R(T &x) { cin >> x; }

void _R(int &x) { scanf("%d", &x); }

void _R(ll &x) { scanf("%lld", &x); }

void _R(double &x) { scanf("%lf", &x); }

void _R(char &x) { scanf(" %c", &x); }

void _R(char *x) { scanf("%s", x); }

void R() {}

template<class T, class... U> void R(T &head, U &... tail) { _R(head); R(tail...); }

template<typename T>

inline T read(T&x){

    x=;int f=;char ch=getchar();

    while (ch<''||ch>'') f|=(ch=='-'),ch=getchar();

    while (ch>=''&&ch<='') x=x*+ch-'',ch=getchar();

    return x=f?-x:x;

}

const int inf = 0x3f3f3f3f;

const int mod = 1e9+;

/**********showtime************/

            const int maxn = 1e5+;

            int b[maxn];

            struct Query{

                int op;

                int le,ri;

                int x;

            } ask[maxn];

            vector<int>v, nv,node;

            int n,k;

            int st[maxn][];

            int Log[maxn];

            void init_st(int n){

                Log[] = -;

                for(int i = ; i<=n; i++) {

                    st[i][] = b[i];

                    Log[i] = Log[i>>] + ;

                }

                for(int j=; (<<j) <= n; j++) {

                    for(int i=; i + ( << j) - <= n; i ++) {

                        st[i][j] = min(st[i][j-], st[i + (<<(j-)) ][j-]);

                    }

                }

            }

            int get(int L, int R) {

                int k = Log[R - L + ];

                return min(st[L][k], st[R-(<<k)+][k]);

            }

            int getmin(int le, int ri) {

                if(ri - le + >= n)

                    return get(, n);

                int l = le % n; if(!l) l = n;

                int r = ri % n; if(!r) r = n;

                if(l <= r) return get(l, r);

                return min(get(l, n), get(, r));

            }

            int getid(int val) {

                return lower_bound(nv.begin(), nv.end(), val) - nv.begin() + ;

            }

            int mn[maxn*],lazy[maxn * ];

            void build(int le, int ri, int rt) {

                if(le == ri) {

                    mn[rt] = node[le-];

                    return;

                }

                int mid = (le + ri) >> ;

                build(le, mid, rt<<);

                build(mid+, ri, rt<<|);

                mn[rt] = min(mn[rt<<], mn[rt<<|]);

            }

            void pushdown(int rt) {

                mn[rt<<] = mn[rt<<|] = lazy[rt];

                lazy[rt<<] = lazy[rt<<|] = lazy[rt];

                lazy[rt] = ;

            }

            void update(int L, int R,int x, int le, int ri, int rt) {

                if(le >= L && ri <= R){

                    lazy[rt] = x;

                    mn[rt] = x;

                    return ;

                }

                if(lazy[rt]) pushdown(rt);

                int mid = (le + ri) >> ;

                if(mid >= L) update(L, R, x, le, mid, rt<<);

                if(mid < R) update(L, R, x, mid+, ri, rt<<|);

                mn[rt] = min(mn[rt<<], mn[rt<<|]);

            }

            int query(int L, int R,int le, int ri, int rt) {

                if(le >= L && ri <= R) {

                    return mn[rt];

                }

                if(lazy[rt]) pushdown(rt);

                int mid = (le + ri) >> ;

                int res = inf;

                if(mid >= L) res = min(res, query(L, R, le, mid, rt<<));

                if(mid < R) res = min(res, query(L, R, mid+,ri, rt<<|));

                mn[rt] = min(mn[rt<<], mn[rt<<|]);

                return res;

            }

int main(){

            scanf("%d%d", &n, &k);

            for(int i=; i<=n; i++) scanf("%d", &b[i]);

            init_st(n);

            int q;  scanf("%d", &q);

            for(int i=; i<=q; i++){

                int op;

                scanf("%d", &op);

                if(op == ) {

                    ask[i].op = op;

                    scanf("%d%d%d", &ask[i].le, &ask[i].ri, &ask[i].x);

                    v.pb(ask[i].le);

                    v.pb(ask[i].ri);

                }

                else{

                    ask[i].op = op;

                    scanf("%d%d", &ask[i].le, &ask[i].ri);

                    v.pb(ask[i].le);

                    v.pb(ask[i].ri);

                }

            }

            sort(v.begin(), v.end());

            v.erase(unique(v.begin(), v.end()), v.end());

            int N = v.size();

            // cout<<" ** "<<endl;

            for(int i=; i<N; i++) {

                node.pb( getmin(v[i],v[i]));

                nv.pb(v[i]);

                if(i+ < N && v[i] +  <= v[i+] - )

                {

                    node.pb(getmin(v[i]+, v[i+]-));

                    nv.pb(v[i]+);

                }

            }   

            /// 把一个开区间当成一个点。

            N = nv.size();

            build(, N, );

            for(int i=; i<=q; i++) {

                if(ask[i].op == ) {

                    update(getid(ask[i].le),getid(ask[i].ri), ask[i].x, , N, );

                }

                else {

                    printf("%d\n",query(getid(ask[i].le), getid(ask[i].ri), , N, ));

                }

            }

            return ;

}

CF803G - Periodic RMQ Problem 动态开点线段树或离线的更多相关文章

Codeforces 803G Periodic RMQ Problem ST表+动态开节点线段树
思路: (我也不知道这是不是正解) ST表预处理出来原数列的两点之间的min 再搞一个动态开节点线段树节点记录ans 和标记 lazy=-1 当前节点的ans可用 lazy=0 没被覆盖过 els ...
【bzoj4999】This Problem Is Too Simple！树链剖分+动态开点线段树
题目描述给您一颗树,每个节点有个初始值. 现在支持以下两种操作: 1. C i x(0<=x<2^31) 表示将i节点的值改为x. 2. Q i j x(0<=x<2^31) ...
[bzoj 3531][SDOI2014]旅行（树链剖分+动态开点线段树）
题目:http://www.lydsy.com:808/JudgeOnline/problem.php?id=3531 分析: 对于每个颜色(颜色<=10^5)都建立一颗线段树什么!那么不是M ...
codeforces 893F - Physical Education Lessons 动态开点线段树合并
https://codeforces.com/contest/893/problem/F 题意: 给一个有根树, 多次查询,每次查询对于$x$i点的子树中,距离$x$小于等于$k$的所有点中权值最小的 ...
HDU 6183 Color it(动态开点线段树)
题目原网址:http://acm.hdu.edu.cn/showproblem.php?pid=6183 题目中文翻译: Time Limit: 20000/10000 MS (Java/Others ...
[2016湖南长沙培训Day4][前鬼后鬼的守护 chen] (动态开点线段树+中位数 or 动规 or 贪心+堆优化)
题目大意给定一个长度为n的正整数序列,令修改一个数的代价为修改前后两个数的绝对值之差,求用最小代价将序列转换为不减序列. 其中,n满足小于500000,序列中的正整数小于10^9 题解(引自mzx神 ...
【BZOJ-4636】蒟蒻的数列动态开点线段树 ||（离散化） + 标记永久化
4636: 蒟蒻的数列 Time Limit: 30 Sec Memory Limit: 256 MBSubmit: 247 Solved: 113[Submit][Status][Discuss ...
codeforces 915E - Physical Education Lessons 动态开点线段树
题意: 最大$10^9$的区间, $3*10^5$次区间修改,每次操作后求整个区间的和题解: 裸的动态开点线段树,计算清楚数据范围是关键... 经过尝试 $2*10^7$会$MLE$ $10^7$会 ...
CF915E Physical Education Lessons 动态开点线段树
题目链接 CF915E Physical Education Lessons 题解动态开点线段树代码 /* 动态开点线段树 */ #include<cstdio> #include&l ...

随机推荐

SpringBoot日志相关
SpringBoot使用的是SLF4j当门面,Logback当实现完成日志级别数字越大,级别越高,框架只会输出大于等于当前日志级别的信息 ERROR 40 WARN 30 INFO 20 DEBU ...
HPU暑期集训积分赛1
A. Nth power of n 单点时限: 1.0 sec 内存限制: 512 MB 求 nn 的个位数. 输入格式多组输入,处理到文件结束.每组数据输入一个 n.(1≤n≤109) 输出格式 ...
应用性能测试神器 Gatling，你用过吗？
在应用程序上线之前,有多少人做过性能测试? 估计大部分开发者更多地关注功能测试,并且会提供一些单元测试和集成测试的用例.然而,有时候性能漏洞导致的影响比未发现的业务漏洞更严重,因为性能漏洞影响的是整个 ...
常用服务部署脚本(nodejs,pyenv,go,redis,)
根据工作总结的常用安装脚本,要求linux-64系统 #!/bin/bash path=/usr/local/src node () { cd $path #wget https://nodejs.o ...
iOS开发 8小时时差问题
今天调试遇到时间计算的问题,发现怎么算都会有差别,后来仔细观察,发现有8小时的时差…… 这篇文章解释的很好,用到了,因此记之. ios有关时间打印出来差8小时的问题
HTML第六章盒子模型
什么是盒子模型: (1)边框: (2)内边距: (3)外边距: (4)元素内容·: (5)背景色·: 边框: 属性: 颜色(border-color),粗细(border-width),样式(bord ...
selenium操作cookies实现免密登录，自动发微博
一直想用selenium实现个小功能,比如发微博之类的,但是有的网站在登录会有验证码,没想到太好的方法解决,于是想到利用cookies来登录网站第一步:获取一个可用的cookies,获取的cooki ...
koa2基于stream(流)进行文件上传和下载
阅读目录一:上传文件(包括单个文件或多个文件上传) 二:下载文件回到顶部一:上传文件(包括单个文件或多个文件上传) 在之前一篇文章,我们了解到nodejs中的流的概念,也了解到了使用流的优点,具 ...
n的阶乘 -牛客
题目描述输入一个整数n,输出n的阶乘(每组测试用例可能包含多组数据,请注意处理) 输入描述: 一个整数n(1<=n<=20) 输出描述: n的阶乘解题思路采用递归求解,也可以使用循环 ...
01、VM安装教程
1.运行下载完成的Vmware Workstation虚拟机软件包,将会看到如图所示,然后点击“下一步”按钮, 2.在最终用户许可协议界面选中“我接受许可协议中的条款”复选框,然后点击“下一步”按钮 ...

CF803G - Periodic RMQ Problem 动态开点线段树 或 离线

题意

思路

CF803G - Periodic RMQ Problem 动态开点线段树 或 离线的更多相关文章

随机推荐

热门专题

CF803G - Periodic RMQ Problem 动态开点线段树或离线

CF803G - Periodic RMQ Problem 动态开点线段树或离线的更多相关文章