A Simple Problem with Integers

本文介绍了一种利用线段树解决区间更新和查询问题的方法。通过延迟更新机制提高效率,实现了对序列进行区间加法操作及求区间和的功能。

线段树之区间求和


这是我的第一篇博客,希望大家能看的清楚明白。
先声明一下我的代码风格吧!我个人喜欢一般先把main()全部写完,然后再定义内部调用的各个函数,先把整个程序的框架写好,再去实现具体的细节,希望不会影响到大家的阅读。头文件就是常用的头文件,预处理我会事先说明的。
这也是我做的第一个区间更新的题目,感觉好辛苦的说。


Question:

Description:
给出了一个序列,你需要处理如下两种询问。
“C a b c”表示给[a, b]区间中的值全部增加c (-10000 ≤ c ≤ 10000)。
“Q a b” 询问[a, b]区间中所有值的和。
Input:
第一行包含两个整数N、 Q(1 ≤ N,Q ≤ 100000)。
第二行包含n个整数,表示初始的序列A (-1000000000 ≤ A[i] ≤ 1000000000)。
接下来Q行询问,格式如题目描述。
Output:
对于每一个Q开头的询问,你需要输出相应的答案,每个答案一行。
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


Answer:

数据预处理一般是数据范围的指定:

#define MAX 100010

自定义数据类型:对于线段树,我只会用结构体数组来表示,结构体包含的是区间范围,以及我们要维护的对象和所需要的标记。
区间更新不同于单点更新,单点更新是即时的,每一次更新都会更新整棵线段树,而区间更新的思想是延迟更新,即在我们需要查询的时候才更新线段树相应的节点,用不到的话自然就不用去维护了,采用这种思想来提高效率。因此在这个题我们需要设立一个延迟标记,表示这个区间是否需要被更新。
所以结构体的设置如下所示(注意数据范围):

struct Tree
{
    int left,right;//区间上下限
    long long delay,sum;//延迟标记和要维护的对象
}tr[MAXN*4];

接下来是main():

int main()
{
    int N,Q;
    while(scanf("%d%d",&N,&Q)==2)//支持多组输入
    {
        int i=0;
        build(1,1,N);//build()用于线段树的初始化,返回值为空,参数需要线段树的根节点(id=1)以及线段树的总长度,我们要在这个函数里面输入我们的初始值哦,也就是每一个叶子节点
        while(Q--)
        {
            char order[2];
            scanf("%s",order);
            if(order[0]=='Q')
            {
                int l,r;
                scanf("%d%d",&l,&r);
                printf("%I64d\n",query(1,l,r));//query()需要返回区间[l,r]内的和,参数是查询的起始节点(id=1)以及需要查询的区间[l,r]
            }
            if(order[0]=='C')
            {
                int l,r;
                long long val;
                scanf("%d%d%I64d",&l,&r,&val);
                updata(1,l,r,val);//updata()用于数据的区间更新,无返回值,参数是查询的根节点(id=1),需要更新的区间[l,r]以及需要被维护的对象的维护方式,这里是对区间的每个数都加上val
            }
        }
    }
    return 0;
}

线段树的初始化函数:

//build()用于线段树的初始化,返回值为空,参数需要线段树的根节点(id)以及线段树的总长度,我们要在这个函数里面输入我们的初始值哦,也就是每一个叶子节点
void build(int id,int l,int r)
{
    tr[id].left=l;
    tr[id].right=r;
    tr[id].delay=0;//初始化时每一个都不需要被标记的
    if(l==r)
        scanf("%I64d",&tr[id].sum);//递归到叶子节点
    else
    {
        int mid=(l+r)/2;
        build(id*2,l,mid);//建立左子树
        build(id*2+1,mid+1,r);//建立右子树
        tr[id].sum=tr[id*2].sum+tr[id*2+1].sum;//每次记得更新父节点的sum值哦
    }
}

线段树数据区间更新函数:

