bzoj1831 AHOI2008

http://www.elijahqi.win/archives/563

1831: [AHOI2008]逆序对

Time Limit: 10 Sec Memory Limit: 64 MB
Submit: 683 Solved: 485
[Submit][Status][Discuss]

Description

小可可和小卡卡想到Y岛上旅游，但是他们不知道Y岛有多远。好在，他们找到一本古老的书，上面是这样说的： 下面是N个正整数，每个都在1~K之间。如果有两个数A和B，A在B左边且A大于B，我们就称这两个数为一个“逆序对”。你数一数下面的数字里有多少个逆序对，你就知道Y岛离这里的距离是多少千米了。 比如说，4 2 1 3 3里面包含了5个逆序对：(4, 2), (4, 1), (4, 3), (4, 3), (2, 1)。 可惜的是，由于年代久远，这些数字里有一部分已经模糊不清了，为了方便记录，小可可用“-1”表示它们。比如说，4 2 -1 -1 3 可能原来是4 2 1 3 3，也可能是4 2 4 4 3，也可能是别的样子。 小可可希望知道，根据他们看清楚的这部分数字，能不能推断出这些数字里最少能有多少个逆序对。
Input

第一行两个正整数N和K。第二行N个整数，每个都是-1或是一个在1~K之间的数。
Output

一个正整数，即这些数字里最少的逆序对个数。
Sample Input

5 4
4 2 -1 -1 3
Sample Output

4

HINT

4 2 4 4 3中有4个逆序对。当然，也存在其它方案得到4个逆序对。

数据范围：
100%的数据中，N<=10000，K<=100。
60%的数据中，N<=100。
40%的数据中，-1出现不超过两次。

Source

Day1

有一些关键性质是从别人那里学习来的 代码中的解释是个人对这题的一个理解吧，蒟蒻经常出错的请指教orz

首先，我们考虑填入的数需要满足的性质。

对于两个空位a和b，分别填入数x和y，且x

#include<cstdio>
#define inf 0x7fffffff
#define N 11000
#define N1 110
inline int read(){
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9') {if(ch=='-')f=-1;ch=getchar();}
    while (ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
inline int min(int x,int y){
    return x<y?x:y;
}
int s[N][N1],t[N][N1],f[N][N1],g[N][N1],ans,n,k,a[N],b[N],num;
int main(){
    freopen("bzoj1831.in","r",stdin);
    n=read();k=read();
    for (int i=1;i<=n;++i) {a[i]=read();if(a[i]==-1) b[++num]=i;}
    //s数组为预处理在1~i位置填j大小的数预处理有多少比j大的数
    //t为预处理n~i位置有多少比j小的数不包含没有填写的-1
    for (int i=1;i<=n;++i) for (int j=1;j<=k;++j) s[i][j]=s[i-1][j]+(a[i]>j) ;
    for (int i=n;i>=1;--i) for (int j=1;j<=k;++j) t[i][j]=t[i+1][j]+(a[i]>-1&&a[i]<j);
    //f[i][j]  表示在第i个空缺的地方填上j这个数值 
    for (int i=1;i<=num;++i){
        //那么在i这个缺口填j数值的时候逆序对是
        f[i][1]=f[i-1][1]+s[b[i]][1];g[i][1]=f[i][1]; 
        for (int j=2;j<=k;++j){
            f[i][j]=g[i-1][j]+s[b[i]][j]+t[b[i]][j];//f[i][j]=上一个坑填J形成的逆序对最小值+这次填j形成的代价 
            g[i][j]=min(g[i][j-1],f[i][j]);//g[i][j]第i个坑，表示1~j这么多填坑方案 哪种方案才最小 
        }
    }
    int ans=inf;
    for (int i=1;i<=k;++i) ans=min(ans,f[num][i]);
    for (int i=1;i<=n;++i) if (a[i]!=-1) ans+=s[i][a[i]];
    printf("%d",ans);
    return 0;
}