Bzoj 3295: [Cqoi2011]动态逆序对 分块,树状数组,逆序对
3295: [Cqoi2011]动态逆序对
Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 2886 Solved: 924
[Submit][Status][Discuss]
Description
Input
Output
Sample Input
1
5
3
4
2
5
1
4
2
Sample Output
2
2
1
样例解释
(1,5,3,4,2)(1,3,4,2)(3,4,2)(3,2)(3)。
HINT
N<=100000 M<=50000
Source
#include<bits/stdc++.h>
using namespace std;
#define MAXN 100010
#define LL long long
int pos[MAXN],BIT[MAXN],a[MAXN],b[MAXN],wz[MAXN],block,n;
LL ans;
int read()
{
int s=,fh=;char ch=getchar();
while(ch<''||ch>''){if(ch=='-')fh=-;ch=getchar();}
while(ch>=''&&ch<=''){s=s*+(ch-'');ch=getchar();}
return s*fh;
}
int Lowbit(int o){return o&(-o);}
void Update(int o,int o1)
{
while(o<=n)
{
BIT[o]+=o1;
o+=Lowbit(o);
}
}
int Sum(int o)
{
int sum=;
while(o>)
{
sum+=BIT[o];
o-=Lowbit(o);
}
return sum;
}
void cl(int k)
{
int l=(k-)*block+,r=min(k*block,n),i;
for(i=l;i<=r;i++)b[i]=a[i];
sort(b+l,b+r+);
}
int ef(int k,int k1)//在第k块中寻找大于k1的第一个位置.
{
int l=(k-)*block+,r=min(k*block,n),ans1;
ans1=-;
while(l<=r)
{
int mid=(l+r)/;
if(b[mid]>k1){ans1=mid;r=mid-;}
else l=mid+;
}
return ans1;
}
int ef1(int k,int k1)
{
int l=(k-)*block+,r=min(k*block,n),ans1;
ans1=-;
while(l<=r)
{
int mid=(l+r)/;
if(b[mid]<k1){ans1=mid;l=mid+;}
else r=mid-;
}
return ans1;
}
void calc1(int l,int r,int D)
{
if(l>r)return;
int L=pos[l],R=pos[r],i,pd;
if(L==R)
{
for(i=l;i<=r;i++)if(a[i]>D&&a[i]!=-)ans--;
return;
}
for(i=l;i<=L*block;i++)if(a[i]>D&&a[i]!=-)ans--;
for(i=L+;i<=R-;i++)
{
pd=ef(i,D);
if(pd!=-)ans-=(LL)(min(i*block,n)-pd+);
}
//ans-=(ef(i,D)-ef(i,0)+1);
for(i=(R-)*block+;i<=r;i++)if(a[i]>D&&a[i]!=-)ans--;
}
void calc2(int l,int r,int D)
{
if(l>r)return;
int L=pos[l],R=pos[r],i,pd,pd1;
if(L==R)
{
for(i=l;i<=r;i++)if(a[i]<D&&a[i]!=-)ans--;
return;
}
for(i=l;i<=L*block;i++)if(a[i]<D&&a[i]!=-)ans--;
for(i=L+;i<=R-;i++)
{
//ans-=(ef1(i,D)-ef(i,0)+1);
pd=ef1(i,D);
if(b[pd]!=-)
{
pd1=ef(i,);
if(pd!=-)
{
if(pd1==-)pd1=;
ans-=(LL)(pd-pd1+);
}
}
}
for(i=(R-)*block+;i<=r;i++)if(a[i]<D&&a[i]!=-)ans--;
}
int main()
{
freopen("inverse.in","r",stdin);
freopen("inverse.out","w",stdout);
int m,i,M,del;
n=read();m=read();
for(i=;i<=n;i++)a[i]=read(),wz[a[i]]=i;
block=(int)sqrt(n);
for(i=;i<=n;i++)pos[i]=(int)(i-)/block+;
if(n%block==)M=n/block;
else M=n/block+;
memset(BIT,,sizeof(BIT));ans=;
for(i=n;i>=;i--)
{
ans+=Sum(a[i]-);
Update(a[i],);
}
for(i=;i<=M;i++)cl(i);
for(i=;i<=m;i++)
{
del=read();
printf("%lld\n",ans);
calc1(,wz[del]-,del);//去寻找 从位置1到要删除的数前一个位置 中有多少个数大于要删除的数,并把个数减去.
calc2(wz[del]+,n,del);//去寻找 从要删除的数后一个位置到位置n 中有多少个数小于要删除的数,并把个数减去.
a[wz[del]]=-;
cl(pos[wz[del]]);
//printf("%d\n",ans);
}
fclose(stdin);
fclose(stdout);
return ;
}
Bzoj 3295: [Cqoi2011]动态逆序对 分块,树状数组,逆序对的更多相关文章
- 【BZOJ 3295】动态逆序对 - 分块+树状数组
题目描述 给定一个1~n的序列,然后m次删除元素,每次删除之前询问逆序对的个数. 分析:分块+树状数组 (PS:本题的CDQ分治解法见下一篇) 首先将序列分成T块,每一块开一个树状数组,并且先把最初的 ...
- Bzoj 3289: Mato的文件管理 莫队,树状数组,逆序对,离散化,分块
3289: Mato的文件管理 Time Limit: 40 Sec Memory Limit: 128 MBSubmit: 1539 Solved: 665[Submit][Status][Di ...
- BZOJ 3744 Gty的妹子序列 (分块+树状数组+主席树)
题面传送门 题目大意:给你一个序列,多次询问,每次取出一段连续的子序列$[l,r]$,询问这段子序列的逆序对个数,强制在线 很熟悉的分块套路啊,和很多可持久化01Trie的题目类似,用分块预处理出贡献 ...
- BZOJ 3744 Gty的妹子序列 分块+树状数组
具体分析见 搬来大佬博客 时间复杂度 O(nnlogn)O(n\sqrt nlogn)O(nnlogn) CODE #include <cmath> #include <cctyp ...
- [树状数组+逆序对][NOIP2013]火柴排队
火柴排队 题目描述 涵涵有两盒火柴,每盒装有n根火柴,每根火柴都有一个高度.现在将每盒中的火柴各自排成一列,同一列火柴的高度互不相同,两列火柴之间的距离定义为:∑ (ai-bi)2,i=1,2,3,. ...
- AcWing 107. 超快速排序(归并排序 + 逆序对 or 树状数组)
在这个问题中,您必须分析特定的排序算法----超快速排序. 该算法通过交换两个相邻的序列元素来处理n个不同整数的序列,直到序列按升序排序. 对于输入序列9 1 0 5 4,超快速排序生成输出0 1 4 ...
- 【bzoj2141】排队 分块+树状数组
题目描述 排排坐,吃果果,生果甜嗦嗦,大家笑呵呵.你一个,我一个,大的分给你,小的留给我,吃完果果唱支歌,大家乐和和.红星幼儿园的小朋友们排起了长长地队伍,准备吃果果.不过因为小朋友们的身高有所区别, ...
- 【bzoj3744】Gty的妹子序列 分块+树状数组+主席树
题目描述 我早已习惯你不在身边, 人间四月天 寂寞断了弦. 回望身后蓝天, 跟再见说再见…… 某天,蒟蒻Autumn发现了从 Gty的妹子树(bzoj3720) 上掉落下来了许多妹子,他发现 她们排成 ...
- 【分块+树状数组】codechef November Challenge 2014 .Chef and Churu
https://www.codechef.com/problems/FNCS [题意] [思路] 把n个函数分成√n块,预处理出每块中各个点(n个)被块中函数(√n个)覆盖的次数 查询时求前缀和,对于 ...
随机推荐
- 数位DP入门Ural1057
CF一战让我觉得很疲倦,所以今天感觉很慢. 昨天解D题时候,因为太累,根本连题目都没看,今天看了之后感觉不会做,听闻是数位DP问题. 有某神说过,DP的功力建立在刷过的题上,我真的毫无功力可言. 介绍 ...
- SGU 解题报告
Volume 1 Volume 2
- PHP接口(interface)和抽象类(abstract)
interface 定义了一个接口类,它里面的方法其子类必须实现.接口是类的一个模板,其子类必须实现接口中定义的所有方法. interface User{ function getHeight ...
- Yii render和renderPartial的区别
以下由我们在信易网络公司开发项目的时候终结出的一些经验 在进行页面输出渲染的时候. 1.render 输出父模板的内容,将渲染的内容,嵌入父模板.|2.renderPartial 则不输出父模板的内容 ...
- MAC 终端快捷建
常用的快捷键: Ctrl + d 删除一个字符,相当于通常的Delete键(命令行若无所有字符,则相当于exit:处理多行标准输入时也表示eof) Ctrl + h 退格删 ...
- thinkphp多表查询
在学习thinkphp 的过程中,需要对多表进行操作,但是在实际过程中,总是遇到各种问题,所以写下这篇博文,作为自己的学习历程 在操作过程中,两表查询都没有问题,但是三表查询就开始出现问题 有以下三张 ...
- WPF 得一些问题汇总
1.CallMethodAction <TextBox Height="30" Name="txtUserName" Width="160&qu ...
- ARM中的PC和AXD的PC
R15 (PC)总是指向“正在取指”的指令,而不是指向“正在执行”的指令或正在“译码”的指令.一般来说,人们习惯性约定将“正在执行的指令作为参考点”,称之为当前第一条指令,因此PC 总是指向第三条指令 ...
- 定时显示div
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" " http://www.w3.org ...
- skiplist 跳表(2)-----细心学习
快速了解skiplist请看:skiplist 跳表(1) http://blog.sina.com.cn/s/blog_693f08470101n2lv.html 本周我要介绍的数据结构,是我非常非 ...