bzoj 3333: 排队计划 解决问题的方法

【原标题】

3333: 排队计划

Time Limit: 20 Sec Memory Limit: 128 MB

Submit: 161 Solved: 71

[Submit][Status]

Description

Input

Output

Sample Input

6 2

160 163 164 161 167 160

2

3

Sample Output

6

3

1

HINT

Source

wyx528命题

【分析】简述一下题目。N个数排成一列。每次指定一个位置P，然后把P~N中全部身高小于等于P的人都拎出来，排一遍序后再放进去。每次操作须要求一遍逆序对。

首先非常easy想到：每次选定一个位置P后，降低的逆序对数量就是P~N中满足要求的人原来构成的逆序对——并且当操作完毕之后。这些人将永远不会产生逆序对了。

考虑到这种性质，我们能够维护一个数的F[i]，表示从I这个位置開始的逆序对数量。（我们并不关心某个位置某一个时刻的值是多少）然后我们用一个线段树维护i~N中的最小值。每次操作寻找P~N的最小值。假设是满足要求的且位置是X，那么我们就把f[x]置为0，然后能够把x位置的数置为无穷大。之后继续寻找，直到找不到为止。由于每一个数仅仅会最多被更新一次，所以效率均摊是NlogN的。

【代码】

#include<cstdio>
#include<algorithm>
#define N 500005
#define INF 1000000005
#define Lo(x) (x&-x)
using namespace std;
typedef long long LL;LL ans=0;
struct Tree{int l,r,min,wh;}a[N*3];
int tree
;//pos P in array to the tree
int pos[N*4];//pos P int tree to the array
int f
,data
,c
,n,cnt,x,i,opt,P,ord,size,now,L,R;
struct array{int x,id;}uni
;
inline int cmp(const array &a,const array &b){return a.x<b.x;}
inline int ask(int x){int s=0;for (;x;x-=Lo(x)) s+=c[x];return s;}
inline void add(int x){for (;x<=n;x+=Lo(x)) c[x]++;}
void build(int k,int l,int r)
{
a[k].l=l;a[k].r=r;
if (l==r) {a[k].min=data[l];pos[k]=l;tree[l]=k;a[k].wh=k;return;}
int mid=(l+r)>>1;
if (l<=mid) build(k<<1,l,mid);
if (r>mid) build(k<<1|1,mid+1,r);
if (a[k<<1].min<a[k<<1|1].min) a[k].min=a[k<<1].min,a[k].wh=a[k<<1].wh;
else a[k].min=a[k<<1|1].min,a[k].wh=a[k<<1|1].wh;
}
int query(int k)
{
if (L<=a[k].l&&a[k].r<=R) return a[k].wh;
int mid=(a[k].l+a[k].r)>>1,resl=0,resr=0;
if (L<=mid) resl=query(k<<1);
if (R>mid) resr=query(k<<1|1);
if (!resl) return resr;if (!resr) return resl;
return (a[resl].min<a[resr].min)?

resl:resr;
}
void update(int k)
{
if (a[k].l==a[k].r) {a[k].min=INF;a[k].wh=k;return;}
int mid=(a[k].l+a[k].r)>>1;
if (ord<=mid) update(k<<1);else update(k<<1|1);
if (a[k<<1].min<a[k<<1|1].min) a[k].min=a[k<<1].min,a[k].wh=a[k<<1].wh;
else a[k].min=a[k<<1|1].min,a[k].wh=a[k<<1|1].wh;
}
int main()
{
read(n);read(opt);//读入优化略去
for (i=1;i<=n;i++) read(x),uni[i]=(array){x,i};
sort(uni+1,uni+n+1,cmp);
for (i=1;i<=n;i++)
data[uni[i].id]=(uni[i].x==uni[i-1].x)?cnt:++cnt;
for (i=n;i;i--)
f[i]=ask(data[i]-1),add(data[i]),ans+=(LL)f[i];
build(1,1,n);printf("%lld\n",ans);
while (opt--)
{
read(P);now=a[tree[P]].min;
while (now<INF)
{
L=P;R=n;ord=query(1);size=a[ord].min;
if (size>now) break;ord=pos[ord];
ans-=(LL)f[ord];f[ord]=0;update(1);
}
printf("%lld\n",ans);
}
return 0;
}