只增莫队学习_4241历史研究_JOI 2013~2014 春季training合宿 竞技1 By PoPoQQQ

最新推荐文章于 2024-04-17 15:51:03 发布
最新推荐文章于 2024-04-17 15:51:03 发布
阅读量370 收藏
点赞数
分类专栏： 莫队/分块 ACM算法之旅
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_39599067/article/details/80747164
版权

（有任何问题欢迎留言或私聊 && 欢迎交流讨论哦

当莫队删数困难时，考虑“只增莫队”

 增数处理很容易，删数处理却较难处理。
 对于[L, R]区间的询问：假设L所处块的右端点为Rs，从Rs向两边增数！
 当L所处块相同时，R是只增的。L的变化肯定也是从Rs增数过来的。
 你可能会说，虽然说L在同一块内，但是他们大小是没有规律的。没错，所以在一次询问完成后，要把L还原到Rs，这样下一次还会是增数处理。

复杂度玄学： O ( 1 e 5 能 过 ) O(1e5能过) O(1e5能过)

AC代码

8002ms

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N = 1e6+7;
const int INF = 0x3f3f3f3f;
const int mod = 1e9+7;
int n,m;
LL ans;
int ar[N],belong[N],br[N],now[N];
LL Ans[N];
int cnt[N],in_left[N];
struct lp{
    int l,r,id;
}cw[N];
bool cmp(const lp &a,const lp &b){
    if(belong[a.l]==belong[b.l]){
        return a.r<b.r;
    }
    return belong[a.l]<belong[b.l];
}
void update_left(int x,int f){
    cnt[now[x]]++;
    if((cnt[now[x]])*1LL*ar[x]>ans){
        ans=(cnt[now[x]])*1LL*ar[x];
    }
}
int main(){
    scanf("%d%d",&n,&m);
    int block=(int)sqrt(n*1.0);
    for(int i=1;i<=n;++i){
        scanf("%d",&ar[i]);
        br[i]=ar[i];
        belong[i]=(i-1)/block+1;
    }
    sort(br+1,br+1+n);
    int k=1;
    for(int i=2;i<=n;++i){
        if(br[i]!=br[i-1]){
            br[++k]=br[i];
        }
    }
    for(int i=1;i<=n;++i){
        now[i]=lower_bound(br+1,br+1+n,ar[i])-br;
    }
    for(int i=1;i<=m;++i){
        scanf("%d%d",&cw[i].l,&cw[i].r);
        cw[i].id=i;
    }
    memset(cnt,0,sizeof(cnt));
    sort(cw+1,cw+m+1,cmp);
    int L=1,R=0,limit;
    LL last_ans=0;
    ans=0;cw[0].l=0;belong[0]=0;
    for(int i=1;i<=m;++i){
        if(belong[cw[i-1].l]!=belong[cw[i].l]){
            memset(cnt,0,sizeof(cnt));
            limit=belong[cw[i].l]*block;
            R=limit;
            ans=last_ans=0;
        }
        L=min(limit+1,cw[i].r+1);
        while(R<cw[i].r){
            R++;ans=max((++cnt[now[R]])*1LL*ar[R],ans);
        }
        last_ans=ans;
        while(L>cw[i].l)update_left(--L,1);
        Ans[cw[i].id]=ans;
        for(int j=L;j<=cw[i].r&&j<=limit;++j){
            cnt[now[j]]--;
        }
        ans=last_ans;
    }
    for(int i=1;i<=m;++i){
        printf("%lld\n",Ans[i] );
    }
    return 0;
}

题目

Description
IOI国历史研究的第一人——JOI教授，最近获得了一份被认为是古代IOI国的住民写下的日记。JOI教授为了通过这份日记来研究古代IOI国的生活，开始着手调查日记中记载的事件。
日记中记录了连续N天发生的时间，大约每天发生一件。
事件有种类之分。第i天(1<=i<=N)发生的事件的种类用一个整数Xi表示，Xi越大，事件的规模就越大。
JOI教授决定用如下的方法分析这些日记：

  1. 选择日记中连续的一些天作为分析的时间段
  2. 事件种类t的重要度为t*(这段时间内重要度为t的事件数)
  3. 计算出所有事件种类的重要度，输出其中的最大值

