Minimum Inversion Number [树状数组]

这也是一道求逆序对的题，不过是求所有数列中的逆序对的最小值，1A了好开心~~~

想了挺久的，首先求升序的每个数的顺序对，即i<j,ai<aj

例如：            1 3 6 9 0 8 5 7 4 2               

 顺序对：        0 1 2 3 0 4 3 5 3 2            p1

在求逆序对： 1 2 4 6 0 4 2 2 1 0             p2

代码如下：

for(int i=0; i<n; i++)     //顺序对
        {
            scanf("%d",&a[i]);
            a[i]++;
            p1[i]=sum(a[i]-1,s1);
            add(a[i],1,s1);
        }
        for(int i=n-1; i>=0; i--)    //逆序对
        {
            ans+=sum(a[i]-1,s);
            p2[i]=sum(a[i]-1,s);
            add(a[i],1,s);
        }

每进行一次数组操作时，数组第一个元素放到最后，原本这位的逆序对为p，换后变成0，除了本位外总体的逆序对增加(n-p-1)。（n-p-1）等于其他元素增加的总和。这样的话，互换之后只要知道第一个元素本身的逆序对是多少就可以了，每次在原来的逆序对总数上加上-p+(n-p-1)。

那么每次互换之后第一个元素的逆序对是多少呢？就是相应的p1+p2，应为互换之后原来的顺序对就变成了逆序对了。

感觉自己说的好乱，表达能力真是我一个软肋~~

不过看别人的讲解，发现自己写的有点多余，想得有点太复杂了，其实在求出原始数组的逆序对数之后套一个公式就可以求出最小值了，对于每一个数组的第一个元素ai，ai之后有（n-ai）个元素比他大，（ai-1）个元素比他小，即他的逆序对是（ai-1），互换之后增加了（n-ai）个逆序对，所以得出sum=sum-(ai-1)+(n-ai)。这种方法比较直接简单，看来还是没有抓住主要特征。

相关链接：

Minimum Inversion Number

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>

using namespace std;
#define N 5010
int a[N];
int lowbit(int x)
{
    return x&-x;
}
void add(int x,int p,int *s)
{
    for(int i=x; i<N; i+=lowbit(i))
    {
        s[i]+=p;
    }
}
int sum(int x,int *s)
{
    int ans=0;
    for(int i=x; i>0; i-=lowbit(i))
    {
        ans+=s[i];
    }
    return ans;
}
int main()
{
    int n,s[N],s1[N],p1[N],p2[N];
    while(~scanf("%d",&n))
    {
        int ans=0;
        memset(s,0,sizeof(s));
        memset(s1,0,sizeof(s1));
        for(int i=0; i<n; i++)     //顺序对
        {
            scanf("%d",&a[i]);
            a[i]++;
            p1[i]=sum(a[i]-1,s1);
            add(a[i],1,s1);
        }
        for(int i=n-1; i>=0; i--)    //逆序对
        {
            ans+=sum(a[i]-1,s);
            p2[i]=sum(a[i]-1,s);
            add(a[i],1,s);
        }
        //for(int i=0;i<n;i++)
        //cout<<sum(a[i]-1,s)<<endl;
        int p=p2[0];
        int maxn=ans;
        //cout<<ans<<" "<<p<<endl;
        for(int i=0;i<n;i++)
        {
            if(ans<maxn)
            maxn=ans;
            ans=ans-p+(n-p-1);
            p=p2[i+1]+p1[i+1];
        }
        cout<<maxn<<endl;
    }
         return 0;
}