LOJ.6281.数列分块入门5(分块 区间开方)

题目链接
int内的数(也不非得是int)最多开方4.5次就变成1了,所以还不是1就暴力,是1就直接跳过。

#include <cmath>
#include <cstdio>
#include <cctype>
#include <algorithm>
#define gc() getchar()
typedef long long LL;
const int N=5e4+5;

int n,size,bel[N],A[N],tm[N];
LL sum[N];

inline int read()
{
    int now=0,f=1;register char c=gc();
    for(;!isdigit(c);c=gc()) if(c=='-') f=-1;
    for(;isdigit(c);now=now*10+c-'0',c=gc());
    return now*f;
}
LL Query_Sum(int l,int r)
{
    LL res=0ll; int t=std::min(bel[l]*size,r);
    for(int i=l; i<=t; ++i) res+=A[i];
    if(bel[l]!=bel[r])
        for(int i=(bel[r]-1)*size+1; i<=r; ++i)
            res+=A[i];
    for(int i=bel[l]+1; i<bel[r]; ++i) res+=sum[i];
    return res;
}
inline void Sqrt(int p)
{
    sum[bel[p]]-=A[p], A[p]=sqrt(A[p]), sum[bel[p]]+=A[p];
}
void Modify(int l,int r)
{
    int t=std::min(bel[l]*size,r);
    for(int i=l; i<=t; ++i)
        if(A[i]!=1) Sqrt(i);
    if(bel[l]!=bel[r])
        for(int i=(bel[r]-1)*size+1; i<=r; ++i)
            if(A[i]!=1) Sqrt(i);
    for(int i=bel[l]+1; i<bel[r]; ++i)
        if(sum[i]!=size)
            for(int j=(i-1)*size+1; j<=i*size; ++j)
                if(A[j]!=1) Sqrt(j);
}

int main()
{
    n=read(), size=sqrt(n);
    for(int i=1; i<=n; ++i) bel[i]=(i-1)/size+1;
    for(int i=1; i<=n; ++i) A[i]=read(),sum[bel[i]]+=A[i];
    for(int opt,l,r,c,i=1; i<=n; ++i)
    {
        opt=read(),l=read(),r=read(),c=read();
        if(opt) printf("%lld\n",Query_Sum(l,r));
        else Modify(l,r);
    }
    return 0;
}

原文地址:https://www.cnblogs.com/SovietPower/p/8454529.html

时间: 2024-11-06 03:48:50

LOJ.6281.数列分块入门5(分块 区间开方)的相关文章

LOJ#6281. 数列分块入门 5

内存限制:256 MiB时间限制:500 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: hzwer 提交提交记录统计讨论 1 测试数据 题目描述 给出一个长为 nnn 的数列,以及 nnn 个操作,操作涉及区间开方,区间求和. 输入格式 第一行输入一个数字 nnn. 第二行输入 nnn 个数字,第 i 个数字为 aia_ia?i??,以空格隔开. 接下来输入 nnn 行询问,每行输入四个数字 opt\mathrm{opt}opt.lll.rrr.ccc,以空格隔开. 若 opt=

#6283. 数列分块入门 7(区间乘法,区间加法,单点询问)

题目链接:https://loj.ac/problem/6283 题目大意:中文题目 具体思路:和线段树的思路相同,注意lazy的下放,对于不完整的区间,我们需要先更新数组a的值,然后再对数组a进行操作.对于完整的操作,我们要注意优先级,如果原来是a*b+c的话,我们要对这个区间乘以e的话,就表示成a*(b*e)+c*e. AC代码: 1 #include<bits/stdc++.h> 2 using namespace std; 3 # define ll long long 4 const

