2021-05-18 acwing线段树，分块练习

最新推荐文章于 2024-08-15 11:25:42 发布
最新推荐文章于 2024-08-15 11:25:42 发布
阅读量107 收藏
点赞数
分类专栏： 每日三题
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_46064382/article/details/116973731
版权

243. 一个简单的整数问题2

1.题意：
2.题解：

简单个鬼鬼

3.ac代码：
#include<bits/stdc++.h>
#include<math.h>
#include<string>
using namespace std;
const int N=5e5+10;
typedef long long ll;
ll a[N];
struct Node{
    int l,r;
    ll sum,add;
}tr[N<<2];
void pushdown(int u){//下放
    if(tr[u].add){
        tr[u<<1].add+=tr[u].add;
        tr[u<<1].sum+=(ll)(tr[u<<1].r-tr[u<<1].l+1)*tr[u].add;
        tr[u<<1|1].add+=tr[u].add;
        tr[u<<1|1].sum+=(ll)(tr[u<<1|1].r-tr[u<<1|1].l+1)*tr[u].add;
        tr[u].add=0;
    }

}
void pushup(int u){
    tr[u].sum=tr[u<<1].sum+tr[u<<1|1].sum;

}
void build(int u,int l,int r){
    if(l==r){
        tr[u]={l,r,a[l],0};
    }else{
        tr[u]={l,r};
        int mid=l+r>>1;
        build(u<<1,l,mid);
        build(u<<1|1,mid+1,r);
        pushup(u);
    }
}
void change(int u,int l,int r,int x){//区间加值
    if(tr[u].l>=l&&tr[u].r<=r){
        tr[u].sum+=(ll)(tr[u].r-tr[u].l+1)*x;
        tr[u].add+=x;
    }else{
        pushdown(u);
        int mid=tr[u].l+tr[u].r>>1;
        if(l<=mid) change(u<<1,l,r,x);
        if(r>mid) change(u<<1|1,l,r,x);
        pushup(u);
    }
}
ll query(int u,int l,int r){//查找区间和
    if(tr[u].l>=l&&tr[u].r<=r){
        return tr[u].sum;
    }
    pushdown(u);
    int  mid=tr[u].l+tr[u].r>>1;
    ll sum=0;
    if(l<=mid) sum+=query(u<<1,l,r);
    if(r>mid) sum+=query(u<<1|1,l,r);
    return sum;

}
int main(){
    int n,m;
    cin>>n>>m;
    for(int i=1;i<=n;i++) cin>>a[i];
    build(1,1,n);
    string op;
    int l,r,d;
    for(int i=1;i<=m;i++){
        cin>>op>>l>>r;
        if(op=="C"){
            cin>>d;
            change(1,l,r,d);;
        }else{
            cout<<query(1,l,r)<<endl;
        }
    }
}

虽然线段树做一点不简单，但是用分块做很简单！！！，数据量1e5就可用，处理分块的总时间复杂度是0（/n)

#include<bits/stdc++.h>
using namespace std;
const int N=5e5+10;
typedef long long ll;
ll a[N];
ll sum[N];
int st[N],ed[N];
int bel[N];
int sz[N];
int mark[N];
void init(int n){
	int sq=sqrt(n);
	for(int i=1;i<=sq;i++){
		st[i]=n/sq*(i-1)+1;//st[i]表示第i号第一个元素的下标 
		ed[i]=n/sq*i;//ed[i]表示第i号元素最后一个元素的下标 
	}
	ed[sq]=n;//把剩余的一小块纳入最后一块中 
	for(int i=1;i<=sq;i++){
		for(int j=st[i];j<=ed[i];j++){
			bel[j]=i;//第j号元素归属i块 
		}
	}
	//求区间块数 
	for(int i=1;i<=sq;i++){
		sz[i]=ed[i]-st[i]+1;
	}
	//求区间和 
	for(int i=1;i<=sq;i++){
		for(int j=st[i];j<=ed[i];j++){
			sum[i]+=a[j];
		}
	}
//	for(int i=1;i<=sq;i++){
//		printf("----%d %d %d\n",st[i],ed[i],sum[i]);
//	}
}
ll query(int x,int y){
	ll s=0;
	if(bel[x]==bel[y]){
		for(int i=x;i<=y;i++){
			s+=a[i]+mark[bel[i]];
		} 
	}
	else{
		for(int i=x;i<=ed[bel[x]];i++){
			s+=a[i]+mark[bel[i]];
		}
		for(int i=st[bel[y]];i<=y;i++){
			s+=a[i]+mark[bel[i]];
		}
		for(int i=bel[x]+1;i<bel[y];i++){
			s+=sum[i]+mark[i]*sz[i];
		}
	}
	return s;
}
void change(int x,int y,int k){
	if(bel[x]==bel[y]){
		for(int i=x;i<=y;i++){
			a[i]+=k;
			sum[bel[i]]+=k;
		}
	}else{
		for(int i=x;i<=ed[bel[x]];i++){
			a[i]+=k;
			sum[bel[i]]+=k;
		}
		for(int i=st[bel[y]];i<=y;i++){
			a[i]+=k;
			sum[bel[i]]+=k;
		}
		for(int i=bel[x]+1;i<bel[y];i++){
			mark[i]+=k;//对中间块打上标记 
		}
	}
}
int main(){
	int n,m;
	cin>>n>>m;
	for(int i=1;i<=n;i++){
		cin>>a[i];
	}
	
	//分块
	 init(n);
	 int sq=sqrt(n);

	int x,y,k;
	string op;
	for(int i=1;i<=m;i++){
		cin>>op>>x>>y;
		if(op=="C"){
			cin>>k;
			change(x,y,k);
		}else{
			cout<<query(x,y)<<endl;
		}
	}
	//区间修改
	 
	
}

