P4396 [AHOI2013]作业 cdq分治

传送门

文章目录

  • 题意:
  • 思路:

题意:

给你一个长度为nnn的数列aaa,有qqq个询问,每次询问[l,r][l,r][l,r]中值域在[a,b][a,b][a,b]中的数出现的次数和在[a,b][a,b][a,b]中出现过的数值个数。

n≤1e5,1≤a≤1e5n\le1e5,1\le a \le 1e5n1e5,1a1e5

思路:

这个题可以值域分块来写,经典模型了。

这里介绍cdqcdqcdq的写法,考虑分的三维是那三维。

比较容易想到前两维是l≤pos≤r,a≤val≤bl\le pos\le r,a\le val\le blposr,avalb,第二个答案的统计是难点,考虑如何统计第二个答案。记pre[val]pre[val]pre[val]代表valvalval之前出现的位置,我们如果能知道有多少0≤pre[val]≤l−10\le pre[val]\le l-10pre[val]l1是不是就可以了,因为[a,b][a,b][a,b]中的每个值的preprepre[0,l−1][0,l-1][0,l1]最多出现一次,所以我们第三维就是这个?

并不是,不能忽略的一个很重要的问题就是cdqcdqcdq需要保证先进行修改,再进行查询,所以第三维应该保证修改在查询之前,与pre[val]pre[val]pre[val]是多少无关,因为第一第二维已经定序。

这里给一组hackhackhack样例

2 1

2 2

1 2 2 2

当然这是对于不离散化来说的,如果你离散化了,那么就另当别论了。

复杂度O(nlog2n)O(nlog^2n)O(nlog2n)

// Problem: P4396 [AHOI2013]作业
// Contest: Luogu
// URL: https://www.luogu.com.cn/problem/P4396
// Memory Limit: 125 MB
// Time Limit: 1000 ms
// 
// Powered by CP Editor (https://cpeditor.org)//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#include<random>
#include<cassert>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid ((tr[u].l+tr[u].r)>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
#define lowbit(x) ((x)&(-x))
using namespace std;//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;const int N=6000010,mod=1e9+7,INF=0x3f3f3f3f;
const double eps=1e-6;int n,m;
int a[N],pre[N],tot;
struct Node {int x,y,z,op,id,val;
}q[N],p[N];
vector<int>v;bool cmp(Node a,Node b) {if(a.x!=b.x) return a.x<b.x;else if(a.y!=b.y) return a.y<b.y;else return a.op<b.op;}int tr[N],ans[N],ans1[N];void add(int x,int c) {for(int i=x;i<N;i+=lowbit(i)) tr[i]+=c;
}int sum(int x) {int ans=0;for(int i=x;i;i-=lowbit(i)) ans+=tr[i];return ans;
}void cdq(int l,int r) {if(l>=r) return;int mid=(l+r)>>1;cdq(l,mid); cdq(mid+1,r);int i=l,j=mid+1;int now=0,cnt=0;while(i<=mid&&j<=r) {// if(l==1&&r==3&&q[i].y==q[j].y) {// cout<<q[i].op<<' '<<q[j].op<<"**"<<endl;// } if(q[i].y<=q[j].y) {p[++now]=q[i];if(q[i].op==0) {add(q[i].z+1,1);cnt++;}i++;} else {p[++now]=q[j];if(q[j].op==1) {ans[q[j].id]+=sum(q[j].z+1)*q[j].val;ans1[q[j].id]+=cnt*q[j].val;}j++;}}while(i<=mid) {p[++now]=q[i];if(q[i].op==0) {add(q[i].z+1,1);cnt++;}i++;}while(j<=r) {p[++now]=q[j];if(q[j].op==1) {ans[q[j].id]+=sum(q[j].z+1)*q[j].val;ans1[q[j].id]+=cnt*q[j].val;}j++;}for(int i=l;i<=mid;i++) if(q[i].op==0) add(q[i].z+1,-1);for(int i=l,j=1;i<=r;i++,j++) q[i]=p[j];
}int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);scanf("%d%d",&n,&m);for(int i=1;i<=n;i++) {int x; scanf("%d",&x);a[i]=x; v.pb(x);}v.pb(0);sort(v.begin(),v.end());for(int i=1;i<=n;i++) {int x; x=a[i];q[++tot]={i,x,pre[x],0,0,0};pre[x]=i;}for(int i=1;i<=m;i++) {int l,r,a,b; scanf("%d%d%d%d",&l,&r,&a,&b);q[++tot]={r,b,l-1,1,i,1};q[++tot]={l-1,a-1,l-1,1,i,1};q[++tot]={l-1,b,l-1,1,i,-1};q[++tot]={r,a-1,l-1,1,i,-1};}sort(q+1,q+1+tot,cmp);cdq(1,tot);for(int i=1;i<=m;i++) printf("%d %d\n",ans1[i],ans[i]);return 0;
}
/*
l<=x<=r
a<=val<=b
0<=pos<=l-1x y z op val0 1 0 1 1
0 2 0 1 -1
1 2 0 0 02 1 0 1 -1
2 2 0 1 1
2 2 1 0 00 1 0 1 1
0 2 0 1 -1
1 2 0 0 0
2 1 0 1 -1
2 2 0 1 1
2 2 1 0 0*/

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/314593.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