//updata()用于数据的区间更新,无返回值,参数是查询的根节点(id),需要更新的区间[l,r]以及需要被维护的对象的维护方式,这里是对区间的每个数都加上val
void updata(int id,int l,int r,long long val)
{
    if(l<=tr[id].left&&r>=tr[id].right)//这个区间这次需要被更新
    {
        tr[id].delay+=val;//表示这个区间需要被更新的方式,但暂时不更新
        tr[id].sum+=val*(tr[id].right-tr[id].left+1);//更新这个区间的和
    }
    else//这个是搜寻需要更新的区间
    {
        if(tr[id].delay!=0)//搜寻到这里发现这个区间的子树还没更新呢,没办法继续搜寻,所以我们先更新它的子树
        {
            pushdown(id);//pushdown()用于更新左右子树,无返回值,参数只需要当前的根节点(id)
            tr[id].delay=0;//记得把延迟标记重新置为0哦,因为它的左右子树已经更新过了
        }
        //我们接着搜寻需要被更新的区间
        int mid=(tr[id].left+tr[id].right)/2;
        if(r<=mid)
        {
            updata(id*2,l,r,val);
        }
        else if(l>mid)
        {
            updata(id*2+1,l,r,val);
        }
        else
        {
            updata(id*2,l,r,val);
            updata(id*2+1,l,r,val);
        }
        tr[id].sum=tr[id*2].sum+tr[id*2+1].sum;//每次更新,父节点也会跟着改变嘛,所以需要回溯更新父节点的的维护对象,这里是区间和
    }
}

这里是查询函数:

//query()需要返回区间[l,r]内的和,参数是查询的起始节点(id)以及需要查询的区间[l,r]
long long query(int id,int l,int r)
{
    if(l<=tr[id].left&&r>=tr[id].right)//这个区间包含在我们要查询的区间,所以返回它的和
    {
        return tr[id].sum;
    }
    else//搜寻我们要查询的区间
    {
        if(tr[id].delay!=0)//搜寻到这里发现之前被标记过,然而我们要继续搜寻它的子树,所以我们需要先把它的子树更新
        {
            pushdown(id);
            tr[id].delay=0;
        }
        int mid=(tr[id].left+tr[id].right)/2;
        if(r<=mid)
            return query(id*2,l,r);
        else if(l>mid)
            return query(id*2+1,l,r);
        else
        {
            return query(id*2,l,r)+query(id*2+1,l,r);//两个子树都需要被搜寻所以两个子树的返回值都需要加起来
        }
    }
}

最后一个需要实现的函数:

//pushdown()用于更新左右子树,无返回值,参数只需要当前的根节点(id)
void pushdown(int id)
{
    //把延迟标记传递给左右子树
    tr[id*2].delay+=tr[id].delay;
    tr[id*2+1].delay+=tr[id].delay;
    //维护左右子树的维护对象,这里是维护区间和
    tr[id*2].sum+=tr[id].delay*(tr[id*2].right-tr[id*2].left+1);
    tr[id*2+1].sum+=tr[id].delay*(tr[id*2+1].right-tr[id*2+1].left+1);
}

这样整个题就完全解决了呢!
时空复杂度就像这样:

MemoryTimeLenth
6556(kb)2657(ms)3818(Bytes)

感觉时空复杂度还能降低,但是进一步的优化我却没什么思路,希望那位大佬看到后可以指教一下。


不知道大家对我的风格还习惯吗?写得不好,欢迎各位大佬们来指点我的不足,如果刚入坑的小伙伴们觉得还有不理解的地方,也请指出来,在能力范围内,我很乐意为你们解答。

