luogu

这题有坨区间最大值,考虑最值分治.分治时每次取出最大值,然后考虑统计跨过这个位置的区间答案,然后两边递归处理.如果之枚举左端点,因为最大值确定,右端点权值要满足\(a_r\le \frac{\max a_k}{a_l}\),所以可以在主席树上询问区间内在一段值域内的数个数.不过如果左半边点数过多会被卡成暴力,这时枚举较小右半边即可.复杂度类似启发式合并,枚举的总点数为为\(O(nlogn)\),所以总复杂度为两个\(log\)

// luogu-judger-enable-o2
#include<bits/stdc++.h>
#define LL long long
#define uLL unsigned long long
#define db double using namespace std;
const int N=1e5+10,inf=1<<30;
int rd()
{
int x=0,w=1;char ch=0;
while(ch<'0'||ch>'9'){if(ch=='-') w=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}
return x*w;
}
int n,m,a[N],b[N],sq[N][2],bk[N];
LL ans;
#define mid ((l+r)>>1)
int s[N*30],ch[N*30][2],rt[N],tt;
void inst(int o1,int o2,int x)
{
int l=1,r=m;
s[o1]=s[o2]+1;
while(l<r)
{
if(x<=mid)
{
ch[o1][0]=++tt,ch[o1][1]=ch[o2][1];
o1=ch[o1][0],o2=ch[o2][0];
r=mid;
}
else
{
ch[o1][0]=ch[o2][0],ch[o1][1]=++tt;
o1=ch[o1][1],o2=ch[o2][1];
l=mid+1;
}
s[o1]=s[o2]+1;
}
}
int quer(int o1,int o2,int l,int r,int ll,int rr)
{
if(ll>rr||!s[o1]) return 0;
if(ll<=l&&r<=rr) return s[o1]-s[o2];
int an=0;
if(ll<=mid) an+=quer(ch[o1][0],ch[o2][0],l,mid,ll,rr);
if(rr>mid) an+=quer(ch[o1][1],ch[o2][1],mid+1,r,ll,rr);
return an;
}
int ma[N<<2];
bool zz=1;
int maxx(int aa,int bb){return a[aa]>=a[bb]?aa:bb;}
void p2(int o){ma[o]=maxx(ma[o<<1],ma[o<<1|1]);}
void bui(int o,int l,int r)
{
if(l==r){ma[o]=l;return;}
bui(o<<1,l,mid),bui(o<<1|1,mid+1,r);
p2(o);
}
int q2(int o,int l,int r,int ll,int rr)
{
if(ll<=l&&r<=rr) return ma[o];
int aa=0,bb=0;
if(ll<=mid) aa=q2(o<<1,l,mid,ll,rr);
if(rr>mid) bb=q2(o<<1|1,mid+1,r,ll,rr);
return maxx(aa,bb);
}
void sv(int l,int r)
{
if(l>r) return;
if(l==r){ans+=b[a[l]]==1;return;}
int md=q2(1,1,n,l,r);
if(md-l+1<r-md+1)
{
for(int i=md;i>=l;--i)
ans+=quer(rt[r],rt[md-1],1,m,1,upper_bound(b+1,b+m+1,b[a[md]]/b[a[i]])-b-1);
}
else
{
for(int i=md;i<=r;++i)
ans+=quer(rt[md],rt[l-1],1,m,1,upper_bound(b+1,b+m+1,b[a[md]]/b[a[i]])-b-1);
}
sv(l,md-1),sv(md+1,r);
} int main()
{
n=rd();
for(int i=1;i<=n;++i) a[i]=b[i]=rd();
sort(b+1,b+n+1),m=unique(b+1,b+n+1)-b-1;
b[++m]=inf;
for(int i=1;i<=n;++i) a[i]=lower_bound(b+1,b+m+1,a[i])-b;
for(int i=1;i<=n;++i) inst(rt[i]=++tt,rt[i-1],a[i]);
bui(1,1,n);
sv(1,n);
printf("%lld\n",ans);
return 0;
}