洛谷P2801 教主的魔法

1.题意：
2.题解：
3.ac代码：

疯掉！，各种问题，这是90分的还有一个数据T了，我真不知道错哪了

#include<bits/stdc++.h>
using namespace std;
const int N=1e6+10;
typedef long long ll;
int a[N];
ll sum[N];
int st[N],ed[N];
int bel[N];
int size[N];
int mark[N];
vector<int> v[N];
void init(int n){
    int sq=sqrt(n);
    for(int i=1;i<=sq;i++){
        st[i]=n/sq*(i-1)+1;//st[i]表示第i号第一个元素的下标
        ed[i]=n/sq*i;//ed[i]表示第i号元素最后一个元素的下标
    }
    ed[sq]=n;//把剩余的一小块纳入最后一块中
    for(int i=1;i<=sq;i++){
        for(int j=st[i];j<=ed[i];j++){
            bel[j]=i;//第j号元素归属i块
            v[i].push_back(a[j]);
        }
        sort(v[i].begin(),v[i].end());
    }
    //求区间块数
    for(int i=1;i<=sq;i++){
        size[i]=ed[i]-st[i]+1;
    }

    //求区间和
//	for(int i=1;i<=sq;i++){
//		for(int j=st[i];j<=ed[i];j++){
//			sum[i]+=a[j];
//
//		}
//	}
//	for(int i=1;i<=sq;i++){
//		printf("----%d %d %d\n",st[i],ed[i],sum[i]);
//	}
}
void update(int b) // 更新排序后的数组
{
    for (int i = 0; i <=size[b]; ++i)
        v[b][i]=a[st[b]+i];
    sort(v[b].begin(), v[b].end());
}
int query(int x,int y,ll c){
    int s=0;
    if(bel[x]==bel[y]){
        for(int i=x;i<=y;i++){
            if(a[i]+mark[bel[x]]>=c) s++;
        }
    }else{
        for(int i=x;i<=ed[bel[x]];i++){
            if(a[i]+mark[bel[x]]>=c) s++;
        }
        for(int i=st[bel[y]];i<=y;i++){
            if(a[i]+mark[bel[y]]>=c) s++;
        }
        for(int i=bel[x]+1;i<bel[y];i++){
            s+=v[i].end()-lower_bound(v[i].begin(),v[i].end(),c-mark[i]);
        }
    }
    return s;
}
void change(int x,int y,int k){
    if(bel[x]==bel[y]){
        for(int i=x;i<=y;i++){
            a[i]+=k;
            //	sum[bel[i]]+=k;
        }
        update(bel[x]);
    }else{
        for(int i=x;i<=ed[bel[x]];i++){
            a[i]+=k;
            //	sum[bel[i]]+=k;
        }
        update(bel[x]);
        for(int i=st[bel[y]];i<=y;i++){
            a[i]+=k;
            //	sum[bel[i]]+=k;
        }
        update(bel[y]);
        for(int i=bel[x]+1;i<bel[y];i++){
            mark[i]+=k;//对中间块打上标记
        }
    }

}
int main(){
    int n,m;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
        scanf("%d",&a[i]);
    }

    //分块
    init(n);

    int x,y,k,c;
    char op;
    for(int i=1;i<=m;i++){
        scanf(" %c %d %d %d",&op,&x,&y,&c);
        if(op=='M'){
            change(x,y,c);
        }else{
            printf("%d\n",query(x,y,c));
           // cout<<query(x,y,c)<<endl;
        }
    }
    //区间修改


}