现在你被要求制作一个帮助教授分析的程序，每次给出分析的区间，你需要输出重要度的最大值。
Input
第一行两个空格分隔的整数N和Q，表示日记一共记录了N天，询问有Q次。
接下来一行N个空格分隔的整数X1…XN，Xi表示第i天发生的事件的种类
接下来Q行，第i行(1<=i<=Q)有两个空格分隔整数Ai和Bi，表示第i次询问的区间为[Ai,Bi]。
Output
输出Q行，第i行(1<=i<=Q)一个整数，表示第i次询问的最大重要度
Sample Input
5 5
9 8 7 8 9
1 2
3 4
4 4
1 4
2 4
Sample Output
9
8
8
16
16
Hint
1<=N<=10^5

1<=Q<=10^5

1<=Xi<=10^9 (1<=i<=N)

Source
JOI 2013~2014 春季training合宿 竞技1 By PoPoQQQ

以前写过的代码1：9479ms

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N = 1e6+7;
const int INF = 0x3f3f3f3f;
const int mod = 1e9+7;
int n,m,ar[N],belong[N],br[N],now[N];
LL ans,Ans[N];
int cnt[N],in_left[N];
struct lp{
    int l,r,id;
}cw[N];
bool cmp(const lp &a,const lp &b){
    if(belong[a.l]==belong[b.l]){
        return a.r<b.r;
    }
    return belong[a.l]<belong[b.l];
}
void update_left(int x,int f){
	in_left[now[x]]++;
	if((in_left[now[x]]+cnt[now[x]])*1LL*ar[x]>ans){
		ans=(in_left[now[x]]+cnt[now[x]])*1LL*ar[x];
	}
}
int main(){
	scanf("%d%d",&n,&m);
		int sz=(int)sqrt(n*1.0);
		for(int i=1;i<=n;++i){
			scanf("%d",&ar[i]);
			br[i]=ar[i];
			belong[i]=(i-1)/sz+1;
		}
		sort(br+1,br+1+n);
		int k=1;
		for(int i=2;i<=n;++i){
			if(br[i]!=br[i-1]){
				br[++k]=br[i];
			}
		}
		for(int i=1;i<=n;++i){
			now[i]=lower_bound(br+1,br+1+n,ar[i])-br;
		}
		for(int i=1;i<=m;++i){
			scanf("%d%d",&cw[i].l,&cw[i].r);
			cw[i].id=i;
		}
		sort(cw+1,cw+m+1,cmp);
		int L=1,R=0,tmd;
		LL tmp=0;
		ans=0;cw[0].l=0;belong[0]=0;
		for(int i=1;i<=m;++i){
			if(belong[cw[i-1].l]!=belong[cw[i].l]){
				memset(cnt,0,sizeof(cnt));
				tmd=belong[cw[i].l]*sz;
				R=tmd;
				ans=tmp=0;
			}
			L=min(tmd+1,cw[i].r+1);
			while(L>cw[i].l)update_left(--L,1);
			while(R<cw[i].r){
				R++;
				tmp=max((++cnt[now[R]])*1LL*ar[R],tmp);
				ans=max(ans,(cnt[now[R]]+in_left[now[R]])*1LL*ar[R]);
			}
			Ans[cw[i].id]=ans;
			for(int j=L;j<=cw[i].r&&j<=tmd;++j){
				in_left[now[j]]--;
			}
			ans=tmp;
		}
		for(int i=1;i<=m;++i){
			printf("%lld\n",Ans[i] );
		}
	return 0;
}

