Luogu3793 由乃救爷爷 分块、ST表

传送门



因为昨天写暴力写挂在UOJ上用快排惨遭卡常,所以今天准备写一个卡常题消遣消遣,然后时间又垫底了QAQ

这道题显然需要支持一个\(O(N)\)预处理\(O(1)\)查询的ST表,显然普通的ST表是做不到的,因为预处理的时间太长了

于是分块优化掉ST表的预处理

约定\(L_i,R_i\)表示第\(i\)个块的左端点和右端点,\(be_i\)表示第\(i\)个数所在的块

对于每一个位置\(i\)预处理\(L_{be_i}\)到\(i\)的所有数的最大值\(lmax_i\)以及\(i\)到\(R_{be_i}\)的所有数的最大值\(rmax_i\);预处理块最大值的ST表

对于一个询问\((l,r)\),如果\(be_i \neq be_r\)查询\(rmax_l , lmax_r\)以及\(be_l+1\)到\(be_r-1\)的所有块的最大值,三者取\(max\);如果\(be_l = be_r\)就暴力

考虑复杂度,设块长为\(S\),预处理复杂度瓶颈为ST表复杂度\(O(\frac{N}{S}log\frac{N}{S})\);询问复杂度上,如果\(be_l \neq be_r\)为\(O(1)\),否则为\(O(S)\)。因为数据随机,出现\(be_l = be_r\)的概率为\(\frac{S}{N}\),所以最坏查询复杂度为\(O(S^2)\),总复杂度为\(O(\frac{N}{S}log\frac{N}{S} + S^2)\)

差不多当\(S = \sqrt[3]{NlogN}\)的时候有最优复杂度,但是因为数据随机所以块长开大一点也没关系

然后可能需要一些奇怪的常数优化比如把max define掉之类的

#include<bits/stdc++.h>
using namespace std;

namespace GenHelper
{
    unsigned z1,z2,z3,z4,b;
    unsigned rand_()
    {
    b=((z1<<6)^z1)>>13;
    z1=((z1&4294967294U)<<18)^b;
    b=((z2<<2)^z2)>>27;
    z2=((z2&4294967288U)<<2)^b;
    b=((z3<<13)^z3)>>21;
    z3=((z3&4294967280U)<<7)^b;
    b=((z4<<3)^z4)>>12;
    z4=((z4&4294967168U)<<13)^b;
    return (z1^z2^z3^z4);
    }
}
void srand(unsigned x)
{
    using namespace GenHelper;
    z1=x; z2=(~x)^0x233333333U; z3=x^0x1234598766U; z4=(~x)+51;
}
int read()
{
    using namespace GenHelper;
    int a=rand_()&32767;
    int b=rand_()&32767;
    return a*32768+b;
}

const int MAXN = 2e7 + 7;
unsigned int arr[MAXN] , maxN1[MAXN] , maxN2[MAXN] , ST[15][26000] , logg2[26000];
unsigned int N , M , S , T , cnt;

#define cmp(a , b) ((a) > (b) ? (a) : (b))

inline unsigned int qST(int x , int y){
    if(x > y)
        return 0;
    int t = logg2[y - x + 1];
    return cmp(ST[t][x] , ST[t][y - (1 << t) + 1]);
}

int main(){
    cin >> N >> M >> S;
    srand(S);
    T = pow(N * log2(N) , 1.0/3);
    cnt = N / T + (N % T ? 1 : 0);
    for(int i = 1 ; i <= N ; ++i){
        arr[i] = read();
        maxN2[i] = arr[i];
        if(i % T != 1)
            maxN2[i] = cmp(maxN2[i] , maxN2[i - 1]);
    }
    for(int i = N ; i ; --i){
        maxN1[i] = arr[i];
        if(i % T)
            maxN1[i] = cmp(maxN1[i] , maxN1[i + 1]);
        if(i % T == 1)
            ST[0][i / T + 1] = maxN1[i];
    }
    for(int i = 2 ; i <= cnt ; ++i)
        logg2[i] = logg2[i >> 1] + 1;
    for(int i = 1 ; 1 << i <= cnt ; ++i)
        for(int j = 1 ; j + (1 << i) - 1 <= cnt ; ++j)
            ST[i][j] = cmp(ST[i - 1][j] , ST[i - 1][j + (1 << (i - 1))]);
    unsigned long long sum = 0;
    while(M--){
        int l = read() % N + 1 , r = read() % N + 1;
        if(l > r)
            swap(l , r);
        if((l - 1) / T == (r - 1) / T){
            unsigned int ans = 0;
            for(int i = l ; i <= r ; ++i)
                ans = cmp(ans , arr[i]);
            sum += ans;
        }
        else
            sum += cmp(cmp(maxN1[l] , maxN2[r]) , qST((l - 1) / T + 2 , (r - 1) / T));
    }
    cout << sum;
    return 0;
}

原文地址:https://www.cnblogs.com/Itst/p/10354018.html

时间: 2024-11-08 10:49:04