这是满分的

#include<bits/stdc++.h>
using namespace std;
const int N=5e5+10;
typedef long long ll;
const int MAXN = 1000005, SQ = 1005;
int st[SQ], ed[SQ], size[SQ], bel[MAXN];
void init_block(int n) // 初始化
{
    int sq = sqrt(n);
    for (int i = 1; i <= sq; ++i)
    {
        st[i] = n / sq * (i - 1) + 1;
        ed[i] = n / sq * i;
    }
    ed[sq] = n;
    for (int i = 1; i <= sq; ++i)
        for (int j = st[i]; j <= ed[i]; ++j)
            bel[j] = i;
    for (int i = 1; i <= sq; ++i)
        size[i] = ed[i] - st[i] + 1;
}
int A[MAXN], mark[SQ];
vector<int> v[SQ]; // 这里用vector存排序后的数组
void update(int b) // 更新排序后的数组
{
    for (int i = 0; i <= size[b]; ++i)
        v[b][i] = A[st[b] + i];
    sort(v[b].begin(), v[b].end());
}
int main()
{
    int n,m;
    scanf("%d %d",&n,&m);
    int sq = sqrt(n);
    init_block(n);
    for (int i = 1; i <= n; ++i)
        scanf("%d",&A[i]);
    for (int i = 1; i <= sq; ++i)
        for (int j = st[i]; j <= ed[i]; ++j)
            v[i].push_back(A[j]);
    for (int i = 1; i <= sq; ++i)
        sort(v[i].begin(), v[i].end());
    while (m--)
    {
        char o;
        scanf(" %c", &o);
        int l,r,k;
        scanf("%d%d%d",&l,&r,&k);
       // int l = read(), r = read(), k = read();
        if (o == 'M')
        {
            if (bel[l] == bel[r]) // 如果同属一块直接暴力
            {
                for (int i = l; i <= r; ++i)
                    A[i] += k;
                update(bel[l]);
                continue;
            }
            for (int i = l; i <= ed[bel[l]]; ++i) // 对零散块暴力处理
                A[i] += k;
            for (int i = st[bel[r]]; i <= r; ++i)
                A[i] += k;
            update(bel[l]);
            update(bel[r]);
            for (int i = bel[l] + 1; i < bel[r]; ++i) // 打上标记
                mark[i] += k;
        }
        else
        {
            int tot = 0;
            if (bel[l] == bel[r])
            {
                for (int i = l; i <= r; ++i)
                    if (A[i] + mark[bel[l]] >= k)
                        tot++;
                printf("%d\n", tot);
                continue;
            }
            for (int i = l; i <= ed[bel[l]]; ++i)
                if (A[i] + mark[bel[l]] >= k)
                    tot++;
            for (int i = st[bel[r]]; i <= r; ++i)
                if (A[i] + mark[bel[r]] >= k)
                    tot++;
            // 二分查找k-mark[i]的位置，因为整块都加上了mark[i]其实就相当于k减去mark[i]
            for (int i = bel[l] + 1; i < bel[r]; ++i) 
                tot += v[i].end() - lower_bound(v[i].begin(), v[i].end(), k - mark[i]);
            printf("%d\n", tot);
        }
    }
    return 0;
}

SPOJ DQUERY - D-query

1.题意：
2.题解：

依旧是外网题不知道过了没

