poj2299UltraQuickSort

http://poj.org/problem?id=2299

求逆序数

第一种做法 利用归并排序 在进行两个升序数组归并时 如果右边指针指向的那个数比左边小 那从左边那个位置到mid的数都比右边那个数大 所以比它大的数的个数就是 左边那个数到mid之间的数 依次类推

#include <stdio.h>
#include <string.h>
long long a[500001],num,x[500001],y[500001];
long long  msort(long long low,long long mid,long long high)
 {
     long long  i,j,n1 = mid-low+1,n2 = high-mid;
     for(i = 1 ; i <= n1 ; i++)
         x[i] = a[low+i-1];
     for(i = 1 ; i <= n2 ; i++)
         y[i] = a[mid+i];
     long long k = low;
     i = 1;
     j = 1;
     while(i<=n1&&j<=n2)
     {
         if(x[i]<=y[j])
             a[k] = x[i++];
         else
         {
             num+=mid-(low+i)+2;//加上i到mid之间的数
             a[k] = y[j++];
         }
         k++;
     }
     while(i<=n1)
         a[k++] = x[i++];
     while(j<=n2)
         a[k++] = y[j++];
 }
 void merge(long long low,long long high)
 {
     long long  mid;
     if(low<high)
     {
         mid = (low+high)/2;
         merge(low,mid);
         merge(mid+1,high);
         msort(low,mid ,high);
     }
 }
int main()
{
    int i,j,k,n,m;
    while(scanf("%d",&n)&&n)
    {
        num = 0;
        for(i = 1; i <= n ; i++)
        scanf("%d",&a[i]);
        merge(1,n);
        printf("%lld\n",num) ;
    }
    return 0;
}

 第二种 树状数组

树状数组时用来求和的 从小到大排好序之后 用一个数组记录原来位置的数现在在什么位置 然后调用更新函数 让原来那个位置标记为1 这样再调用getsum函数 就可以得到原来这个位置前面有多少比它小的数 依次加起来

#include <iostream>
#include <stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
typedef struct node
{
    int a,p;
}st;
st q[500001];
int re[500001],n,po[500001];
bool cmpp(st x, st y)
{
    return x.a<y.a;
}
void add(int i,int da)
{
    while(i<=n)
    {
       re[i]+=da;
       i += i&(-i);
    }
}
int getsum(int i)
{
    long long s = 0;
    while(i)
    {
        s+=re[i];
        i-=i&(-i);
    }
    return s;
}
int main()
{
    int i,j,m;
    long long s1;
    while(scanf("%d",&n)&&n)
    {
        s1 = 0;
        memset(re,0,sizeof(re));
        for(i = 1 ; i <= n ; i++)
        {
            scanf("%d",&q[i].a);
            q[i].p = i;
        }
        sort(q+1,q+n+1,cmpp);
        for(i = 1 ; i <= n ; i++)
            po[q[i].p] = i;
        for(i = 1 ; i <= n ; i++)
        {
            add(po[i],1);
            s1+=(i-getsum(po[i]));
        }
        printf("%lld\n",s1);
    }
    return 0;
}