LiberOJ 6278 数列分块入门 2(分块)

 题解:非常高妙的分块,每个块对应一个桶,桶内元素全部sort过,加值时,对于零散块O(sqrt(n))暴力修改,然后暴力重构桶.对于大块直接整块加.查询时对于非完整块O(sqrt(n))暴力遍历.对于完整的大块用lower_bound或者手写二分log(sqrt(n)查找,总复杂度O(n*sqrt(n)*log(sqrt(n))) 代码如下: #include<cmath> #include<vector> #include<cstdio> #include<c

loj 6278 6279 数列分块入门 2 3

参考:「分块」数列分块入门1 – 9 by hzwer 2 Description 给出一个长为\(n\)的数列,以及\(n\)个操作,操作涉及区间加法,询问区间内小于某个值\(x\)的元素个数. 思路 每个块内保持升序排列. 则块外暴力统计,块内二分查找分界点. 一些注意点,如: 要记录下标: 块外暴力修改完之后需要再排序: 在块内二分查找的值是\(c-tag[i]\)而非\(c\). Code #include <bits/stdc++.h> #define maxn 50010 #def

loj 6277 6280 数列分块入门 1 4

参考:「分块」数列分块入门1 – 9 by hzwer 1 Description 给出一个长为\(n\)的数列,以及\(n\)个操作,操作涉及区间加法,单点查值. 思路 用\(tag\)记录每个块整体的增量. Code #include <bits/stdc++.h> #define maxn 50010 #define F(i, a, b) for (int i = (a); i < (b); ++i) #define F2(i, a, b) for (int i = (a); i

LOJ#6284. 数列分块入门 8

#6284. 数列分块入门 8 内存限制:256 MiB时间限制:500 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: hzwer 提交提交记录统计讨论 1 测试数据 题目描述 给出一个长为 nnn 的数列,以及 nnn 个操作,操作涉及区间询问等于一个数 ccc 的元素,并将这个区间的所有元素改为 ccc. 输入格式 第一行输入一个数字 nnn. 第二行输入 nnn 个数字,第 i 个数字为 aia_ia?i??,以空格隔开. 接下来输入 nnn 行询问,每行输入三个数字 ll

数列分块入门1-9 LibreOJ

数列分块入门1-9 LibreOJ 我也不知道为什么一个大二的ACM选手没学分块. 我怎么记得大一的时候,学长教给我的分块就只有 block 和 num 两个变量来着...好吧,应该是我没认真学.正好前两天朋友给学弟开课,乘机去蹭了一节课.然后...我还是不会哇,菜的一逼塌糊涂. 还是卿学姐好哇,多听几遍,睡得贼香. 分块原理 分块嘛,其实就是优雅的暴力,和莫队(不会)有点异曲同工的赶脚.通过将数组分成小块以降低复杂度. 通常情况下: 每个块的大小(block)为 \(\sqrt{n}\) 块数

LibreOJ6279. 数列分块入门 3 题解

题目链接:https://loj.ac/problem/6279 题目描述 给出一个长为 \(n\) 的数列,以及 \(n\) 个操作,操作涉及区间加法,询问区间内小于某个值 \(x\) 的前驱(比其小的最大元素). 输入格式 第一行输入一个数字 \(n\). 第二行输入 \(n\) 个数字,第 \(i\) 个数字为 \(a_i\),以空格隔开. 接下来输入 \(n\) 行询问,每行输入四个数字 \(opt\).\(l\).\(r\).\(c\),以空格隔开. 若 \(opt=0\),表示将位于

LOJ6277~6285 数列分块入门

Portals 分块需注意的问题 数组大小应为\(N+\sqrt N\),因为最后一个块可能会超出\(N\)的范围.改成记录\(blk,fr,to\)就不用担心这个了 当操作的区间在一个块内时,要特判成暴力修改. 要清楚什么时候应该+tag[t] 最后一个块是越界的,注意是否有影响 数列分块入门 1 给出一个长为\(n\)的数列,以及\(n\)个操作,操作涉及区间加法,单点查值. //数列分块入门 1 #include <cstdio> #include <cmath> inlin