WebClient, HttpClient, HttpWebRequest ,RestSharp之间的区别与抉择

NETCore提供了三种不同类型用于生产的REST API&#xff1a; HttpWebRequest;WebClient;HttpClient&#xff0c;开源社区创建了另一个名为RestSharp的库。如此多的http库&#xff0c;该怎样选择呢&#xff1f;01HttpWebRequest这是.NET创建者最初开发用于使用HTTP请求的标准类。…

树的直径与重心

树的直径与重心 或许更好的阅读体验 树的直径求解方法一 思路 先选取一个点rt作为根节点&#xff0c;dfs去找到一个最长路径的点U&#xff0c;然后通过这个点去dfs&#xff0c;找到路径最长的点V&#xff0c;U->V就是这课树的直径。 证明正确性&#xff1a; 假如rt在直…

Codeforces Round #446 (Div. 1) B. Gluttony 构造 + 补集思想

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个数组aaa&#xff0c;保证aaa中每个数都互不相同&#xff0c;让你构造一个数组bbb&#xff0c;满足对于任意的Sx1,x2,...,xk,1≤xi≤n,0≤k<nS{x_1,x_2,...,x_k},1\le x_i\le n,0\le k< nSx1​…

.NetCore 2.1以后的HttpClient最佳实践

ASP.NET Core 2.1中出现一个新的HttpClientFactory功能&#xff0c;它有助于解决开发人员在使用HttpClient实例从其应用程序发出外部Web请求时可能遇到的一些常见问题。01介绍在.NETCore平台的2.1新增了HttpClientFactory&#xff0c;虽然HttpClient这个类实现了disposable&…

D. Multiset(树状数组 + 二分)

Multiset 可能更好的阅读体验 思路 二分 树状数组做法 我们发现每个数的范围是$ < 1e6$的&#xff0c;所以可以直接在线操作&#xff0c;不用离散化离线操作。 这个时候我们的treetreetree数组就相当与一个桶&#xff0c;每个桶里统计的是值为其下标的个数&#xff0c…

Codeforces Round #381 (Div. 1) A. Alyona and mex 构造

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 你需要确定一个长度为nnn的数组aaa&#xff0c;满足给定的mmm个[l,r][l,r][l,r]限制&#xff0c;需要保证构造出来的aaa数组对于每个[l,r][l,r][l,r]的mexmexmex最小值最大。 n,m≤1e5n,m\le1e5n,m≤1e5 思…

点分治(简要讲解 + 模板)

树上点分治 思想 两个点之间的距离无非就是两种关系&#xff1a;我们约定dis[i]dis[i]dis[i]表示这个点到当前根节点的距离 dis[u]dis[v]dis[u] dis[v]dis[u]dis[v]&#xff0c;在同一个根节点的不同子树上。dis[u]dis[v]dis[u] dis[v]dis[u]dis[v]&#xff0c;在同一个棵子…

在 alpine 中使用 NPOI

在 alpine 中使用 NPOIIntro在 .net 中常使用 NPOI 来做 Excel 的导入导出&#xff0c;NPOI 从 2.4.0 版本开始支持 .netstandard2.0&#xff0c;在 dotnet core 应用也可以用 DotNetCore.NPOI。对于 .NET Core 应用来说&#xff0c;如果没有特殊的需求&#xff0c;alpine 是最…

Educational Codeforces Round 16 C. Magic Odd Square 矩阵构造

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个奇数nnn&#xff0c;让你构造一个n∗nn*nn∗n的矩阵&#xff0c;矩阵的每个位置依次填上[1,n∗n]之内的数[1,n*n]之内的数[1,n∗n]之内的数&#xff0c;满足每行、每列、以及主对角线的和都是奇数。…

