求逆序数之分治排序

winston

求逆序数
时间限制：2000 ms | 内存限制：65535 KB
难度：5
描述
在一个排列中，如果一对数的前后位置与大小顺序相反，即前面的数大于后面的数，那么它们就称为一个逆序。一个排列中逆序的总数就称为这个排列的逆序数。

现在，给你一个N个元素的序列，请你判断出它的逆序数是多少。

比如 1 3 2 的逆序数就是1。

输入
第一行输入一个整数T表示测试数据的组数（1<=T<=5)
每组测试数据的每一行是一个整数N表示数列中共有N个元素（2〈=N〈=1000000）
随后的一行共有N个整数Ai(0<=Ai<1000000000)，表示数列中的所有元素。

数据保证在多组测试数据中，多于10万个数的测试数据最多只有一组。
输出
输出该数列的逆序数
样例输入
2
2
1 1
3
1 3 2样例输出
0
1
思路：先分治后归并
AC代码：#include<iostream>

2. #include<string.h>
3. #include<string>
4. #define N 1000001
5. using namespace std;
6. long  a[N],b[N];
7. long long num;//尼玛的这里让我错了好几次。。。
8. void merge(int p,int q,int r,int s)
9. {   int a1=p,k=p,c=r;
10.     while(a1<=q&&c<=s)
11.     {  if(a[a1]<=a[c]) b[k++]=a[a1++];
12.        else{b[k++]=a[c++];num+=q-a1+1;}
13.     }
14.     if(a1>q)
15.     for(;c<=s;++c)
16.     b[k++]=a[c];
17.     if(c>s)
18.     for(;a1<=q;++a1)
19.     b[k++]=a[a1];
20.     for(;p<k;++p)
21.     a[p]=b[p];
22. }
23. void merge_sort(int p,int r)
24. { if(p<r)
25.   {  int q=(p+r)/2;
26.       merge_sort(p,q);
27.       merge_sort(q+1,r);
28.       merge(p,q,q+1,r);
29.   }
30. }
31. int  main()
32. {  int Case;
33.    cin>>Case;
34.    while(Case--)
35.    { memset(a,0,sizeof(a));
36.      memset(b,0,sizeof(b));
37.     num=0;
38.      int n;
39.       cin>>n;
40.       for(int i=0;i<n;++i)
41.       cin>>a[i];
42.       merge_sort(0,n-1);
43.       cout<<num<<endl;
44.    }return 0;
45. }

复制代码

作者：smallacmer 发表于2012-2-17 8:42:53 原文链接