hdu1394(暴力/线段树/归并排序求逆序对的个数)

题目链接：Minimum Inversion Number

题目大意：给出一个序列a1, a2, ..., an-1, an （大小为0<=ai-1<=n-1），

如下所示，

a1, a2, ..., an-1, an 
a2, a3, ..., an, a1 
a3, a4, ..., an, a1, a2 
...
an, a1, a2, ..., an-1

一个逆序对表示为ai>aj&&i<j

求上述排列中逆序对最小的数量

分析 ：

对某个数x,当它排在第一位的时候，后面比她小的数有x个即（0，1，...x-1），设此时逆序对数为tot，若最前面一个数为x，序列中比x 小的数为[0, x-1], 共x个，比x大的数为[x+1, n-1],  共tot-x-1个

当把它放到末尾，原先的逆序对变为顺序的，那么逆序对数要-x，并+tot-x-1；

由以上分析得，主要每次将a[i]移动到序列头，每次移动-想-x+tot-x-1

1、暴力

因为n<=5000，暴力解法n^2满足要求

#include <stdio.h>

const int maxn=5005;
int a[maxn];
int main()
{
    //freopen("in.txt","r",stdin);
    int n,tot,ans;
    while(scanf("%d",&n)!=EOF){
        tot=0;;
        for(int i=0;i<n;i++)scanf("%d",&a[i]);
        for(int i=0;i<n;i++){
            for(int j=i+1;j<n;j++){
                if(a[i]>a[j])tot++;
            }
        }
        ans=tot;
        for(int i=0;i<n;i++){
            tot=tot+n-1-a[i]-a[i];
            if(ans>tot)ans=tot;
        }
        printf("%d\n",ans);
    }
    return 0;
}

2、线段树

     

叶子节点代表的左右端点a[i]，节点中num维护是否插入了a[i]，按顺序插入a[i]，每次插入查询区间[a[i]，n-1]中的点数，因为这个区间内的点都比seq[i]先插入且比seq[i]大，所以，这个 区间内的点的个数就等于a[i]的逆序数。

代码

#include <stdio.h>

const int maxn=5005;
struct node
{
    int l,r;
    int sum;
}segTree[maxn*3];

void push_up(int i)
{
    segTree[i].sum=segTree[i<<1].sum+segTree[(i<<1)|1].sum;
}

void build(int i,int l,int r)
{
    segTree[i].l=l;
    segTree[i].r=r;
    if(l==r){
        segTree[i].sum=0;
        return ;
    }
    int mid=(l+r)/2;
    build(i<<1,l,mid);
    build((i<<1)+1,mid+1,r);
    push_up(i);
}

void update(int i,int val,int cnt)
{
    if(segTree[i].l==val&&segTree[i].r==val){
        segTree[i].sum+=cnt;
        return ;
    }
    int mid=(segTree[i].l+segTree[i].r)/2;
    if(val<=mid)update(i<<1,val,cnt);
    else update((i<<1)|1,val,cnt);
    push_up(i);
}

int query(int i,int l,int r)
{
    if(segTree[i].l==l&&segTree[i].r==r){
        return segTree[i].sum;
    }
    int mid=(segTree[i].l+segTree[i].r)/2;
    if(l>mid)return query((i<<1)|1,l,r);
    else if(r<=mid)return query(i<<1,l,r);
    else return query(i<<1,l,mid)+query((i<<1)|1,mid+1,r);
}

int main()
{
    //freopen("in.txt","r",stdin);
    int a[maxn],n,tot,ans;
    while(scanf("%d",&n)!=EOF){
        build(1,0,n-1);
        tot=0;
        for(int i=0;i<n;i++)scanf("%d",&a[i]);
        for(int i=0;i<n;i++){
            tot+=query(1,a[i],n-1);
            update(1,a[i],1);
        }
        ans=tot;
        for(int i=0;i<n;i++){
            tot+=n-1-2*a[i];
            if(ans>tot)ans=tot;
        }
        printf("%d\n",ans);
    }
    return 0;
}

3、归并排序

        将待排序序列R[0...n-1]看成是n个长度为1的有序序列，将相邻的有序表成对归并，得到n/2个长度为2的有序表；将这些有序序列再次归并，得到n/4个长度为4的有序序列；如此反复进行下去，最后得到一个长度为n的有序序列。

步骤：
     （1）“分解”——将序列每次折半划分。
     （2）“合并”——将划分后的序列段两两合并后排序。

       在每次合并过程中，都是对两个有序的序列段进行合并，然后排序。这两个有序序列段分别为 a[low, mid] 和 a[mid+1, high]。先将他们合并到一个局部的暂存数组tmp中，带合并完成后再将tmp复制回a中。为了方便描述，我们称 a[low, mid] 第一段，a[mid+1, high] 为第二段。每次从两个段中取出一个记录进行关键字的比较，将较小者放入tmp中。最后将各段中余下的部分直接复制到tmp中。经过这样的过程，tmp已经是一个有序的序列，再将其复制回a中，一次合并排序就完成了。

    while(i<=mid&&j<=r){  
        if(a[i]>a[j]){  //关键代码
            tmp[pos++]=a[j++];  
        }  
        else{  
            tmp[pos++]=a[i++];  
        }  
    }  
    while(i<=mid)tmp[pos++]=a[i++];  
    while(j<=r)tmp[pos++]=a[j++];  
    for(i=0;i<pos;i++)a[l+i]=tmp[i];  

注意上面的关键代码，正式我们要找的逆序对

对于

a[l],a[l+1],...a[i],....a[mid]   和 a[mid+1],...a[j]...a[r]

因为左右两个序列都有序，当a[i]>a[j]，即a[i]<a[i+1]<...<a[mid-1]<a[mid]，共mid-i+1个数

故tot+=mid-i+1

代码

#include <stdio.h>

const int maxn=5005;
int a[maxn],b[maxn],tmp[maxn],tot;
void merge(int l,int r,int mid)
{
    int pos=0;
    int i=l,j=mid+1;
    while(i<=mid&&j<=r){
        if(a[i]>a[j]){
            tot+=mid+1-i;
            tmp[pos++]=a[j++];
        }
        else{
            tmp[pos++]=a[i++];
        }
    }
    while(i<=mid)tmp[pos++]=a[i++];
    while(j<=r)tmp[pos++]=a[j++];
    for(i=0;i<pos;i++)a[l+i]=tmp[i];
}

void merge_sort(int l, int r)
{
    if(l==r)return;
    int mid=(l+r)/2;
    merge_sort(l,mid);
    merge_sort(mid+1,r);
    merge(l,r,mid);
}

int main()
{
    //freopen("in.txt","r",stdin);
    int n,ans;
    while(scanf("%d",&n)!=EOF){
        for(int i=0;i<n;i++){
            scanf("%d",&a[i]);
            b[i]=a[i];
        }
        tot=0;
        merge_sort(0,n-1);
        ans=tot;
        for(int i=0;i<n;i++){
            tot+=n-1-2*b[i];
            if(ans>tot)ans=tot;
        }
        printf("%d\n",ans);
    }
    return 0;
}