luogu P4062 [Code+#1]Yazid 的新生舞会(线段树+套路)
今天原来是平安夜啊
感觉这题是道好题。
一个套路枚举权值\(x\),把权值等于\(x\)的设为1,不等于的设为-1,然后问题转化为多少个区间权值和大于。
发现并不是很好做,还有一个套路,用前缀和查分来表示区间。然后就是
\]
\]
然后树状数组可以做\(a[i]\leq7\)的数据。
那么\(a[i]\)那么大该怎么办?
考虑我们构建的\(1,-1\)数列中连续-1的数列很多。
然后这些连续-1不会互相影响的贡献,然后我们考虑直接算出这些连续-1的贡献。
设这段连续-1的两端的sum为\(l,r\)。贡献为:
\]
然后我们在线段树上维护\(i\)和\(cnt[i]\),然后用线段树可以把这个-1序列\(logn\)处理掉。
然后每一个1单独处理,最后总复杂度\(O(nlogn)\)
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
using namespace std;
#define int long long
#define mid ((l+r)>>1)
#define ls now<<1
#define rs now<<1|1
const int N=501000;
vector<int> vec[N];
int sum[N*2*5],sumi[N*2*5],lazy[N*2*5],a[N],b[N],last,lastsum,ans,n;
int read(){
int sum=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){sum=sum*10+ch-'0';ch=getchar();}
return sum*f;
}
void update(int now){
sum[now]=sum[ls]+sum[rs];
sumi[now]=sumi[ls]+sumi[rs];
}
void pushdown(int l,int r,int now){
if(lazy[now]==0)return;
if(l==r)return;
sum[ls]+=(mid-l+1)*lazy[now];
sumi[ls]+=(mid+l)*(mid-l+1)/2*lazy[now];
sum[rs]+=(r-(mid+1)+1)*lazy[now];
sumi[rs]+=(r+mid+1)*(r-(mid+1)+1)/2*lazy[now];
lazy[ls]+=lazy[now];
lazy[rs]+=lazy[now];
lazy[now]=0;
}
void add(int l,int r,int L,int R,int k,int now){
if(R<l||L>r)return;
pushdown(l,r,now);
if(L<=l&&r<=R){
sum[now]+=(r-l+1)*k;
sumi[now]+=(l+r)*(r-l+1)/2*k;
lazy[now]+=k;
return;
}
add(l,mid,L,R,k,ls);
add(mid+1,r,L,R,k,rs);
update(now);
}
int check(int l,int r,int L,int R,int now){
if(R<l||L>r)return 0;
pushdown(l,r,now);
if(L<=l&&r<=R)return sum[now];
return check(l,mid,L,R,ls)+check(mid+1,r,L,R,rs);
}
int checki(int l,int r,int L,int R,int now){
if(R<l||L>r)return 0;
pushdown(l,r,now);
if(L<=l&&r<=R)return sumi[now];
return checki(l,mid,L,R,ls)+checki(mid+1,r,L,R,rs);
}
signed main(){
n=read();int hh=read();
for(int i=1;i<=n;i++)a[i]=read(),b[i]=a[i];
sort(b+1,b+1+n);
int tot=unique(b+1,b+1+n)-b-1;
for(int i=1;i<=n;i++)a[i]=lower_bound(b+1,b+1+tot,a[i])-b;
for(int i=1;i<=n;i++)vec[a[i]].push_back(i);
add(1,2*n,n,n,1,1);
for(int i=1;i<=tot;i++){
vec[i].push_back(n+1);
last=0;lastsum=0;
for(int j=0;j<vec[i].size();j++){
int L=vec[i][j]-last-1;
int l=lastsum-L;
int r=lastsum-1;
if(vec[i][j]!=last+1){
ans+=L*check(1,2*n,1,l+n-1,1);
ans+=(r+n)*check(1,2*n,l+n,r+n,1);
ans-=checki(1,2*n,l+n,r+n,1);
add(1,2*n,l+n,r+n,1,1);
}
if(vec[i][j]==n+1)continue;
ans+=check(1,2*n,1,l+n,1);
add(1,2*n,l+1+n,l+1+n,1,1);
last=vec[i][j];lastsum=l+1;
}
last=0;lastsum=0;
for(int j=0;j<vec[i].size();j++){
int L=vec[i][j]-last-1;
int l=lastsum-L;
int r=lastsum-1;
if(vec[i][j]!=last+1)add(1,2*n,l+n,r+n,-1,1);
if(vec[i][j]==n+1)continue;
add(1,2*n,l+1+n,l+1+n,-1,1);
last=vec[i][j];lastsum=l+1;
}
}
printf("%lld\n",ans);
return 0;
}
luogu P4062 [Code+#1]Yazid 的新生舞会(线段树+套路)的更多相关文章
- 【BZOJ5110】[CodePlus2017]Yazid 的新生舞会 线段树
[BZOJ5110][CodePlus2017]Yazid 的新生舞会 Description Yazid有一个长度为n的序列A,下标从1至n.显然地,这个序列共有n(n+1)/2个子区间.对于任意一 ...
- P4062 [Code+#1]Yazid 的新生舞会
思路:分治 提交:2次 错因:数组开小 题解: 我们枚举一下众数\(x\). 设\(s[n]=\sum_{i=1}^n [a[i]==x]\) 那么对于区间\((l,r]\),有\(s[r]-s[l] ...
- 洛谷 P4062 - [Code+#1]Yazid 的新生舞会 的线性做法
洛谷题面传送门 一个线性做法. \(n\log n\) 解法可以戳这里查看 首先回顾一下 \(n\log n\) 解法的过程:我们对于每一个数 \(x\),考察其出现位置,设为 \(t_1,t_2,t ...
- 洛谷 P4062 - [Code+#1]Yazid 的新生舞会(权值线段树)
题面传送门 题意: 给出一个序列 \(a\),求 \(a\) 有多少个子区间 \([l,r]\),满足这个区间中出现次数最多的数出现次数 \(>\dfrac{r-l+1}{2}\) \(1 \l ...
- BZOJ.5110.[CodePlus2017]Yazid 的新生舞会(线段树/树状数组/分治)
LOJ BZOJ 洛谷 又来发良心题解啦 \(Description\) 给定一个序列\(A_i\).求有多少个子区间,满足该区间众数出现次数大于区间长度的一半. \(n\leq5\times10^5 ...
- 【线段树】【P4062】 [Code+#1]Yazid 的新生舞会
Description 给定一个长度为 \(n\) 的序列,求有多少子区间满足区间众数严格大于区间长度的一半.如果区间有多个出现次数最多且不同的数则取较小的数为众数. Limitation 对于全部的 ...
- [题解] [Code+#1]Yazid 的新生舞会
题面 题解 upd : \(cnt_i\) 代表值为 \(i\) 的个数 我们可以暴力枚举众数 \(k\) 把等于 \(k\) 的赋值成 1 , 不等于 \(k\) 的赋值成 -1 这样原序列就变成了 ...
- 「CodePlus 2017 11 月赛」Yazid 的新生舞会(树状数组/线段树)
学习了新姿势..(一直看不懂大爷的代码卡了好久T T 首先数字范围那么小可以考虑枚举众数来计算答案,设当前枚举到$x$,$s_i$为前$i$个数中$x$的出现次数,则满足$2*s_r-r > 2 ...
- luogu P3919 [模板]可持久化数组(可持久化线段树/平衡树)(主席树)
luogu P3919 [模板]可持久化数组(可持久化线段树/平衡树) 题目 #include<iostream> #include<cstdlib> #include< ...
随机推荐
- 蓝桥杯_基础训练_Sine之舞
基础练习 Sine之舞 时间限制:1.0s 内存限制:512.0MB 问题描述 最近FJ为他的奶牛们开设了数学分析课,FJ知道若要学好这门课,必须有一个好的三角函数基本功.所以他准备和奶牛们做 ...
- 【XSY2692】杨柳 - 网络流
题目来源:2018冬令营模拟测试赛(十) 题解: 继续鬼畜网络流…… 首先这题有个显然的做法:bfs预处理出每个起点到每个终点的最短步数,然后直接建边加超级源汇跑费用流即可: 但是这样边数是$n^2$ ...
- IOS - 退出程序
- (void)exitApplication { OAAppDelegate *app = [UIApplication sharedApplication].delegate; UIWindow ...
- JS中的异步
Hello,日常更新的我“浪”回来了!!! JS中有三座高山:异步和单线程.作用域和闭包.原型原型链 今天“浪”的主题是JS中的异步和单线程的问题. 主要从这三个方面入手 一.什么是异步(与同步作比较 ...
- iOS 开发仿网易云音乐歌词海报
使用网易云音乐也是一个巧合,我之前一直使用QQ音乐听歌,前几天下 app 手机内存告急.于是就把QQ音乐给卸载掉了,正好晚上朋友圈里有一个朋友用网易云音乐分享了一首歌曲,于是我也就尝试下载了网易云音乐 ...
- 【BZOJ3270】博物馆 概率DP 高斯消元
链接: #include <stdio.h> int main() { puts("转载请注明出处[辗转山河弋流歌 by 空灰冰魂]谢谢"); puts("网 ...
- MonkeyRunner源代码分析之启动
在工作中由于要追求完毕目标的效率,所以很多其它是强调实战.注重招式.关注怎么去用各种框架来实现目的.可是假设一味仅仅是注重招式.缺少对原理这个内功的了解,相信自己非常难对各种框架有更深入的理解. 从几 ...
- VS 2013+Qt 5.4.1
Qt应用能够用Qt Creator开发,也能够使用Visual Studio. 我之前一直用Qt Creator.也始终认为这是最好的选择.只是有人偏爱Visual Studio,我也由于工作须要,要 ...
- 英语影视台词---七、THE GREAT GATSBY QUOTES
英语影视台词---七.THE GREAT GATSBY QUOTES 一.总结 一句话总结:了不起的盖茨比 1.“So we beat on, boats against the current, b ...
- [LeetCode] Longest Substring Without Repeating Characters 最长无重复字符的子串 C++实现java实现
最长无重复字符的子串 Given a string, find the length of the longest substring without repeating characters. Ex ...