P4062 [Code+#1]Yazid 的新生舞会 树状数组维护三阶差分

传送门

文章目录

  • 题意:
  • 思路:

题意:

给你一个序列aaa,让你求有多少个子区间满足存在一个数是这个区间的绝对众数,绝对众数指该数在区间内出现的次数严格大于r−l+12\frac{r-l+1}{2}2rl+1
n≤5e5,0≤ai≤n−1n\le5e5,0\le a_i\le n-1n5e5,0ain1

思路:

考虑将每个数单独拿出来考虑贡献,将选出来的数原来位置置为111,其余位置置为−1-11,设sumsumsum为修改后数组的前缀和,那么如果[l,r][l,r][l,r]满足要求的话需要sumr−suml−1>r−l+1−(sumr−suml−1)sum_r-sum_{l-1}>r-l+1-(sum_r-sum_{l-1})sumrsuml1>rl+1(sumrsuml1),令l=l+1l=l+1l=l+1并且移项得2∗sumr−r>2∗suml−l2*sum_r-r>2*sum_l-l2sumrr>2sumll,所以我们就可以维护一个树状数组,每次求一下[−n,2∗sumr−r)[-n,2*sum_r-r)[n,2sumrr)有多少数即可。
但是字符集很大,与nnn同阶,直接计算是n2lognn^2lognn2logn的。
我们某个字符发现有xxx个的时候,可以分为x+1x+1x+1段,而且段内都是公差为−1-11的等差数列,单调递减,所以段内数之间没有贡献,我们考虑将整个区间当成整体来看,这样整体复杂度就从n2−>nn^2->nn2>n了。
怎么将整个区间看成整体呢?
考虑值有可能为负数,我们加一个偏移量即可。设cic_ici是值为iii的个数,那么每次查询完区间[x,y][x,y][x,y],那么给x,yx,yx,y区间对应的权值加111即可,设fi=∑j=1icjf_i=\sum_{j=1}^ic_jfi=j=1icj,代表权值的前缀和,那么[x,y][x,y][x,y]的答案就是∑i=x−1y−1fi\sum_{i=x-1}^{y-1}f_ii=x1y1fi,又是一个前缀和的形式,我们考虑设gi=∑j=1ifjg_i=\sum_{j=1}^{i}f_jgi=j=1ifj,那么答案就是gy−1−gx−2g_{y-1}-g_{x-2}gy1gx2
到目前位置我们就可以用线段树实现这个二阶前缀和即可,但是树状数组更快,所以考虑将[x,y][x,y][x,y]区间操作也打个差分,所以变成了三阶前缀和,推导如下(侵删)
在这里插入图片描述

所以维护一下di,di∗i,di∗i∗id_i,d_i*i,d_i*i*idi,dii,diii即可。

// Problem: P4062 [Code+#1]Yazid 的新生舞会
// Contest: Luogu
// URL: https://www.luogu.com.cn/problem/P4062
// Memory Limit: 500 MB
// Time Limit: 4000 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//in.in","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=4000010,mod=1e9+7,INF=0x3f3f3f3f;
const double eps=1e-6;int n,type;
LL tr1[N],tr2[N],tr3[N];
int a[N];
vector<int>pos[N];void add(int x,int c) {for(int i=x;i<=n*2+2;i+=lowbit(i)) {tr1[i]+=c;tr2[i]+=1ll*c*x;tr3[i]+=1ll*c*x*x;}
}LL sum(int x) {if(x<=0) return 0;LL ans=0;for(int i=x;i;i-=lowbit(i)) ans+=tr1[i]*(x+2)*(x+1)-tr2[i]*(2*x+3)+tr3[i];return ans/2;
}int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);// rd_ac();int _; scanf("%d",&_);while(_--) {scanf("%d",&n);for(int i=1;i<=n;i++) scanf("%d",&a[i]),pos[a[i]].pb(i);int rc=n+1; LL ans=0;for(auto t:pos) {t.pb(n+1);int pre=0,cnt=0;for(auto now:t) {int x=2*cnt-(now-1)+rc,y=2*cnt-pre+rc;ans+=sum(y-1)-sum(x-2);add(x,1); add(y+1,-1);cnt++; pre=now;}pre=0; cnt=0;for(auto now:t) {int x=2*cnt-(now-1)+rc,y=2*cnt-pre+rc;add(x,-1); add(y+1,1);cnt++; pre=now;}}printf("%lld\n",ans);for(int i=1;i<=n;i++) pos[a[i]].clear();}return 0;
}
/**/

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

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

