bzoj 3295 [Cqoi2011]动态逆序对

3295: [Cqoi2011]动态逆序对

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 4785  Solved: 1613
[Submit][Status][Discuss]

Description

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

Input

输入第一行包含两个整数n和m，即初始元素的个数和删除的元素个数。以下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

【分析】

吼题啊...

位置(w)一个维度，权值(a)一个维度，把位置排序，权值CDQ一下，时间线(t)用树状数组统计一下

题目就是要求对于一个t0，满足t<t0,w<w0,a>a0的个数和t<t0,w>w0,a<a0的个数和...这两个做两遍CDQ分治就好了...(类似陌上花开QAQ)

昨晚肝了半晚上没肝出来...诶蠢

对了我写的CDQ自带大常数（QAQ 2000ms+）

【代码】

#include<algorithm>#include<iostream>#include<cstring>#include<cstdio>#define ll long long#define fo(i,j,k) for(i=j;i<=k;i++)using namespace std;const int mxn=100005;int n,m,T,num;int to[mxn],c[mxn],tim[mxn];ll res[mxn];struct node{    int a,w,t;}p[mxn],tmp[mxn];inline bool comp(node x,node y){    return x.w<y.w;}inline int lowbit(int x) {return x&-x;}inline void add(int x,int v){for(int i=x;i<=n;i+=lowbit(i))  if(tim[i]!=T) c[i]=v,tim[i]=T;  else c[i]+=v; } inline int getsum(int x){int sum=0;for(int i=x;i>=1;i-=lowbit(i))  if(tim[i]==T) sum+=c[i];return sum;}inline void CDQ(int l,int r){    if(l==r) return;    int i,j,mid=l+r>>1,l1=l,l2=mid+1;    fo(i,l,r)      if(p[i].a<=mid) tmp[l1++]=p[i];      else tmp[l2++]=p[i];    memcpy(p+l,tmp+l,sizeof(p[0])*(r-l+1));    CDQ(l,mid);T++;j=l;    fo(i,mid+1,r)    {    for(j;j<=mid && p[j].w<p[i].w;j++)      add(p[j].t,1);    res[p[i].t]+=getsum(p[i].t);}CDQ(mid+1,r);l1=l,l2=mid+1;fo(i,l,r){if((comp(p[l1],p[l2]) || l2>r) && l1<=mid)  tmp[i]=p[l1++];else   tmp[i]=p[l2++];}memcpy(p+l,tmp+l,sizeof(p[0])*(r-l+1));}int main(){    int i,j,c;    scanf("%d%d",&n,&m);    for(i=n;i>=1;i--)      scanf("%d",&p[i].a),p[i].w=i,to[p[i].a]=i;    for(i=m;i>=1;i--)     {        scanf("%d",&c);        p[to[c]].t=i-m+n;    }    fo(i,1,n) if(!p[i].t) p[i].t=++num;    sort(p+1,p+n+1,comp);    CDQ(1,n);    fo(i,1,n) p[i].w=n-p[i].w+1,p[i].a=n-p[i].a+1;    sort(p+1,p+n+1,comp);    CDQ(1,n);    fo(i,1,n) res[i]+=res[i-1];    for(i=n;i>=n-m+1;i--) printf("%lld\n",res[i]);    return 0;}