【GDOI2018模拟7.12】A 线段树+树状数组

这题用光了我所有的脑细胞，还没A，成功导致我今天比赛瞬间爆炸。。

题意略。

首先肯定不能用归并无脑求。。。否则你怎么处理后面一问，这里明显提示我们第一问的某些东西可以用到第二问中。
先用树状数组求出总逆序对数，同时算出每一位的贡献。
然后对于第二问我们离线处理，先把所有的处理记录一下，记录每一位被更新的时间。
然后我们按照序列的顺序处理。
对于一个数来说，他的贡献变为0的时刻就是在他前面比他大的而且最早进行的操作对应的时刻。
那么我们用线段树维护一下一段区间内被修改的最小时间，查询的时候查询当前被修改的最小时间，那么明显这个位上的贡献都可以加到这个时间点去，那么按照序列的顺序做下去就好了，最后按照时间输出每个时间点上的贡献值。

#include<cstdio>
#include<algorithm>
#include<cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=1e6+10;
const int inf=1e9;
int n,m;
typedef long long ll;
ll ans;
int val[N],tim[N],f[N];
ll value[N];
int mn[N*4];
struct node
{
    int x,id;
}a[N];
int c[N],tot,b[N];
bool cmp(node a,node b)
{
    return a.x<b.x;
}
inline int lowbit(int x)
{
    return x&(-x);
}
inline void add(int x,int v)
{
    while (x<N)
    {
        val[x]+=v;
        x+=lowbit(x);
    }
}
inline int find(int x)
{
    int ans=0;
    while (x>0)
    {
        ans+=val[x];
        x-=lowbit(x);
    }
    return ans;
}

inline void change(int x,int l,int r,int pos,int v)
{
    if (l==r)
    {
        mn[x]=min(mn[x],v);
        return;
    }
    int mid=(l+r)>>1;
    if (pos<=mid)change(x<<1,l,mid,pos,v);
    else change(x<<1|1,mid+1,r,pos,v);
    mn[x]=min(mn[x<<1],mn[x<<1|1]);
}
inline int query(int x,int l,int r,int l1,int r1)
{
    if (l1==l&&r1==r)return mn[x];
    int mid=(l+r)>>1;
    if (r1<=mid)return query(x<<1,l,mid,l1,r1);
    else if (l1>mid)return query(x<<1|1,mid+1,r,l1,r1);
    else return min(query(x<<1,l,mid,l1,mid),query(x<<1|1,mid+1,r,mid+1,r1));
}
int main()
{
    scanf("%d%d",&n,&m);
    fo(i,1,n)
    {
        scanf("%d",&a[i].x);
        a[i].id=i;
        tim[i]=inf;
    }
    fo(i,0,m-1)
    {
        int x;
        scanf("%d",&x);
        if (tim[x]==inf)
        {
            tim[x]=i;
        }
    }
//  printf("fuckk%d\n",tim[1]);
    fd(i,n,1)
    {
        f[i]=find(a[i].x);
        add(a[i].x+1,1);
        ans+=f[i];
        //printf("%d\n",f[i]);
    }
    printf("%lld\n",ans);
    sort(a+1,a+1+n,cmp);
    tot=0;
    fo(i,1,n)
    if (a[i].x!=a[i-1].x)b[a[i].id]=++tot;
    else b[a[i].id]=tot;
    //fo(i,1,n)printf("%d\n",b[i]);
    memset(mn,127,sizeof(mn));
    fo(i,1,n)
    {
        change(1,1,tot,b[i],tim[i]);
        int time=query(1,1,tot,b[i],tot);
        if (time<=m)value[time]+=f[i];
    }
    fo(i,0,m-1)
    {
        ans=ans-value[i];
        printf("%lld\n",ans);
    }
}