A Simple Problem with Integers----线段树

最新推荐文章于 2023-07-02 16:47:46 发布
最新推荐文章于 2023-07-02 16:47:46 发布
阅读量3.1k 收藏 5
点赞数 3
分类专栏： ACM_数据结构
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/bless924295/article/details/51569413
版权
C - A Simple Problem with Integers
Time Limit:5000MS     Memory Limit:131072KB     64bit IO Format:%I64d & %I64u

Description

You have N integers, A1, A2, ... , AN. You need to deal with two kinds of operations. One type of operation is to add some given number to each number in a given interval. The other is to ask for the sum of numbers in a given interval.

Input

The first line contains two numbers N and Q. 1 ≤ N,Q ≤ 100000.
The second line contains N numbers, the initial values of A1, A2, ... , AN. -1000000000 ≤ Ai ≤ 1000000000.
Each of the next Q lines represents an operation.
"C abc" means adding c to each of Aa, Aa+1, ... , Ab. -10000 ≤ c ≤ 10000.
"Q ab" means querying the sum of Aa, Aa+1, ... , Ab.

Output

You need to answer all Q commands in order. One answer in a line.

Sample Input

10 5
1 2 3 4 5 6 7 8 9 10
Q 4 4
Q 1 10
Q 2 4
C 3 6 3
Q 2 4

Sample Output

4
55
9
15


题目链接：http://acm.hust.edu.cn/vjudge/contest/view.action?cid=66989#problem/C


kuangbin大神的专题，这两天比赛也多，考试也多，好久都没写点东西了，今晚没有什么太大的事，正值集训队休息，敲篇博客，看看手生没生，然后准备准备今晚的codeforces，哎，好晚的说---------||||||||


这个题是裸的区间线段树更新的题，在讲区间线段树的时候我想应该先讲一下单点线段树的更新，比较好理解。

放一个裸的单点线段树更新的题：http://acm.sdut.edu.cn/sdutoj/problem.php?action=showproblem&problemid=3397

这个是我刚开始学代码有一次比赛学长们出的防AK的题，恩，楠神出的，卡了我好久。。。刚开始学你就出线段树，楠神你也是绝了。。。

我是用静态方法构建线段树，即用数组来做（毕竟下标比较好用），左子树下标为2×i，右子树为2×i+1，不要问我为什么，完全二叉树的性质，保险起见，数组一般开4×n，具体为什么仔细一想就明白了，因为是完全二叉树，所以第一个叶子节点的下标为大于等于n的第一个2^i，举例来说n==33，所以第一个叶子节点的下标为64（2^6），所以左孩子为128，右孩子是129，所以是4倍。其实实际用到的是2×n+1个，有一些空间用不到。

有几个函数在这里解释一下是什么意思PushUp（由下往上更新修改后的数据），Build（初始化线段树），Update（修改叶子结点的数据），Query（查询区间数据）


百度图片里找了张图片，很清晰，如果你能够在写代码的时候脑子里一直想着树的形状会很好写

跟着代码走一边过程，不用我讲你就会单点更新了

#include <cstdio>  
#include <iostream>  
#include <cstring>  
#include <algorithm>  
using namespace std;  
int a[6000000];  
void PushUp(int st)  
{  
    a[st]=a[st<<1]+a[st<<1|1];  
}  
void Build(int L,int R,int st)  
{  
    a[st]=0;  
    if(L==R)  
    {  
        scanf("%d",&a[st]);  
        return ;  
    }  
    int mid=(L+R)>>1;  
    Build(L,mid,st<<1);  
    Build(mid+1,R,st<<1|1);  
    PushUp(st);  
}  
void Updata(int L,int R,int st,int x,int d)  
{  
    if(L==R)  
    {  
        a[st]+=d;  
        return ;  
    }  
    int mid=(L+R)>>1;  
    if(x<=mid)  
        Updata(L,mid,st<<1,x,d);  
    else  
        Updata(mid+1,R,st<<1|1,x,d);  
    PushUp(st);  
}  
  
int Query(int L,int R,int st,int l,int r)  
{  
    if(R<l||L>r)  
        return 0;  
    if(L>=l&&R<=r)  
        return a[st];  
    int mid=(L+R)>>1;  
    int ans=0;  
    if(l<=mid)  
        ans+=Query(L,mid,st<<1,l,r);  
    if(r>mid)  
        ans+=Query(mid+1,R,st<<1|1,l,r);  
    return ans;  
}  
int main()  
{  
    int n,m;  
    char str[100];  
    int x,d,l,r;  
    while(~scanf("%d%d",&n,&m))  
    {  
        Build(1,n,1);  
        while(m--)  
        {  
            scanf("%s",str);  
            if(!strcmp(str,"UPDATE"))  
            {  
                scanf("%d%d",&x,&d);  
                Updata(1,n,1,x,d);  
            }  
            else  
            {  
                scanf("%d%d",&l,&r);  
                printf("%d\n",Query(1,n,1,l,r));  
            }  
        }  
        cout<<endl;  
    }  
    return 0;  
}

