bzoj1831 逆序对 (dp+树状数组)

注意到，所有的-1应该是一个不降的序列，否则不会更优
那就先求出来不是-1的的逆序对个数，然后设f[i][j]表示第i个-1放成j的前i个-1带来的最小逆序对数量
这个可以树状数组来求

#include<bits/stdc++.h>
#define pa pair<int,int>
#define CLR(a,x) memset(a,x,sizeof(a))
using namespace std;
typedef long long ll;
const int maxn=1e4+10,maxk=110;

inline ll rd(){
    ll x=0;char c=getchar();int neg=1;
    while(c<'0'||c>'9'){if(c=='-') neg=-1;c=getchar();}
    while(c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
    return x*neg;
}


int f[maxn][maxk],tr[maxn],a[maxn],b[maxn],N,K,M;

inline int lowbit(int x){return x&(-x);}
inline void add(int x,int y){
    for(;x<=K;x+=lowbit(x)) tr[x]+=y;
}
inline int query(int x){
    int re=0;for(;x;x-=lowbit(x)) re+=tr[x];return re;
}

int main(){
    //freopen("","r",stdin);
    int i,j,k;
    N=rd(),M=0,K=rd();
    int ans=0;
    for(i=1;i<=N;i++){
        a[i]=rd();
        if(a[i]==-1) b[++M]=i;
        else{
            ans+=query(K)-query(a[i]);
            add(a[i],1);
        }
    }
    CLR(tr,0);
    for(i=1,j=1;i<=M;i++){
        for(;j<=b[i];j++){
            if(a[j]!=-1) add(a[j],1);
        }
        for(k=1;k<=K;k++)
            f[i][k]+=query(K)-query(k);
    }
    CLR(tr,0);
    for(i=M,j=N;i;i--){
        for(;j>=b[i];j--){
            if(a[j]!=-1) add(a[j],1);
        }
        for(k=1;k<=K;k++)
            f[i][k]+=query(k-1);
    }
    for(i=1;i<=M+1;i++){
        for(j=1;j<=K;j++){
            int mi=1e9;
            for(k=1;k<=j;k++){
                mi=min(mi,f[i-1][k]);
            }
            f[i][j]+=mi;
        }
    }
    printf("%d
",f[M+1][K]+ans);
    return 0;
}