Luogu3793 由乃救爷爷 分块、ST表的相关文章

Luogu 3793 由乃救爷爷

\(\verb|Luogu 3793 由乃救爷爷|\) rmq,数据随机 \(n,\ m\leq 2\times10^7\) lxl ST表 分块,大小设为 \(x\) 预处理每个块两端到块内每个点的前缀 \(\max\) 和后缀 \(\max\) 预处理块间ST表 数据随机 就成了期望 \(O(n)\) 的rmq 重点是没人卡你,卡还不一定卡的住,还要冒着被暴力AC的风险 然后就愉悦地~ #include <bits/stdc++.h> using namespace std; #defin

[BZOJ1012] [JSOI2008] 最大数maxnumber (ST表)

Description 现在请求你维护一个数列,要求提供以下两种操作:1. 查询操作.语法:Q L 功能:查询当前数列中末尾L个数中的最大的数,并输出这个数的值.限制:L不超过当前数列的长度.2. 插入操作.语法:A n 功能:将n加上t,其中t是最近一次查询操作的答案(如果还未执行过查询操作,则t=0),并将所得结果对一个固定的常数D取模,将所得答案插入到数列的末尾.限制:n是非负整数并且在长整范围内.注意:初始时数列是空的,没有一个数. Input 第一行两个整数,M和D,其中M表示操作的个

RMQ问题 - ST表的简单应用

2017-08-26 22:25:57 writer:pprp 题意很简单,给你一串数字,问你给定区间中最大值减去给定区间中的最小值是多少? 用ST表即可实现 一开始无脑套模板,找了最大值,找了最小值,分别用两个函数实现,实际上十分冗余 所以TLE了 之后改成一个函数中同时处理最大值和最小值,就可以了 AC代码如下: /* @theme:poj 3264 @writer:pprp @declare:ST表(sparse table)稀疏表,用动态规划的思想来解决RMQ问题: @date:2017

[模板]ST表浅析

ST表,稀疏表,用于求解经典的RMQ问题.即区间最值问题. Problem: 给定n个数和q个询问,对于给定的每个询问有l,r,求区间[l,r]的最大值.. Solution: 主要思想是倍增和区间dp. 状态:dp[i][j] 为闭区间[i,i+2^j-1]的最值. 这个状态与转移方程的关系很大,即闭区间的范围涉及到了转移方程的简便性. 转移方程:dp[i][j]=max(dp[i][j-1],dp[i+2^(j-1)][j-1]). 这是显然的,但这里有个细节:第一个项的范围为[i,i+2^

【模板】ST表 洛谷P1816 忠诚

P1816 忠诚 题目描述 老管家是一个聪明能干的人.他为财主工作了整整10年,财主为了让自已账目更加清楚.要求管家每天记k次账,由于 管家聪明能干,因而管家总是让财主十分满意.但是由于一些人的挑拨,财主还是对管家产生了怀疑.于是他决定用一种特别的方法来判断管家的忠诚,他把每次的 账目按1,2,3…编号,然后不定时的问管家问题,问题是这样的:在a到b号账中最少的一笔是多少?为了让管家没时间作假他总是一次问多个问题. 输入输出格式 输入格式: 输入中第一行有两个数m,n表示有m(m<=100000

GCD(st表)

GCD Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)Total Submission(s): 3432    Accepted Submission(s): 1227 Problem Description Give you a sequence of N(N≤100,000) integers : a1,...,an(0<ai≤1000,000,000). There a

UVA 11475 Extend to Palindrome(后缀数组+ST表)

[题目链接] http://acm.hust.edu.cn/vjudge/problem/27647 [题目大意] 给出一个字符串,要求在其后面添加最少的字符数,使得其成为一个回文串.并输出这个回文串. [题解] 用拼接符将原字符串倒序相接,做一遍后缀数组,查询两串相应位置的LCP就是以该点为中心的回文串长度的一半分,奇偶求出所有后缀回文串,保留最长的,则补充部分为剩下的前缀的倒置.至于查询两串的LCP我们可以在height数组建立ST表查询. [代码] #include <cstdio> #

lca最近公共祖先(st表)

大体思路 1.求出每个元素在树中的深度 2.用st表预处理的方法处理出f[i][j],f[i][j]表示元素i上方第2^j行对应的祖先是谁 3.将较深的点向上挪,直到两结点的深度相同 4.深度相同后,祖先可能就在上方,再走几步就到了,于是两个点同时向上移 具体的方法和代码贴在下面 ↓ 具体来看 1.求出每个元素在树中的深度 //求每个节点在树中的深度 void dfs(int pos,int pre)//pre是pos的父节点 { for(int i=0;i<v[pos].size;i++)//

【BZOJ-3956】Count ST表 + 单调栈

3956: Count Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 173  Solved: 99[Submit][Status][Discuss] Description Input Output Sample Input 3 2 0 2 1 2 1 1 1 3 Sample Output 0 3 HINT M,N<=3*10^5,Ai<=10^9 Source CH Round#64 MFOI杯水题欢乐赛day1 By Gromah So