bzoj 3295: [Cqoi2011]动态逆序对(树套树)
来源:互联网 发布:网络推广类的书籍 编辑:程序博客网 时间:2024/06/05 00:56
3295: [Cqoi2011]动态逆序对
Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 3046 Solved: 970
[Submit][Status][Discuss]
Description
对于序列A,它的逆序对数定义为满足i<j,且Ai>Aj的数对(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
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)。
2
2
1
样例解释
(1,5,3,4,2)(1,3,4,2)(3,4,2)(3,2)(3)。
HINT
N<=100000 M<=50000
Source
题解:树状数组套线段树。
树状数组中维护的是逆序对的前缀和。
线段树中维护的是权值在树状数组当前点控制范围内区间中点的个数。
把数据离线,倒着搞,原本的删除操作就变成了插入操作。插入一个点那么这个新增节点的逆序对数就是数列中他前面的点中值比他大的(用前面的总点数-前面比他小的个数)+数列中他后面的节点中值比他小的。(统计的时候利用树状数组,就可以把比他小的所有权值通过lowbit的方式全部统计到)
#include<iostream>#include<cstdio>#include<cstring>#include<algorithm>#include<cmath>#define N 10000003#define M 100003#define LL long long using namespace std;int n,m,a[M],pos[M],pd[M],b[M],sum[N],c[M];int rs[N],ls[N],tot,p,s[M],sz,root[M];LL ans[M];int lowbit(int x){return (-x)&x;}int query(int now,int l,int r,int ll,int rr){//if (delta[now]) pushdown(now,l,r);if (!now) return 0;if (l>=ll&&r<=rr) return sum[now];int mid=(l+r)/2;int ans=0;if (ll<=mid) ans+=query(ls[now],l,mid,ll,rr);if (rr>mid) ans+=query(rs[now],mid+1,r,ll,rr);return ans;}void update(int x){sum[x]=sum[ls[x]]+sum[rs[x]];}void pointchange(int &k,int l,int r,int x){if (!k) k=++sz;if (l==r) { sum[k]++; return ; }int mid=(l+r)/2;if (x<=mid) pointchange(ls[k],l,mid,x);else pointchange(rs[k],mid+1,r,x);update(k);}int add(int x){int ans=0;for (int i=x;i;i-=lowbit(i)) ans+=s[i];return ans;}int solve(int x,int l,int r){int ans=0;for (int i=x;i;i-=lowbit(i)) ans+=query(root[i],1,n,l,r);return ans;}void insert(int x){int t=x; int q=add(pos[t]); int a=0; int b=0;a=solve(x,1,pos[t]-1);if (pos[t]+1<=n) b=solve(x,pos[t]+1,n);int k=(q-a)+b;for (int i=x;i<=n;i+=lowbit(i)){ pointchange(root[i],1,n,pos[t]); c[i]+=k; }}LL sum1(int x){LL ans=0;for (int i=x;i;i-=lowbit(i)) ans+=(LL)c[i];return ans;}void change(int x){for (int i=x;i<=n;i+=lowbit(i)) s[i]++;}int main(){scanf("%d%d",&n,&m);for (int i=1;i<=n;i++) scanf("%d",&a[i]),pos[a[i]]=i;for (int i=1;i<=m;i++) scanf("%d",&b[i]),pd[b[i]]=1;int maxn=0;for (int i=n;i>=1;i--) if (!pd[a[i]]) { insert(a[i]); change(pos[a[i]]); maxn=max(maxn,a[i]); }for (int i=m;i>=1;i--) { maxn=max(maxn,b[i]); insert(b[i]); change(pos[b[i]]); ans[i]=(LL)sum1(maxn); }for (int i=1;i<=m;i++) printf("%I64d\n",ans[i]);}
0 0
- bzoj 3295: [Cqoi2011]动态逆序对(树套树)
- [BZOJ 3295] Cqoi2011 动态逆序对
- 【BZOJ 3295】 [Cqoi2011]动态逆序对
- BZOJ 3295 [Cqoi2011]动态逆序对
- BZOJ 3295 [CQOI2011] 动态逆序对
- BZOJ 3295: [Cqoi2011]动态逆序对
- 【BZOJ 3295】[Cqoi2011]动态逆序对
- bzoj 3295 [Cqoi2011]动态逆序对
- BZOJ 3295 [Cqoi2011]动态逆序对
- [bzoj]3295: [Cqoi2011]动态逆序对 主席树
- 主席树初探 & bzoj 3295: [Cqoi2011] 动态逆序对 题解
- bzoj 3295: [Cqoi2011]动态逆序对 【cdq分治】
- bzoj 3295: [Cqoi2011]动态逆序对 cdq分治+树状数组
- BZOJ 3295: [Cqoi2011]动态逆序对 分块大法好
- BZOJ 3295 [Cqoi2011] 动态逆序对 CDQ分治题解
- 【BZOJ】3295 [Cqoi2011]动态逆序对 树状数组+线段树
- 【BZOJ】【P3295】【CQOI2011】【动态逆序对】【题解】【树套树】
- 3295: [Cqoi2011]动态逆序对
- Vim插件推荐
- 剑指offer:两个链表的第一个公共结点
- Mysql中文乱码问题完美解决方案
- Linux hostname主机名配置文件/etc/hosts详解
- 删除弹出确认框,成功后有提示
- bzoj 3295: [Cqoi2011]动态逆序对(树套树)
- 上机实现Linux命令整理(1)
- Android屏幕适配-----屏幕尺寸
- JVM之垃圾回收
- java观察者模式
- 设计模式之Adapter
- Java Web的数据库操作
- eclipse, Log4j配置(真心的详细~)
- JS实现监听事件,添加删除显示控件