以前写过的代码2： 10946ms

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N = 1e6+7;
const int INF = 0x3f3f3f3f;
const int mod = 1e9+7;
int n,m,cnt[N],ar[N],br[N],now[N],belong[N],cnt_L,cnt_R;
LL ans,tmp_ans,Ans[N];
struct lp{
    int l,r,id;
}cw[N],liu[N];
bool cmp(const lp &a,const lp &b){
    if(belong[a.l]==belong[b.l]){
        return a.r<b.r;
    }
    return belong[a.l]<belong[b.l];
}
bool cmp2(const lp &a,const lp &b){
	if(a.l!=b.l)return a.l<b.l;
	return a.r<b.r;
}
void update(int x,int f){
	cnt[now[x]]+=f;
	cnt_L=min(cnt_L,now[x]);
	cnt_R=max(cnt_R,now[x]);
	if(ans<(cnt[now[x]]*1LL*ar[x])&&f==1){
		ans=(cnt[now[x]]*1LL*ar[x]);
	}
}
int main(){
	scanf("%d%d",&n,&m);
		int sz=(int)sqrt(n*1.0);
		for(int i=1;i<=n;++i){
			scanf("%d",&ar[i]);
			br[i]=ar[i];
			belong[i]=(i-1)/sz+1;
		}
		sort(br+1,br+1+n);
		int k=1;
		for(int i=2;i<=n;++i){
			if(br[i]!=br[i-1]){
				br[++k]=br[i];
			}
		}
		for(int i=1;i<=n;++i){
			now[i]=lower_bound(br+1,br+1+k,ar[i])-br;
		}
		for(int i=1;i<=m;++i){
			scanf("%d%d",&cw[i].l,&cw[i].r);
			cw[i].id=i;
		}
		sort(cw+1,cw+m+1,cmp);
		int L=0,R,liu_k=0;
		cnt_L=0,cnt_R=0;
		ans=tmp_ans=0;cw[0].l=liu[0].l=0;
		memset(cnt,0,sizeof(cnt));
		for(int i=1;i<=m;++i){
			if(belong[cw[i].l]==belong[cw[i].r]){
				liu[++liu_k].l=cw[i].l;
				liu[liu_k].r=cw[i].r;
				liu[liu_k].id=cw[i].id;
				continue;
			}
			if(belong[cw[i].l]!=belong[L]){
				cnt_L=max(cnt_L-5,0);cnt_R+=5;
				for(int j=cnt_L;j<=cnt_R;++j)cnt[j]=0;
				ans=0;R=cw[i].r;L=belong[cw[i].l]*sz;
				cnt_L=INF,cnt_R=0;
				for(int j=L;j<=R;++j)update(j,1);
				tmp_ans=ans;
			}
			ans=tmp_ans;L=belong[cw[i].l]*sz;
			while(R<cw[i].r)update(++R,1);
			tmp_ans=ans;
			while(L>cw[i].l)update(--L,1);
			Ans[cw[i].id]=ans;
			for(int j=L;j<belong[cw[i].l]*sz;++j)update(j,-1);
		}
		sort(liu+1,liu+liu_k+1,cmp2);
		L=0;
		for(int i=1;i<=liu_k;++i){
			if(liu[i].l!=L){
				ans=0;
				for(int j=0;j<=k;++j)cnt[j]=0;
				L=liu[i].l,R=liu[i].r;
				for(int j=L;j<=R;++j)update(j,1);
			}else{
				for(int j=R+1;j<=liu[i].r;++j){
					update(j,1);
				}
				R=liu[i].r;
			}
			Ans[liu[i].id]=ans;
		}
		for(int i=1;i<=m;++i){
			printf("%lld\n",Ans[i] );
		}
	return 0;
}

bzoj4241