luogu P4755 Beautiful Pair的更多相关文章

  1. 【题解】P4755 Beautiful Pair(启发式合并的思路+分治=启发式分治)

    [题解]P4755 Beautiful Pair upd: 之前一个first second烦了,现在AC了 由于之前是直接抄std写的,所以没有什么心得体会,今天自己写写发现 不知道为啥\(90\) ...

  2. 洛谷 P4755 - Beautiful Pair(主席树+分治+启发式优化)

    题面传送门 wssb,我紫菜 看到这类与最大值统计有关的问题可以很自然地想到分治,考虑对 \([l,r]\) 进行分治,求出对于所有 \(l\le x\le y\le r\) 的点对 \((x,y)\ ...

  3. Luogu 4755 Beautiful Pair

    分治 + 主席树. 设$solve(l, r)$表示当前处理到$[l, r]$区间的情况,我们可以找到$[l, r]$中最大的一个数的位置$mid$,然后扫一半区间计算一下这个区间的答案. 注意,这时 ...

  4. P4755 Beautiful Pair

    题目 洛谷 做法 \(i≤x≤j,a[i]<\frac{a[x]}{a[j]}\) 考虑\(a[x]\)的贡献,单调栈预处理\(L,R\)能作为最大值的区间 枚举一端点,仅需另一端点满足条件即可 ...

  5. 洛谷$P4755\ Beautiful\ Pair$ 最大值分治

    正解:最大值分治 解题报告: 传送门$QwQ$ 昂考虑如果已经钦定了点$x$是这个$max$了,然后现在要求有多少对$[l,r]$满足$a_x=max\left\{a_i\right\},i\in[l ...

  6. P4755 Beautiful Pair (分治 + 主席树)

    题意:1e5的数组 计算有多少对 ai * aj <= max(ai ai+1...aj-1 aj) 题解:在处理这种涉及到区间极值的题时 好像是个套路分治 从级值中间分成两个区间 从区间短的那 ...

  7. 「LGR-049」洛谷7月月赛 D.Beautiful Pair

    「LGR-049」洛谷7月月赛 D.Beautiful Pair 题目大意 : 给出长度为 \(n\) 的序列,求满足 \(i \leq j\) 且 $a_i \times a_j \leq \max ...

  8. [luogu4755]Beautiful Pair

    [luogu4755]Beautiful Pair luogu 第一次写最大值分治感觉有点丑 每次找到最大值mid,扫小的一边,主席树查大的一边小于等于\(\frac{a[mid]}{a[i]}\)的 ...

  9. luoguP4755 Beautiful Pair

    https://www.luogu.org/problemnew/show/P4755 考虑分治,在 [l, r] 区间中用线段树找到最大的一个点,处理经过它的可行数对的个数,统计个数可以离线树状数组 ...

随机推荐

  1. Java实现字串统计

    对字符串的操作,无论再难的算法题,只要时间充足,相信每个同学都可以搞定. 但是浪费太多时间去搞一个逻辑算法没太大意义,学会学习,不但可以增长自己的知识,更可以节省时间,俗话说,一寸光阴一寸金,寸金难买 ...

  2. NetScaler VPX configration

    境搭建概述 本文主要介绍Netscaler的安装配置,以及与StoreFront相结合,实现外网访问内网资源.当用户访问时,Netscaler Gateway Virtual Server将把请求转给 ...

  3. Warning: setcookie() expects parameter 3 to be long, string given

    Warning: setcookie() expects parameter 3 to be long, string given 这个是我用php7.0会报这个错误, 切换低版本php5.6就ok

  4. Redis内存碎片率

    一. 内存碎片率mem_fragmentation_ratio = used_memory_rss / used_memoryused_memory :Redis使用其分配器分配的内存大小used_m ...

  5. Android新项目GBSS:第1篇 搭建开发环境

    最近接手一个Android新项目,之前也没做过这方面的开发,算是边学边干,这两天看了一下Android开发的书,大致入门了一点,今天把所需要的软件都下了下来,准备开工,先列一下开发环境: 所有的软件都 ...

  6. log4j日志配置文件

    log4j.properties: ### 设置### log4j.rootLogger = debug,stdout,D,E ### 输出信息到控制抬 ### log4j.appender.stdo ...

  7. phpfpm和nginx设置开机自动启动

    ## 添加执行权限 chmod a+x /etc/init.d/nginx chmod a+x /etc/init.d/php-fpm ## 加入服务 chkconfig --add nginx ch ...

  8. DDoS攻击的几种类型

    DDoS攻击的几种类型   随着网络攻击的简单化,如今DDoS攻击已经不止出现在大型网站中,就连很多中小型网站甚至是个人网站往往都可能面临着被DDoS攻击的的风险.或许很多站长对DDoS攻击并不是很了 ...

  9. CTF—攻防练习之HTTP—SQL注入(SSI注入)

    主机:192.168.32.152 靶机:192.168.32.161 ssI是赋予html静态页面的动态效果,通过ssi执行命令,返回对应的结果,若在网站目录中发现了.stm .shtm .shtm ...

  10. 直方图匹配原理与python、matlab实现

    直方图匹配本质上是让两幅图像的累积直方图尽量相似,累积直方图相似了,直方图也就相似了. 把原图像img的直方图匹配到参考图像ref的直方图,包括以下几个步骤: 1. 求出原图像img的累积直方图img ...