3.ac代码：
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e6+10;
int sq;
struct query{
	int l,r,id;
	bool operator<(const query&o) const{
		if(l/sq!=o.l/sq){
			return l<o.l;
		}
		if(l/sq&1){
			return r<o.r;
		}
		return r>o.r;
	}
}Q[N];
int cur;
int a[N],ans[N],cnt[N];
inline void add(int p){
	if(cnt[a[p]]==0){
		cur++;
	}
	cnt[a[p]]++;
}
inline void del(int p){
	cnt[a[p]]--;
	if(cnt[a[p]]==0)
		cur--;
}
int main(){
	int n;
	cin>>n;
	sq=sqrt(n);
	for(int i=1;i<=n;i++){
		scanf("%d",&a[i]);
	}
	int q;
	cin>>q;
	for(int i=1;i<=q;i++){
		cin>>Q[i].l>>Q[i].r;
		Q[i].id=i;
	}
	sort(Q+1,Q+q+1);
	int l=1,r=0;
	for(int i=1;i<=q;i++){
		while(l>Q[i].l){
			add(--l);
		}
		while(r<Q[i].r){
			add(++r);
		}
		while(l<Q[i].l){
			del(l++);
		}
		while(r>Q[i].r){
			del(r--);
		}
		ans[Q[i].id]=cur;
	}
	for(int i=1;i<=q;i++){
		cout<<ans[i]<<endl;
	}
}
数据结构——分块入门—＞扩展
沫惘的博客
02-24 828
数据结构——分块入门前言分块详解Loj 6277. 数列分块入门 1传送门题目分析CodeLoj 6278. 数列分块入门 2传送门题目分析Code练习预留扩展 前言 \quad首先，我们来考虑这样一个模型：有一段连续的序列a[1]∼a[n]a[1]\sim a[n]a[1]∼a[n]，现在我们需要执行几类操作： 求出其中一段区间的和  \,这时我们可以使用前缀和，记录a[1]∼a[i]的和为sum[i]a[1]\sim a[i]的和为sum[i]a[1]∼a[i]的和为sum[i]，显然sum[i]=s
【机试刷题经验】清北保研机试刷题指南
程序员阿浪
09-21 2423
这篇博客主要是整理了一下我在保研过程中看到的机试经验分享。 推荐书目：《王道考研机试指南》 清华： 主要考察动态规划，数据结构，算法： a. 动态规划考题的主要考点包括：背包问题、区间DP、环形DP、判定性DP、棋盘分割、最长公共子序列、最长上升子序列、编辑距离。 b. 数据结构考题的主要考点包括：队列、链表、单调队列、二分、哈希表、快排、归并排、堆排、并查集、霍夫曼树、线段树、二叉树。...
1000：入门测试题目[不一样的题解][85种写法]【A+B问题】
lljloimjo的博客
02-20 1342
算法一、DFS一号 #include using namespace std; int n = 2, a[5], s; int dfs(int x, int sum) { if (x > n) return sum; int i = dfs(x + 1, sum); int j = dfs(x + 1, sum + a[x]); if (i == s) return i; if (j == s) return j; return
算法笔记(个人用)(不定期更新)
cuber-lotus的博客
03-23 2694
文章目录前言说明Note文档编辑相关语言特性相关 (C/C++)输出数组最后一个元素后面无空格二分查找函数rand() 范围算法相关排序区间求和题单基础数学素数素数判断素数筛朴素筛选 (埃式筛)朴素筛选优化线性筛 (欧拉筛)闰年判断蔡勒（Zeller）公式GCD10转R进制卡特兰数完美数分解因数欧几里得-欧拉定理数根[1~n] 中完全平方数的个数计算Π pia二分查找顺序序列找唯一值逼近类二分容斥原理位运算整形存储大小intlong long`__builtin_函数`取模自动取模工具类除法取模判断奇偶统计
算法の笔记
最新发布
weixin_45419526的博客
08-15 760
包括基本算法、数据结构、图论、动态规划、数学知识、字符串
蓝桥杯_贪心
violet的博客
04-06 492
#include <iostream> using namespace std; 由数据范围反推算法复杂度以及算法内容 链接：https://www.acwing.com/blog/content/32/ 一般ACM或者笔试题的时间限制是1秒或2秒。 在这种情况下，C++代码中的操作次数控制在 107～108为最佳。 n≤30, 指数级别, dfs+剪枝，状态压缩dp n≤100 => O(n3)，floyd，dp，高斯消元 n≤1000 => O(n2)，O(n2logn)，
2021-2022收藏资料练习-7.1与三角形有关的线段习题.doc
09-26
2021年09月16日
2021-2022收藏资料七年级数学--与三角形有关的线段.doc
09-26
2021年09月16日
2021-2022收藏资料复习精品 与三角形有关的线段.doc
09-25
- 练习题1和2：计算等腰三角形的周长，并判断哪些线段组合能构成三角形。 - 练习题3：分析给定线段长度是否能满足三角形的条件，从而决定能否组成三角形。 通过以上讲解，我们可以看到，本节课主要涵盖了三角形的...
线段树1代码资源（普通/指针线段树,分块,树状数组）
01-31
本资源包包含了多种线段树的实现方式，包括普通线段树、指针线段树以及分块技术与树状数组的结合，适用于C++编程环境。 1. **普通线段树**：普通线段树是线段树的基础形态，通过数组存储节点，每个节点代表一个区间...
2021 ICPC Asia EC网络预选赛I、II
weixin_46064382的博客
10-06 535
1.题意： 2.题解： 3.ac代码： 在这里插入代码片 1.题意： 2.题解： 3.ac代码： 在这里插入代码片 1.题意： 2.题解： 3.ac代码： 在这里插入代码片 1.题意： 2.题解： 3.ac代码： 在这里插入代码片 1.题意： 2.题解： 3.ac代码： 在这里插入代码片 1.题意： 2.题解： 3.ac代码： 在这里插入代码片 1.题意： 2.题解： 3.ac代码： 在这里插入代码片 1.题意： 2.题解： 3.ac代码： 在这里插入代码片 1.题意： 2.
2021-01-13Winter Camp 2021 Training 3训练题解
weixin_46064382的博客
01-16 270
contest 一.A. Phoenix and Balance 1.大意： 有n个数，分别是从2到2^n，将他们随机分成2堆（每堆都是n/2个），保证n是偶数，让两堆n/2个数的和的差最小 2.题解： 找规律，发现分成一堆为：2n/2,2n/2+1,…,2n-1,(共n/2个数）， 另一堆为·：21,22,…2n/2-1,2n, 这种情况两边和相差最小 推公式得：2n/2+1-2 用pow会T，用了快速幂 看了别人的方法，发现也可以是前一个数*2+2，实质差不多 3.ac代码： #include
2021-05-09 acwing算法基础课 求组合数
weixin_46064382的博客
06-01 259
AcWing 885. 求组合数 I 1.题意： 2.题解： 3.ac代码： 在这里插入代码片 AcWing 886. 求组合数 II 1.题意： 2.题解： 3.ac代码： 在这里插入代码片 AcWing 887. 求组合数 III 1.题意： 2.题解： 3.ac代码： 在这里插入代码片 1.题意： 2.题解： 3.ac代码： 在这里插入代码片 1.题意： 2.题解： 3.ac代码： 在这里插入代码片 1.题意： 2.题解： 3.ac代码： 在这里插入代码片 1.题意： 2.题解：
2021-05-24异或线段树题1
weixin_46064382的博客
05-28 249
在 [first, last) 区域内查找 lower_bound(,x) 第一个>=x int a[5] = { 1,2,3,4,5 }; //从 a 数组中找到第一个不小于 3 的元素 int *p = lower_bound(a, a + 5, 3); *p=3 upper bound(,x) 第一个>x int a[5] = { 1,2,3,4,5 }; //从 a 数组中找到第一个大于 3 的元素 int *p = uppe
2021-01-19Winter Camp 2021 Training 6训练
weixin_46064382的博客
02-07 242
一. 1.题面： 2.题解： 3.ac代码： #include<stdio.h> #include<bits/stdc++.h> #include<iostream> #include<string.h> #include<stdlib.h> #define N 300010 using namespace std; typedef long long ll; int dp[N]; int a[N]; int main(){ int t,n,i,
2021-05-31 2021牛客寒假算法基础集训营1
weixin_46064382的博客
06-02 198
I 限制不互素对的排列 1.题意： 2.题解： 3.ac代码： #include<bits/stdc++.h> const int N=1e5+10; using namespace std; int a[N]; int tot; int main(){ int n,k; cin>>n>>k; if(n<6&&k==n/2){ cout<<-1<<endl; retur
2021-04-25 acwing算法基础课 前缀和与差分 795，796，797，798
weixin_46064382的博客
04-25 197
795 前缀和 1.题意： 2.题解： 3.ac代码： #include<bits/stdc++.h> using namespace std; int n,m,l,r; const int N=1e5+10; typedef long long ll; int a[N],d[N]; int main(){ cin>>n>>m; for(int i=1;i<=n;i++){ cin>>a[i]; d[i]=
2021-05-12 acwing算法基础课 快速幂，扩展欧几里得，中国剩余定理，高斯消元
weixin_46064382的博客
06-01 194
AcWing 875. 快速幂 1.题意： 2.题解： 3.ac代码： #include<bits/stdc++.h> using namespace std; const int N=1e6+10; typedef long long ll; ll quick_pow(ll a,ll k,ll p){ ll res=1; while(k){ if(k&1) res=res*a%p; k>>=1; a=
2021-2022年计算机图形学考试精华题目解析
本资源是一份针对2021-2022年计算机图形学考试的精选题目集，包含了五个主要的考察点。首先，题目要求考生实现全象限整型Bresenham直线段生成算法的C语言版本，该算法用于高效地在屏幕上描绘像素化线条，通过控制...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值