相关文章

一句话概括4本管理著作

之前618活动的时候买了不少书&#xff0c;最近一段时间集中看了四本管理相关的书籍&#xff0c;收获颇丰。在这里分享给大家。如果你不是管理者也没关系&#xff0c;可以换个视角来看看管理者眼中的世界&#xff0c;毕竟&#xff0c;“赤兔马”也需要让“伯乐”看到才行嘛。虽然…

[译].Net中的内存-什么分配在了哪里

原文链接&#xff1a;https://jonskeet.uk/csharp/memory.html人们在理解值类型和引用类型之间的差异时因为“值类型在栈上分配&#xff0c;引用类型在堆上分配”这句话造成了很多混乱。这完全是不对的&#xff0c;本文试图澄清这个问题。变量中有什么&#xff1f;理解.NET中内…

上车时机已到--.NETCore是适应时代发展的雄鹰利剑

要起飞了.NET Core 3.0-prevew7&#xff1a;https://dotnet.microsoft.com/download/dotnet-core/3.0随着 .NET Core 3.0-prevew7 的发布&#xff0c;开源社区的一支重要力量重要即将起飞&#xff1b;官方指出&#xff0c;该预览版为可能为最终版本&#xff0c;在很长一段时间内…

HDU - 7028 Decomposition 无向完全图构造欧拉回路

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一张无向完全图&#xff0c;让你构造kkk个长度分别为lil_ili​的路径&#xff0c;这些路径不相交&#xff0c;且∑lin∗(n−1)2\sum l_i\frac{n*(n-1)}{2}∑li​2n∗(n−1)​。 n≤1000,nmod21n\le1000,…

「Azure」数据分析师有理由爱Azure之二-立即申请帐号开始学习之旅

目前关于Azure的学习资料不多&#xff0c;除了官方的文档和Microsoft Learn频道外&#xff0c;几乎没有什么中文性资料可学习&#xff0c;就算有&#xff0c;也是以IT的思维方式来展开介绍&#xff0c;对没有IT背景的数据分析师来说&#xff0c;非常难于适应。本篇给大家带来一…

矩阵快速幂的最简单用法

矩阵快速幂 链接&#xff1a;https://ac.nowcoder.com/acm/contest/1168/K 来源&#xff1a;牛客网 题目描述 这个勇者明明超强却过分慎重&#xff0c;勇者龙宫院圣哉与n名冒险者一起去讨伐神秘魔物&#xff0c;龙宫院圣哉十分谨慎&#xff0c;他只会在最后一刻出手&#xff…

DevOps案例研究|史上最能“拜客户教”的公司,是如何做到持续交付的?(第1趴)...

内容来源&#xff1a;DevOps案例深度研究 –Amazon持续交付之道战队&#xff08;本文只展示部分PPT及研究成果&#xff0c;更多细节请关注案例分享会&#xff0c;及本公众号。&#xff09;本案例内容贡献者&#xff1a;单冰 (Topic Leader)、 赵栋、梁兴龙、李杰、毛艳清、牛恒…

P2261 [CQOI2007]余数求和 整除分块

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; n,k≤1e9n,k\le1e9n,k≤1e9 思路&#xff1a; 考虑转换式子&#xff0c;∑i1nkmodi∑i1n(k−⌊ki⌋∗i)n∗k−∑i1n⌊ki⌋∗i\sum_{i1}^{n}k\bmod i\sum_{i1}^n(k-\left \lfloor \frac{k}{i} \right \rfloo…

背包问题——第一篇

一&#xff0c;01背包 最简单也是最经典的背包问题。 首先我们知道背包问题是一种d问题&#xff0c;最重要的就是要去找到他的状态转移方程。而在01背包中转移方程就比较简单了&#xff0c;这里用一个二维数组进行标表示。 ans[i][j]max(ans[i-1][j],ans[i-1][j-v[i]w[i]); 在…

小白开学Asp.Net Core 《十》

小白开学Asp.Net Core 《十》 — — Session、Cookie、Cache&#xff08;老生常谈&#xff09;一、背景在常谈Session和Cookie之前我们先来简单的了解下Http&#xff08;可以说这是必须的&#xff0c;默认大家都了解&#xff09;总结一句话&#xff1a;HTTP是一种无状态的协议&…