好吧，我承认上面的代码有点啰嗦，不过好入门，入了门以后来看我简化版的

其实原理是一样的，看这个题http://acm.hust.edu.cn/vjudge/contest/view.action?cid=66989#problem/A

经典线段树单点更新

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <cmath>
#include <stack>
#include <deque>
#define eps 0.00000001
using namespace std;
int a[600000];
int c,b;
void PushUp(int st)
{
    while(st--)
        a[st]=a[st<<1]+a[st<<1|1];
}
void Updata(int x,int y)
{
    while(x)
    {
        a[x]+=y;
        if(a[x]<0)
            a[x]=0;
        x>>=1;
    }
}
int Query(int k,int x,int y)
{
    int mid=x+((y-x)>>1),ans=0;
    if(c<=x&&b>=y) return a[k];
    if(c<=mid) ans+=Query(k<<1,x,mid);
    if(mid<b) ans+=Query(k*2+1,mid+1,y);
    return ans;
}
int main()
{
    int t;
    int x,y;
    scanf("%d",&t);
    int n;
    int top=0;
    while(t--)
    {
        memset(a,0,sizeof(a));
        top++;
        scanf("%d",&n);
        int h=1;
        while(1)
        {
            if(h>=n)
                break;
            h<<=1;
        }
        int i;
        for(i=0;i<n;i++)
        {
            scanf("%d",&a[h+i]);
        }
        printf("Case %d:\n",top);
        PushUp(h);
        char s[1000];
        while(1)
        {
            scanf("%s",s);
            if(!strcmp(s,"End"))
                break;
            else if(!strcmp(s,"Query"))
            {
                scanf("%d%d",&c,&b);
                cout<<Query(1,1,h)<<endl;
            }
            else if(!strcmp(s,"Add"))
            {
                scanf("%d%d",&x,&y);
                Updata(h+x-1,y);
            }
            else if(!strcmp(s,"Sub"))
            {
                scanf("%d%d",&x,&y);
                y=-y;
                Updata(h+x-1,y);
            }
        }
    }
    return 0;
}
这个就简洁好多了吧，懂了下标后好好利用就行。

这里我放一个题，就当练练手吧http://acm.hust.edu.cn/vjudge/contest/view.action?cid=66989#problem/B

好了，讲了这么多，单点线段树算是讲完了，我们开看区间更新线段树。

原理是差不多的，有点不太一样，每个节点增加一个add变量来记录增量的累加，因为每一个数都更新到叶子节点是很不明智的选择，所以，我们找到那个区间，记录下来。

我记得网上有一张神图清楚明了，但是我找了好久没找到，当时我也是看到了那张神图恍然大悟。。。

上代码吧

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#define eps 0.00000001
using namespace std;
int num[400000];
struct node
{
    int l,r;
    long long nsum;
    long long add;
}seg[400000];
void Build(int i,int l,int r)
{
    seg[i].l=l;
    seg[i].r=r;
    seg[i].add=0;
    if(l==r)
    {
        seg[i].nsum=num[l];
        return ;
    }
    int mid=l+((r-l)>>1);
    Build(i<<1,l,mid);
    Build(i<<1|1,mid+1,r);
    seg[i].nsum=seg[i<<1].nsum+seg[i<<1|1].nsum;
}
void Add(int i,int l,int r,long long c)
{
    if(seg[i].l==l&&seg[i].r==r)
    {
        seg[i].add+=c;
        return ;
    }
    seg[i].nsum+=c*(r-l+1);
    int mid=(seg[i].l+seg[i].r)>>1;
    if(r<=mid) Add(i<<1,l,r,c);
    else if(l>mid) Add(i<<1|1,l,r,c);
    else
    {
        Add(i<<1,l,mid,c);
        Add(i<<1|1,mid+1,r,c);
    }
}
long long Query(int i,int a,int b)//查询a-b的总和
{
    if(seg[i].l==a&&seg[i].r==b)
    {
        return seg[i].nsum+(b-a+1)*seg[i].add;
    }
    seg[i].nsum+=(seg[i].r-seg[i].l+1)*seg[i].add;
    int mid=(seg[i].l+seg[i].r)>>1;
    Add(i<<1,seg[i].l,mid,seg[i].add);
    Add(i<<1|1,mid+1,seg[i].r,seg[i].add);
    seg[i].add=0;
    if(b<=mid)  return Query(i<<1,a,b);
    else if(a>mid)  return Query(i<<1|1,a,b);
    else return Query(i<<1,a,mid)+Query(i<<1|1,mid+1,b);
}
int main()
{
    int n,m,i;
    int a,b,c;
    char s[10];
    while(~scanf("%d%d",&n,&m))
    {
        for(i=1;i<=n;i++)
        {
            scanf("%d",&num[i]);
        }
        Build(1,1,n);
        for(i=0;i<m;i++)
        {
            scanf("%s",s);
            if(s[0]=='C')
            {
                scanf("%d%d%d",&a,&b,&c);
                Add(1,a,b,c);
            }
            else
            {
                scanf("%d%d",&a,&b);
                cout<<Query(1,a,b)<<endl;
            }
        }
    }
    return 0;
}

