题意

给出n个数,要求支持单点修改和区间逆序对,强制在线。
n,m<=50000

题解

和不带修改差不多,预处理出smaller[i][j]代表前i块小于j的数的数量,但不能用f[i][j]代表第i块到第j块逆序对的数量,这样不好维护。

我们用f[i][j]代表从第i块选出一个元素与从第j块选出一个元素组成逆序对的数量,维护时最多修改根号n个f数组,查询时用前缀和起到与不带修改时f数组的作用。

其他部分和不带修改时差不多。

然后问题就解决了。

 #include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
const int N=;
int n,Block,a[N],block[N],L[],R[],tr1[][],tr[][N],ans,m;
int lowbit(int x){
return x&-x;
}
void add(int id,int x,int w){
for(int i=x;i<=n;i+=lowbit(i)){
tr[id][i]+=w;
}
}
int getsum(int id,int x){
int tmp=;
for(int i=x;i;i-=lowbit(i)){
tmp+=tr[id][i];
}
return tmp;
}
void add1(int id,int x,int w){
for(int i=x;i<=block[n];i+=lowbit(i)){
tr1[id][i]+=w;
}
}
int getsum1(int id,int x){
int tmp=;
for(int i=x;i;i-=lowbit(i)){
tmp+=tr1[id][i];
}
return tmp;
}
int main(){
scanf("%d",&n);Block=sqrt(n);
for(int i=;i<=n;i++){
scanf("%d",&a[i]);
block[i]=(i-)/Block+;
if(!L[block[i]])L[block[i]]=i;
R[block[i]]=i;
}
for(int i=;i<=block[n];i++)
for(int j=L[i];j<=n;j++){
add(i,a[j],);
}
for(int i=;i<=block[n];i++){
for(int j=L[i];j<=R[i];j++){
add(,a[j],);
add1(i,i,getsum(,n)-getsum(,a[j]));
}
for(int j=R[i]+;j<=n;j++){
add1(i,block[j],getsum(,n)-getsum(,a[j]));
}
for(int j=L[i];j<=R[i];j++){
add(,a[j],-);
}
}
scanf("%d",&m);
while(m--){
int k,x,y;
scanf("%d%d%d",&k,&x,&y);
x^=ans;y^=ans;
if(k==){
for(int i=;i<=block[x]-;i++){
int size=R[i]-L[i]+;
add1(i,block[x],size-(getsum(i,y)-getsum(i+,y))-(size-(getsum(i,a[x])-getsum(i+,a[x]))));
}
for(int i=block[x]+;i<=block[n];i++){
add1(block[x],i,getsum(i,y-)-getsum(i+,y-)-(getsum(i,a[x]-)-getsum(i+,a[x]-)));
}
for(int i=L[block[x]];i<=R[block[x]];i++){
add(,a[i],);
add1(block[x],block[x],-(getsum(,n)-getsum(,a[i])));
}
for(int i=L[block[x]];i<=R[block[x]];i++){
add(,a[i],-);
}
for(int i=;i<=block[x];i++){
add(i,a[x],-);add(i,y,);
}
a[x]=y;
for(int i=L[block[x]];i<=R[block[x]];i++){
add(,a[i],);
add1(block[x],block[x],getsum(,n)-getsum(,a[i]));
}
for(int i=L[block[x]];i<=R[block[x]];i++){
add(,a[i],-);
}
}
else{
if(block[x]+>=block[y]){
ans=;
for(int i=x;i<=y;i++){
add(,a[i],);
ans+=getsum(,n)-getsum(,a[i]);
}
for(int i=x;i<=y;i++){
add(,a[i],-);
}
printf("%d\n",ans);
}
else{
ans=;
for(int i=block[x]+;i<=block[y]-;i++){
ans+=getsum1(i,block[y]-);
}
for(int i=x;i<=R[block[x]];i++){
add(,a[i],);
ans+=getsum(,n)-getsum(,a[i]);
ans+=getsum(block[x]+,a[i]-)-getsum(block[y],a[i]-);
}
for(int i=L[block[y]];i<=y;i++){
add(,a[i],);
ans+=getsum(,n)-getsum(,a[i]);
ans+=getsum(block[x]+,n)-getsum(block[y],n)-(getsum(block[x]+,a[i])-getsum(block[y],a[i]));
}
for(int i=x;i<=R[block[x]];i++){
add(,a[i],-);
}
for(int i=L[block[y]];i<=y;i++){
add(,a[i],-);
}
printf("%d\n",ans);
}
}
}
return ;
}