P3935 Calculating 整除分块

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 我们设s(x)∑i1nf(x)s(x)\sum_{i1}^nf(x)s(x)∑i1n​f(x)&#xff0c;那么答案就是s(r)−s(l−1)s(r)-s(l-1)s(r)−s(l−1)。 容易发现&#xff0c;我们要求的f(x)f(x)f(x)实际上就是xxx的…

Sticks-hdu-1455深度搜索dfs

Sticks 题目传送门。 题目大意是&#xff0c;给你一个数字n代表总共的棍子数量&#xff0c;要做的就是&#xff0c;把这几根棍子拼成长度相同的棍子&#xff0c;并且让所拼成的棍子的长度尽可能地小&#xff0c;也就是拼成的棍子的数量尽可能的多。 在这里很简单的想到要给棍…

gRPC in ASP.NET Core 3.0 -- Protocol Buffer(1)

开发环境&#xff1a;IDE: VSCodeVSCode的扩展插件&#xff1a;vscode-proto3和Clang-Format这两个扩展Windows还需要安装Clang&#xff0c;Windows 64位系统的地址如下&#xff1a;Clang for Windows (64-bit)(http://llvm.org/releases/3.7.0/LLVM-3.7.0-win64.exe)&#xff…

[蓝桥杯2016初赛]方格填数

正如题目&#xff0c;我们要保证相邻的数不相连就行了&#xff0c;一共要填十个数。 下面是我的思路&#xff0c;供大家参考一下&#xff0c; 我把上面的各子放进了这个5*6的矩形&#xff0c;其中标记1–10的分别是要填数的各格子&#xff0c;一开始我先初始化这个二维数组为全…

2021牛客暑期多校训练营7 xay loves monotonicity 线段树区间合并

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 题面挺绕口的&#xff0c;还是看原题比较好。 大概的意思就是让你从给定的区间中选择一个以左端点为起点的一个上升子序列&#xff0c;让后将这些下标存下来&#xff0c;在bbb中将这些位置拿出来后&#xff…

容器化单页面应用中Nginx反向代理与Kubernetes部署

在《容器化单页面应用中RESTful API的访问》一文中&#xff0c;我介绍了一个在容器化环境中单页面应用访问后端服务的完整案例。这里我将继续使用这个案例&#xff0c;介绍一下容器化单页面应用部署的另一个场景&#xff1a;将Nginx的职责独立出来。注&#xff1a;这里单页面应…

[蓝桥杯2019初赛]迷宫-DFS、BFS两种方法

迷宫问题的最短路&#xff0c;加最小字典序 迷宫文件maze.txt传送门 作者写的2019年B组蓝桥杯解集 . . . DFS的版本 #include<iostream> #include<cstring> using namespace std; const int ax[4]{0,0,1,-1}; const int ay[4]{1,-1,0,0}; const char dir[5]{R,L…

2021牛客暑期多校训练营7 xay loves trees dfs序 + 主席树

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你两棵树&#xff0c;让你在第一棵树选一条最长的深度递增的链&#xff0c;链上每一个点在第二棵树上都不互为祖先。 n≤3e5n\le3e5n≤3e5 思路&#xff1a; 之前做过差不多的题传送门。 之前哪个题是不…

xxl-job dotnet core executor执行器开源

DotXxlJob[(github)https://github.com/xuanye/DotXxlJob][https://github.com/xuanye/DotXxlJob] xxl-job的dotnet core 执行器实现&#xff0c;支持XXL-JOB 2.01 XXL-JOB概述[XXL-JOB][1]是一个轻量级分布式任务调度平台&#xff0c;其核心设计目标是开发迅速、学习简单、轻量…

2019年第十届蓝桥杯 C / C ++省赛 B 组真题题解

A: 组队 输入数据 编号 号位1 2 3 4 51 97 90 0 0 0 2 92 85 96 0 0 3 0 0 0 0 93 4 0 0 0 80 86 5 89 83 97 0 0 6 82 86 0 0 0 7 0 0 0 87 90 8 0 97 96 0 0 9 0 0 89 0 0 10 95 99 0 0 0 11 0 0 96 97 0 12 0 0 0 93 98 13 94 91 0 0 0 14 0 83 87 0 0 15 …