其实线段树看网上讲的再多，不如顺着代码来一便的快！

bless295
关注
  • 3
    点赞
  • 踩
  • 5
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
专栏目录
Printing Floating-Point Numbers Quickly and Accurately with Integers - 2010 (dtoa-pldi2010)-计算机科学
04-22
Printing Floating-Point Numbers Quickly and Accurately with IntegersFlorian Loitsch Inria Sophia Antipolis2004 Route des Lucioles - BP 93 - 06902 Sophia Antipolis Cedex florian.loitsch@inria....
可能是最好的A*算法简介，没有之一！
weixin_43619346的博客
08-17 1662
关于A∗A^*A∗算法，中文资源有点少，能够通俗解释A∗A^*A∗算法的也很少，之前比较不错的是这里的英文博客，斯坦福大佬十几年前写的博客，也有一些中文翻译。然后我去大佬博客主页翻了翻，看到了最近的一篇关于A∗A^*A∗的文章，详情请看这里的原文连接。本人粗略搜了一下没发现有对应的中文翻译。因此在这里翻译一份，也可以算作是一个重新的学习吧！由于本人水平有限，难免有纰漏，欢迎随时交流指出，转载请注明出处。 写在前面 本文翻译自 Red Blob Games，强烈建议看一看大佬的原博客，有很多互动的可视化演示，
‘A’与”A”的区别
电脑管家
12-04 1万+
‘A’是单独的字符 “A”则代表了一个字符的数组，该数组包含两个字符：一个是字母A，另一个是空字符。
a += a-= a*a
子衿_青青的博客
10-18 3万+
int a = 3; a += a-= a*a; 最终a的值是多少？ 这要分语言了。 1）在c语言中，结果是-12。原因是先算a*a（结果为9，此时a的值没变，还是3）；然后算 a-=a*a，等效于 a = a - 9;（结果为-6，因为赋值符号，此时a的值为-6，发生了改变）；最后算a = a + a ，结果为-12 2）在JAVA中，结果为-3。好像中间的那次赋值对它没影响，
北京大学OpenJudge 3243 A Simple Problem with Integers （线段树区间求和）
qq_37618760的博客
08-24 181
描述 You have N integers, A1, A2, ... , AN. You need to deal with two kinds of operations. One type of operation is to add some given number to each number in a given interval. The other is to ask for ...
20180521 数组时间复杂度
weixin_30773135的博客
05-21 174
A为整数数组， N为A的数组长度，请问执行以下代码，最坏情况下的时间复杂度为____。 1 2 3 4 5 6 7 8 9 10 11 void fun(int A[],int n) { ...
Improved Division by Invariant Integers-计算机科学
04-22
Niels Möller and Torbj̈orn GranlundAbstract—This paper considers the problem of dividing a two-word integer by a single-word integer, together with a few extensions and applications. Due to lack ...
LanguageDay:Objective-C 语言简介
06-28
语言日Objective-C 语言简介##什么是计算机程序？... ####Integers 整数是从 -4294967295 到 4294967295（或 -2^32 -1 到 2^32 -1）的正数或负数。 您可以像这样创建一个名为“myInteger”的整数。 int myInte
Palindrome Sub-Array
08-28
Palindrome Sub-Array,Problem Description 　A palindrome sequence is a sequence which is as same as its reversed order. For example, 1 2 3 2 1 is a palindrome sequence, but 1 2 3 2 2 is not. Given a 2-...
numbers-to-words:https
04-27
Numbers_to_words ... 规模短。 从这里缩放： : 缩放为10 ^ N，其中N从0到3003（步骤3）。 支持负数-已添加。 支持数字加上浮点数和十分之一的分数。... 大写的第一个符号，最后加点-添加。... 不需要BigInteger，数字...
1453: A Simple Problem
Warren's Blog
06-23 313
1453: A Simple Problem 时间限制: 1 Sec 内存限制: 128 MB 提交: 411 解决: 179 [提交][状态][讨论版] 题目描述 一个数字序列从起始位置开始的部分序列为：1,-1,2,-1,3,-2,5,-3,8,-5,13,-8,21,-13,34,-21,55,……问题是：输入一个值n，请将第1项到第n项的数字序列输出。输入 输入一个n值（1<=n<=88
全藏字【藏文字】的生成算法设计
qq_41139540的博客
04-17 6162
藏文字符是拼音性文字，现代藏字由30个辅音字母和4个元音符号（简称为元音）拼写组合而成，既可以前后书写，还可以上下叠加，构成二维的平面文字。现代藏字均以一个称为“基字”的辅音字母为核心，其前后添加和上下叠加，组成一个完整的字符结构。
A Simple Problem with Integers
qq_46126537的博客
01-21 177
原题传送门 代码 #include <bits/stdc++.h> using namespace std; typedef long long ll; const int N = 1e5 + 10; int n, q,t; ll a[N],sum[N], add[N]; int l[N], r[N],pos[N]; void divide() { t = sqrt((double)n); for (int i = 1; i <= t; i++) {
2020CSP初赛试题答案解析
热门推荐
zsjzliziyang的博客
10-13 4万+
目录 1.CSP-J试题（C++A卷） 2.CSP-J答案（C++A卷） 3.CSP-J试题（C++B卷） 4.CSP-J答案（C++B卷） 5.CSP-S试题（C++A卷） 6.CSP-S答案（C++A卷） 7.CSP-S试题（C++B卷） 8.CSP-S答案（C++B卷） 9.CSP-S解析 1.CSP-J试题（C++A卷） 2.CSP-J答案（C++A卷） 3.CSP-J试题（C++B卷） 4.CSP-J答...
A simple problem
七刀的专栏
07-12 780
A simple problem 题目描述 Dreamone has a lovely cat. The cat can comfort you if you are very dejected, she can also play with you if  you are very bored. So Dreamone loves her very much. Of cours
Unity Shader数学基础——矩阵
九本才的博客
01-18 2893
矩阵的定义 由m×n个数a，排成的m行n列的数表称为m行n列的矩阵，简称m×n矩阵。记作： 这m×n个数称为矩阵A的元素，简称为元，数a位于矩阵A的第i行第j列，称为矩阵A的(i,j)元，以数a为(i,j)元的矩阵可记为(a)或(a)m×n，m×n矩阵A也记作Amn。 元素是实数的矩阵称为实矩阵，元素是复数的矩阵称为复矩阵。而行数与列数都等于n的矩阵称为n阶矩阵或n阶方阵。 矩阵基本运算 加法 矩阵的加法满足下列运算律(A，B，C都是同型矩阵)： 应该注意的是只有同..
【海量数据挖掘/数据分析】 之 关联规则挖掘 Apriori 算法 （数据集、事务、频繁项集、关联规则、支持度、置信度）
最新发布
仙魁XAN
07-02 1万+
Apriori 算法 是关联规则挖掘算法 ,关联规则 反映了对象之间 相互依赖关系 ,可以通过 一个对象 的行为或属性 预测 其它对象的行为或属性;关联规则不是 因果关系, 有可能有因果关系 , 有可能没有;如 : 购买商品时 , 啤酒 与 尿布 就有关联关系 , 这两个之间肯定没有因果关系 , 有一种未知的关联关系;关联规则挖掘步骤 :① 步骤一 : 找出 支持度 ≥最小支持度阈值的 频繁项集;② 步骤二 : 根据 频繁模式 生成 满足可信度阈值的 关联规则;
第二章 软件过程 单元测试练习
yangjie_xin的博客
06-28 664
10. A公司拟开发一款应用软件，由于客户需求经常发生变化，软件开发人员拟采用把软件开发分为多个短周期，采用多次迭代的方式开发软件，以下最合适的软件过程模型是( )。5. A公司计划开发一款商业软件，由于项目风险较大，可能需要经常项目变更，故开发人员希望在软件开发过程中加入风险分析，以下最适合的软件过程模型是( )。1. A公司打算开发一个大型软件，拟采用标准化流程开发，将项目分为若干阶段，每个阶段用标准的UML语言进行编写，以下最适合的软件过程模型是（ ）。A. 计划、分析、设计、实现、调试。
系统自带输入法➋➌➍➎➏➐➑➒问题
weixin_33912638的博客
09-29 1万+
2019独角兽企业重金招聘Python工程师标准>>> ...
Cannot specify more than two integers for --resize
05-25
这个错误是因为你在使用命令行工具时，在 --resize 参数中指定了超过两个整数。请检查你的命令行语法并确保 --resize 参数只包含两个整数。例如： ``` python script.py --resize 256 256 ``` 这样就指定了一个大小为 256x256 的图片。如果你需要调整更多的维度，你需要使用其他参数或方法来实现。
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值