3295: [Cqoi2011]动态逆序对

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 2886  Solved: 924
[Submit][Status][Discuss]

Description

对于序列A,它的逆序对数定义为满足i<j,且Ai>Aj的数对(i,j)的个数。给1到n的一个排列,按照某种顺序依次删除m个元素,你的任务是在每次删除一个元素之前统计整个序列的逆序对数。

Input

输入第一行包含两个整数nm,即初始元素的个数和删除的元素个数。以下n行每行包含一个1到n之间的正整数,即初始排列。以下m行每行一个正整数,依次为每次删除的元素。
 

Output

 
输出包含m行,依次为删除每个元素之前,逆序对的个数。

Sample Input

5 4
1
5
3
4
2
5
1
4
2

Sample Output

5
2
2
1

样例解释
(1,5,3,4,2)(1,3,4,2)(3,4,2)(3,2)(3)。

HINT

N<=100000 M<=50000

Source

 题解:
分块+树状数组+逆序对
好像正解是CDQ分治 QAQ
然而不会。。。
胡搞了一个分块,记录一下每个数出现的位置。然后,考虑到每次删除,对答案有影响的只有当前这个数前头比它大的个数,和后头比他小的个数。所以就分块统计有多少个。但是,如何删除呢?
我是在更新完答案后,把当前的这个数放成-1。然后对于统计这个数前头比它大的数的个数没有影响(因为-1比任何数都小,不会被统计上)。但对于后头比它小的个数就不一样了。但,-1比任何数都小,所以它一定在每个排好序的块的最前端的一段。我们只需要二分去找第一个不是-1的位置,然后就可以统计了。(各种二分。。。)
至于速度嘛,rank榜上倒数。。。QAQ赶快去学CDQ分治。。。
 #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]动态逆序对 分块,树状数组,逆序对的更多相关文章

  1. 【BZOJ 3295】动态逆序对 - 分块+树状数组

    题目描述 给定一个1~n的序列,然后m次删除元素,每次删除之前询问逆序对的个数. 分析:分块+树状数组 (PS:本题的CDQ分治解法见下一篇) 首先将序列分成T块,每一块开一个树状数组,并且先把最初的 ...

  2. Bzoj 3289: Mato的文件管理 莫队,树状数组,逆序对,离散化,分块

    3289: Mato的文件管理 Time Limit: 40 Sec  Memory Limit: 128 MBSubmit: 1539  Solved: 665[Submit][Status][Di ...

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

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

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

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

  5. [树状数组+逆序对][NOIP2013]火柴排队

    火柴排队 题目描述 涵涵有两盒火柴,每盒装有n根火柴,每根火柴都有一个高度.现在将每盒中的火柴各自排成一列,同一列火柴的高度互不相同,两列火柴之间的距离定义为:∑ (ai-bi)2,i=1,2,3,. ...

  6. AcWing 107. 超快速排序(归并排序 + 逆序对 or 树状数组)

    在这个问题中,您必须分析特定的排序算法----超快速排序. 该算法通过交换两个相邻的序列元素来处理n个不同整数的序列,直到序列按升序排序. 对于输入序列9 1 0 5 4,超快速排序生成输出0 1 4 ...

  7. 【bzoj2141】排队 分块+树状数组

    题目描述 排排坐,吃果果,生果甜嗦嗦,大家笑呵呵.你一个,我一个,大的分给你,小的留给我,吃完果果唱支歌,大家乐和和.红星幼儿园的小朋友们排起了长长地队伍,准备吃果果.不过因为小朋友们的身高有所区别, ...

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

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

  9. 【分块+树状数组】codechef November Challenge 2014 .Chef and Churu

    https://www.codechef.com/problems/FNCS [题意] [思路] 把n个函数分成√n块,预处理出每块中各个点(n个)被块中函数(√n个)覆盖的次数 查询时求前缀和,对于 ...

随机推荐

  1. 数位DP入门Ural1057

    CF一战让我觉得很疲倦,所以今天感觉很慢. 昨天解D题时候,因为太累,根本连题目都没看,今天看了之后感觉不会做,听闻是数位DP问题. 有某神说过,DP的功力建立在刷过的题上,我真的毫无功力可言. 介绍 ...

  2. SGU 解题报告

    Volume 1 Volume 2

  3. PHP接口(interface)和抽象类(abstract)

    interface 定义了一个接口类,它里面的方法其子类必须实现.接口是类的一个模板,其子类必须实现接口中定义的所有方法. interface User{     function getHeight ...

  4. Yii render和renderPartial的区别

    以下由我们在信易网络公司开发项目的时候终结出的一些经验 在进行页面输出渲染的时候. 1.render 输出父模板的内容,将渲染的内容,嵌入父模板.|2.renderPartial 则不输出父模板的内容 ...

  5. MAC 终端快捷建

    常用的快捷键: Ctrl + d        删除一个字符,相当于通常的Delete键(命令行若无所有字符,则相当于exit:处理多行标准输入时也表示eof) Ctrl + h        退格删 ...

  6. thinkphp多表查询

    在学习thinkphp 的过程中,需要对多表进行操作,但是在实际过程中,总是遇到各种问题,所以写下这篇博文,作为自己的学习历程 在操作过程中,两表查询都没有问题,但是三表查询就开始出现问题 有以下三张 ...

  7. WPF 得一些问题汇总

    1.CallMethodAction <TextBox Height="30" Name="txtUserName" Width="160&qu ...

  8. ARM中的PC和AXD的PC

    R15 (PC)总是指向“正在取指”的指令,而不是指向“正在执行”的指令或正在“译码”的指令.一般来说,人们习惯性约定将“正在执行的指令作为参考点”,称之为当前第一条指令,因此PC 总是指向第三条指令 ...

  9. 定时显示div

    <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" " http://www.w3.org ...

  10. skiplist 跳表(2)-----细心学习

    快速了解skiplist请看:skiplist 跳表(1) http://blog.sina.com.cn/s/blog_693f08470101n2lv.html 本周我要介绍的数据结构,是我非常非 ...