BZOJ 3787 Gty的文艺妹子序列(分块+树状数组+前缀和)的更多相关文章

  1. BZOJ 3787: Gty的文艺妹子序列 [分块 树状数组!]

    传送门 题意:单点修改,询问区间内逆序对数,强制在线 看到加了!就说明花了不少时间.... 如果和上题一样预处理信息,用$f[i][j]$表示块i到j的逆序对数 强行修改的话,每个修改最多会修改$(\ ...

  2. BZOJ3787:Gty的文艺妹子序列(分块,树状数组)

    Description Autumn终于会求区间逆序对了!Bakser神犇决定再考验一下他,他说道: “在Gty的妹子序列里,某个妹子的美丽度可也是会变化的呢.你还能求出某个区间中妹子们美丽度的逆序对 ...

  3. BZOJ3787 gty的文艺妹子序列 【树状数组】【分块】

    题目分析: 首先这种乱七八糟的题目就分块.然后考虑逆序对的统计. 一是块内的,二是块之间的,三是一个块内一个块外,四是都在块外. 令分块大小为$S$. 块内的容易维护,单次维护时间是$O(S)$. 块 ...

  4. BZOJ 3787: Gty的文艺妹子序列

    3787: Gty的文艺妹子序列 Time Limit: 50 Sec  Memory Limit: 256 MBSubmit: 186  Solved: 58[Submit][Status][Dis ...

  5. 【bzoj3744】Gty的妹子序列 分块+树状数组+主席树

    题目描述 我早已习惯你不在身边, 人间四月天 寂寞断了弦. 回望身后蓝天, 跟再见说再见…… 某天,蒟蒻Autumn发现了从 Gty的妹子树(bzoj3720) 上掉落下来了许多妹子,他发现 她们排成 ...

  6. BZOJ 3744 Gty的妹子序列 (分块+树状数组+主席树)

    题面传送门 题目大意:给你一个序列,多次询问,每次取出一段连续的子序列$[l,r]$,询问这段子序列的逆序对个数,强制在线 很熟悉的分块套路啊,和很多可持久化01Trie的题目类似,用分块预处理出贡献 ...

  7. 【BZOJ3744】Gty的妹子序列 分块+树状数组

    [BZOJ3744]Gty的妹子序列 Description 我早已习惯你不在身边, 人间四月天 寂寞断了弦. 回望身后蓝天, 跟再见说再见…… 某天,蒟蒻Autumn发现了从 Gty的妹子树(bzo ...

  8. BZOJ 3744 Gty的妹子序列 分块+树状数组

    具体分析见 搬来大佬博客 时间复杂度 O(nnlogn)O(n\sqrt nlogn)O(nn​logn) CODE #include <cmath> #include <cctyp ...

  9. BZOJ 4765 普通计算姬 dfs序+分块+树状数组+好题!!!

    真是道好题...感到灵魂的升华... 按dfs序建树状数组,拿前缀和去求解散块: 按点的标号分块,分成一个个区间,记录区间子树和 的 总和... 具体地,需要记录每个点u修改后,对每一个块i的贡献,记 ...

随机推荐

  1. keepalived+双主架构部署

    在高可用集群环境中,keepalived使用的是VIP,利用keepalived自带的服务监控功能和自定义脚本来实现MYSQL故障时自带切换. Keepalived基于VRRP协议,虚拟冗余路由协议, ...

  2. BZOJ 3165 李超线段树

    思路: 李超线段树 我是把线段转成斜率的形式搞得 不知道有没有更简单的方法 //By SiriusRen #include <cmath> #include <cstdio> ...

  3. 内存文件系统:tachyon(现在叫Alluxio)

    此文于2015 年 8 月 10 日发布 Tachyon 是什么 Tachyon 是 AMPLab 开发的一款内存分布式文件系统.它介于计算层和存储层之间,可以简单的理解为存储层在内存内的一个 Cac ...

  4. NPOI简单的给某个单元格字体设置颜色

    参考文档有: https://www.cnblogs.com/gossip/p/4307486.html https://bbs.csdn.net/topics/391042064?page=1 效果 ...

  5. 前端那些事之----jQuery

    1.jquery是什么     一个js的框架,可以方便的使用js 2 什么是jQuery对象     是由jQuery封装后的DOM对象     注意:与DOM对象的方法不同,不可以混用,但是可以相 ...

  6. FFT&NTT学习笔记

    具体原理就不讲了qwq,毕竟证明我也不太懂 FFT(快速傅立叶变换)&NTT(快速数论变换) FFT //求多项式乘积 //要求多项式A和多项式B的积多项式C //具体操作就是 //DFT(A ...

  7. 你可能需要了解下Laravel集合

    前言 集合通过 Illuminate\Support\Collection 进行实例,Laravel的内核大部分的参数传递都用到了集合,但这并不代表集合就是好的.Laravel作为快捷并优雅的开发框架 ...

  8. Openstack API 开发 快速入门

    Openstack 做为流行的开源云计算平台,其最大特性是利用其提供的基础设施API,让我们可以以软件的方式来动态管理IAAS资源.Openstack 提供的api是流行的Rest API.     ...

  9. CentOS 7.2 (mini) 里iptables防火墙怎么关闭?

    centos从7开始默认用的是firewalld,这个是基于iptables的,虽然有iptables的核心,但是iptables的服务是没安装的.所以你只要停止firewalld服务即可:sudo ...

  10. code vs 1216 跳马问题

    1216 跳马问题  时间限制: 1 s  空间限制: 128000 KB  题目等级 : 黄金 Gold 题解       题目描述 Description 题目 输入描述 Input Descri ...