<think>好的,我现在需要帮用户翻译一个关于图论的题目。首先,用户之前让我翻译过汉明距离的问题,现在又来一个更复杂的关于森林的最小删边问题。用户可能是正在准备编程竞赛或者学习算法,需要准确理解题目要求。 首先,我要仔细阅读用户提供的英文内容,确保理解每一个术语和条件。题目要求的是将给定的无向图变成森林所需要删除的最小边数。森林的定义是不含任何环的图,也就是所有连通分量都是树的结构。树的性质是边数等于顶点数减一,所以每个连通分量的边数超过顶点数减一时,就需要删除多余的边。 接下来,我需要将问题准确翻译成中文,同时注意保持技术术语的正确性。例如,"simple undirected graph"要翻译为“简单无向图”,"forest"对应“森林”。还要注意数学符号和格式,比如变量用斜体,公式使用LaTeX,并且符合之前用户指定的系统级指令中的格式要求,比如行内公式用$...$,独立公式用$$...$$。 然后,题目中的约束条件也需要准确处理。比如顶点数N的范围是1到2e5,边数M的范围是0到min(N(N-1)/2, 2e5),并且每条边的两个顶点u_i和v_i满足1≤u_i<v_i≤N。这些约束条件需要清晰列出,方便解题时考虑算法的时间复杂度和数据结构的选择。 另外,用户可能希望翻译后的内容结构清晰,包含问题描述、关键术语解释、示例和实现建议,就像之前处理汉明距离的问题一样。因此,我需要组织内容,分点说明,可能还需要给出解题思路,比如使用并查集(Union-Find)来检测环,从而计算需要删除的边数。 需要注意用户提到的输入约束中的格式问题,比如原英文中的数学表达式可能有排版错误,例如“0≤M≤min( 2 N(N−1) ,2×10 5 )”应该正确翻译为“0 ≤ M ≤ min( N(N-1)/2 , 2×10^5 )”。这里要确保所有的数学公式都正确无误,使用合适的LaTeX语法。 最后,确保整个翻译符合用户之前提供的格式要求,比如使用中文回答,正确应用行内和独立公式的格式,结构清晰,分点明确。可能还需要添加示例和实现建议,帮助用户更好地理解如何解决这个问题。</think>**题目描述** 给定一个简单无向图,包含$N$个顶点(编号为$1$到$N$)和$M$条边。其中第$i$条边连接顶点$u_i$和$v_i$。 求需要删除的**最小边数**,使得剩余的图成为一个**森林**。 --- ### 关键术语说明 1. **森林** 定义为不含任何**环**的简单无向图。其等价于由若干棵树(连通无环子图)组成的集合。 - 树的性质:含$k$个顶点的树有且仅有$k-1$条边。 2. **简单无向图** - 无自环(边连接相同顶点) - 无重复边(任意两顶点间最多一条边) --- ### 输入约束 - $N$:顶点数,满足$1 \leq N \leq 2 \times 10^5$ - $M$:边数,满足$0 \leq M \leq \min\left(\frac{N(N-1)}{2}, 2 \times 10^5\right)$ - 边定义:$1 \leq u_i < v_i \leq N$ - 输入保证图是简单图,且所有值为整数。 --- ### 解题思路 1. **森林的边数限制** 若图最终成为森林,其总边数为$N - C$,其中$C$为连通分量个数。因此需删除的边数为: $$ \text{最小删除边数} = M - (N - C) $$ 2. **连通分量检测** 使用**并查集(Union-Find)**算法统计连通分量数$C$,并记录合并失败次数(即检测到的环数)。 - 每检测到一次合并失败(形成环),需删除至少一条边。 --- ### 示例 **输入** $N=4$, $M=5$,边集合: $$ \{(1,2), (1,3), (2,3), (2,4), (3,4)\} $$ **输出** $2$(原图含2个环,需删除2条边) --- ### 实现代码(Python) ```python def min_edges_to_forest(): import sys input = sys.stdin.read data = input().split() N = int(data[0]) M = int(data[1]) parent = list(range(N+1)) def find(u): while parent[u] != u: parent[u] = parent[parent[u]] u = parent[u] return u cycles = 0 ptr = 2 for _ in range(M): u = int(data[ptr]) v = int(data[ptr+1]) ptr += 2 pu, pv = find(u), find(v) if pu == pv: cycles += 1 else: parent[pv] = pu # 计算连通分量数C C = sum(1 for i in range(1, N+1) if parent[i] == i) return M - (N - C) # 等价于 cycles + (C_initial - C) ``` --- ### 时间复杂度 - **并查集操作**:近似$O(M \cdot \alpha(N))$,其中$\alpha$为阿克曼函数的反函数,可视为常数。 - 总复杂度:$O(N + M)$,适用于大规模输入。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值