关于WCF、WebAPI、WebService之间的区别总结 分布式通信技术

早在1996年Gartner就前瞻性地提出了面向服务架构的思想(SOA)&#xff0c;SOA 的走红在很大程度上归功于 Web Service 标准的成熟和应用的普及。Service Oriented Ambiguity 中文一般理解为&#xff1a;面向服务架构&#xff0c;简称SOA&#xff0c;这个概念算得上微服务的鼻祖了…

LCA求解的四种模板

LCA求解的四种模板 或许更好的阅读体验 树剖在线求解LCA 思想 树剖这里就不多解释了&#xff0c;求解LCA的过程就是轻重链的跳转&#xff0c;跟树剖求任意两点间的距离一样的操作&#xff0c;只不过不用线段树去维护disdisdis了&#xff0c;那就直接上代码吧。 代码 #inc…

HDU - 1998 奇数阶魔方

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个奇数nnn&#xff0c;构造一个nnn阶幻方。 3≤n≤193\le n\le 193≤n≤19 思路&#xff1a; 模板题了&#xff0c;直接构造一个幻方即可。 首先在第一行中间放一个111&#xff0c;之后每次看一下(…

Azure 物联网开发者体验 7 月更新:边缘计算开发工具,ARM64 设备开发,VS Code 容器化开发工具...

欢迎来到 Azure 物联网开发者体验的 7 月更新&#xff01;在本次发布中&#xff0c;微软为物联网开发人员提供了许多新的功能和改进&#xff01;物联网边缘计算开发工具正式发布Azure IoT Edge 于 2017 年发布。随着近两年的对边缘计算开发工具的不断开发与改善&#xff0c;Azu…

Tarjan缩点

Tarjan缩点 或许更好的阅读体验 P3387 【模板】缩点 思路 既然时缩点的模板&#xff0c;那么缩点自然少不了了&#xff0c;缩点后我们的到新的有向无环图&#xff0c;然后再利用这个无环图去找一条最大权值的路径&#xff0c;路径和即为答案。 我们改如何选取起点来避免不…

.NET中的值类型与引用类型

.NET中的值类型与引用类型这是一个常见面试题&#xff0c;值类型(Value Type)和引用类型(Reference Type)有什么区别&#xff1f;他们性能方面有什么区别&#xff1f;TL;DR&#xff08;先看结论&#xff09;值类型引用类型创建位置栈托管堆赋值时复制值复制引用动态内存分配无需…

Canada Cup 2016 C. Hidden Word 字符串构造

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个长度为272727的字符串sss&#xff0c;其中262626个字母每种都至少出现一次&#xff0c;让你构造一个213213213的矩阵&#xff0c;使得每个字母都出现一次&#xff0c;并且存在一条路径&#xff0c;经…

E:Johnny and Grandmaster

Johnny and Grandmaster 或许更好的阅读体验 思路 这道题就是把一组数分成两个集合&#xff0c;使这两个集合的对p的次方的和的差的最小值&#xff0c;也就是求sum1−sum2sum1 - sum2sum1−sum2得最小值&#xff0c; 由于结果过大我们可能需要对结果取模。那么这题得关键在于…

跨语言调用Hangfire定时作业服务

背景Hangfire允许您以非常简单但可靠的方式执行后台定时任务的工作。内置对任务的可视化操作。非常方便。但令人遗憾的是普遍都是业务代码和hagnfire服务本身聚合在一个程序中运行&#xff0c;极大的限制了hangfire的扩展和跨语言调用。所以萌生了开发一个支持restful api调用的…

Codeforces Round #453 (Div. 1) D. Weighting a Tree 构造 + dfs树

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一颗nnn个点的图&#xff0c;每个点都有一个点权cic_ici​&#xff0c;要求你给每个边赋一个权值kik_iki​&#xff0c;要求对于每个点与他相连的边的权值之和等于这个点的点权cic_ici​。 n≤1e5,n−1≤…

Codeforces Round #648 (Div. 2)(A, B, C, D)

Codeforces Round #648 (Div. 2) 或许更好的阅读体验 A:Matrix Game 思路 题意可以说是非常简单的&#xff0c;我们选定的格子的行列都不能存在1&#xff0c;可以发现我们可以放的格子一定是固定的&#xff0c;然后这题就变成了技术总共可以放多少个棋子了&#xff0c;所以…