Cwolf9
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
BZOJ JOI 2013~2014 春季training合宿 系列题解1
08-04
BZOJ JOI 2013~2014 春季training合宿 系列题解1 本文主要讲述哈希算法在OJ系统中的应用，通过对BZOJ JOI 2013~2014 春季training合宿 系列题解1的分析，来阐述哈希算法的原理和实现。 首先，我们来看一下哈希算法...
「UOJ224」「NOI2016」旷野大计算
infinity_edge的博客
06-30 1187
题目描述 随着人类计算机技术的发展，计算机的能力不断提升，让跳蚤国王非常羡慕。 终于有一天,跳蚤国王发布政令：大力发展跳蚤国的计算机产业！然而，跳蚤国尚未进行工业革命，无法制造出电子计算机所需的元器件。但是跳蚤国王想出了一个绝妙的想法：把每只跳蚤作为一个计算节点，每只跳蚤只完成一个特定的小任务。 跳蚤国王带领 nnn 只跳蚤来到了一片旷野上，把跳蚤作为计算节点在旷野上排列好，并编号为 111...
一维二维前缀和差分与莫队总结
litxiaoyao的博客
08-16 297
前缀和差分与莫队总结 一.前缀和与差分 1.基础知识 a.一维 前缀和与差分运算为互逆运算，任意一个数组a的前缀和数组的差分数组是它本身。 //为了避免数组越位，下表从1开始用 for (int i = 1; i <= n; ++i) { s[i] = s[i - 1] + a[i]; } 0 1 2 3 4 5 a 2 6 8 4 3 s 0 2 8 16 20..
二维莫队复杂度分析及块长计算
最新发布
PLUS_CAT的博客
04-17 1012
设一个矩阵长N，宽M，题目询问次数为Q，块长Bn，块宽Bm，指针（矩形端点坐标）为xa​ya​xb​yb​。那么，xa​xb​Bn，单次更新时间复杂度OM​，同理，ya​Bm，单次更新时间复杂度ON。而yb​对于Xidxa​​Xidxb​​Yidya​​不变的情况下，最大更新次数N，单次更新时间复杂度ON。
【BZOJ2639】矩形计算（二维莫队）
C20200324的博客
11-15 385
题目大意： 和我一维莫队（复杂度计算在内）（刚写的，很烂，但还是想要点击量）用的例题一模一样，只不过是二维。 思路： 好像是雅礼的题，正解恶心人，分两个方法分别算出现次数多和少的数。 但是呢，我懒，还是二维莫队（真的就是把一维莫队搬上二维，两个左右下标）了。 排序是先排第一维，在相同时排第二维。 复杂度证明和一维莫队基本一致。 再在上面乘以一个n（就是改变一维下标，对另一维进行O（n）查询）， 但这道题的n只有200。 复杂度为O（n*（q*sqrt(n)+n*sqrt(n)）） 2秒应
Important matches Gym - 101840I —— 二维莫队
天翼之城的博客
11-11 435
In the qualification rounds of the World Cup, there were two disjoint sets of teams S1 = [1, · · · , N] and S2 = [1, · · · , M]. Given a matrix X that decides the importance of the match between team ...
[NOI2016]旷野大计算
df4516的博客
04-22 658
Subtask0 造计算机神题。给一个忠告：珍爱生命，远离旷野大计算...... 代码在这里：戳我 Subtask1 给定\(a,b\)；求\(-2a-2b\)。 熟悉操作环境：\([-(a+b)]<<1\)，共\(6\)次操作。 Subtask2 给定\(a\)；求\(\frac{1}{1+e^{17a}}\)。 熟悉关键函数及操作：\(s(-((a<<...
joi_221202_v7.2.0_32231.apk
12-05
joi_221202_v7.2.0_32231.apk
joi_230707_7.7.0_44416_share-2ERL05.apk
07-09
joi_230707_7.7.0_44416_share-2ERL05.apk
第19回日本情報オリンピック春合宿オンラインコンテストRound2 1日目 注销.zip
03-20
Everybody is welcome to attend JOI Open Contest! The duration of the contest is 5 hours. The same contests are held twice (Round 1 and Round 2). Contestants can participate in one of Round 1 or ...
express-joi-validation：使用joi验证Express应用程序的输入和参数
02-05
快速Joi验证 使用Joi模式验证快速输入的中间件。 功能包括： TypeScript支持。 指定验证请求输入的顺序。 将传入的req.body ， req.query等替换为经过验证的结果将原始req.body保留在名为req.originalBody的新属性中...
ACM专用模板，包括数据结构，图论，字符串，数论，几何计算，KD树，ST表，二叉搜索树，二维ST表（任意矩形，正方形），二维树状数组模板，分块ACM模板.rar
07-26
ACM专用模板，包括数据结构，图论，字符串，数论，几何计算，KD树，ST表，二叉搜索树，二维ST表（任意矩形，正方形），二维树状数组模板，分块，权值线段树 区间最早出现问题，树状数组模板，线段树模板，线段树扫描线，主席树模板
UOJ224 NOI2016 旷野大计算 构造、造计算机
weixin_30505751的博客
12-25 156
传送门——UOJ 传送门——Luogu 这段时间请不要找Itst聊天，Itst已经做疯了 事实证明大模拟题不可做 query 1 送分，加起来一起乘即可 I I + 1 2 < 3 1 - 4 O 5 query 2 送分 I - 1 < 2 4 + 2 3 S 4 O 5 query 3 送分，把原数扩大$INF$倍，再...
【jzoj4763】【旷野大计算】【莫队】
inklutcuah’s blog
09-22 998
题目大意给出长度为n的序列a，有一些询问a[l,r]中出现过的数字乘以出现的次数的最大值。解题思路可以考虑使用莫队，但是我们发现加一个数简单，删一个数就有点难了，可以使用线段树维护，然而这会tle。我们可以尽量使用加，按l所在区间降序为第一关键字，r升序为第二关键字。维护一个桶，先不考虑l所在区间的部分，考虑完整的块和右边，可以轻松维护。再考虑l所在区间的部分，先存下当前的答案，暴力加上l所在区间的
UOJ224/洛谷P1737 【NOI2016】旷野大计算 造计算机
litble的成(tui)长(fei)史
05-31 2545
论一条咸鱼是有多闲系列。 浑浑噩噩做了两天，感觉自己造计算机的水平得到了本质上的提升。 官方题解看不懂啊喂，民间题解在哪找啊喂。 test1 非常简单。 I I + 1 2 - 3 + 4 4 O 5 test2 非常简单。 I &amp;amp;lt; 1 4 + 1 2 - 3 S 4 O 5 test3 可以利用一下精度的性质。 我们可以将这个a乘以一个很大很大的数，这...
莫队算法及各种莫队扩展
niiick
05-15 457
了解莫队之前先看一下这样一个问题 Q：有一个长为N序列，有M个询问：在区间[L,R]内，出现了多少个不同的数字。(数字范围为0到1000000之间的整数)，N ≤ 50000，M ≤ 200000。 对于这个问题不难想到一个简单的暴力 int L=1,R=0; void add(int x) { cnt[x]++; if(cnt[x]==1) sum++; } void d...
【NOIP模拟】旷野大计算
Facico的博客
09-09 1756
Description**Solution这是一道套路题，不过我现在才学会套路。 对于莫队算法，进队统计答案很快，出队统计答案很慢的情况有一个套路：把所有的出队变成进队，每次询问的[l,r]，把l放到l所在的块的末尾为l’，然后把l到l’-1的都入队，更新答案，因为右端点是递增的，所以右边就逐个进队就好了，然后结束的时候就把l’到l-1退队。每次l换块的时候，就把所有的数据暴力重构。Code#in
平衡规划后的华丽暴力——莫队算法
weixin_30588907的博客
07-12 195
Perface 莫队算法是一种用于处理一类支持离线的区间查询问题（基本上都和颜色有关）的暴力算法。 它的主要思想是暴力查询，每次以尽量少的步数移动答案区间位置，从而达到一个可以接受的时间复杂度。 普通莫队 例题：给定一个长度为\(N(N\in[1,10^5])\)的数列，有\(M(M\in[1,10^5])\)个询问，每次询问区间\([a_i,b_i]\)中有多少种不同的数字。 考虑离线。我...
[bzoj2639]矩形计算——二维莫队or四维数点
ylsoi的博客
10-09 1090
题目大意： 输入一个n*m的矩阵，矩阵的每一个元素都是一个整数，然后有q个询问，每次询问一个子矩阵的权值。矩阵的权值是这样定义的，对于一个整数x，如果它在该矩阵中出现了p次，那么它给该矩阵的权值就贡献p2p^2p2。 思路： 二维莫队的模板题，但是复杂度好像不太明确。。。 于是可以换一种思路，考虑平方的几何意义，即两个颜色相同的点可以构成一个点对，每一次询问在某矩阵中出现的点对个数。 对于每一个矩...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值