2018北大暑校acm算法训练课程求排列的逆序数树状数组/分治

最新推荐文章于 2021-07-17 23:23:43 发布

空想昼夜

最新推荐文章于 2021-07-17 23:23:43 发布

阅读量286

点赞数

分类专栏： pku_acm

本文链接：https://blog.csdn.net/qq_38842456/article/details/81291277

版权

pku_acm 专栏收录该内容

17 篇文章 1 订阅

订阅专栏

总时间限制: 1000ms 内存限制: 65536kB
描述
在Internet上的搜索引擎经常需要对信息进行比较，比如可以通过某个人对一些事物的排名来估计他（或她）对各种不同信息的兴趣，从而实现个性化的服务。

对于不同的排名结果可以用逆序来评价它们之间的差异。考虑1,2,…,n的排列i1，i2，…，in，如果其中存在j,k，满足 j < k 且 ij > ik，那么就称(ij,ik)是这个排列的一个逆序。

一个排列含有逆序的个数称为这个排列的逆序数。例如排列 263451 含有8个逆序(2,1),(6,3),(6,4),(6,5),(6,1),(3,1),(4,1),(5,1)，因此该排列的逆序数就是8。显然，由1,2,…,n 构成的所有n!个排列中，最小的逆序数是0，对应的排列就是1,2,…,n；最大的逆序数是n(n-1)/2，对应的排列就是n,(n-1),…,2,1。逆序数越大的排列与原始排列的差异度就越大。

现给定1,2,…,n的一个排列，求它的逆序数。

输入
第一行是一个整数n，表示该排列有n个数（n <= 100000)。
第二行是n个不同的正整数，之间以空格隔开，表示该排列。
输出
输出该排列的逆序数。
样例输入
6
2 6 3 4 5 1
样例输出
8
提示
1. 利用二分归并排序算法（分治）；
2. 注意结果可能超过int的范围，需要用long long存储。

当初这个题是学树状数组的板子题,所以就直接用树状数组做了.

#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<stack>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<sstream>
#include<iostream>
#include<algorithm>
const int MAX_N=101010;
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int INF = 0x3f3f3f3f;
#define eps 1e-7
//#define Kongxiangzhouye
#define _for(i,a,b) for( int i=(a); i<(b); ++i)
#define _rep(i,a,b) for( int i=(a); i<=(b); ++i)
inline int read(){
        int a=0;
        int fh=1;
        char c=getchar();
        while(!(c<='9'&&c>='0')){
            if(c=='-')
                fh*=-1;
            c=getchar();
        }
        while(c<='9'&&c>='0'){
            a=a*10+c-'0';
            c=getchar();
        }
        return a*fh;
}
int tree[MAX_N];
int c[MAX_N];
int n;
struct Node{
    int v;
    int order;
}a[MAX_N];
struct cmpFunctor{
    inline bool operator ()(const Node& t1,const Node& t2){
        return t1.v<t2.v;
    }
};

int lowbit(int k){return k&(-k);}

void update(int t,int value){ 
    for(int i=t;i<=n;i+=lowbit(i))
        c[i]+=value;  
}

int getsum(int t){  
    int sum=0;
    for(int i=t;i>=1;i-=lowbit(i))
        sum += c[i];
    return sum;
}

int main(){
    while (scanf("%d",&n)!=EOF&&n){
        _rep(i,1,n){
            a[i].v=read();
            a[i].order = i;
        }
        sort(a+1,a+n+1,cmpFunctor());
        memset(c,0,sizeof(c));
        _rep(i,1,n)
            tree[a[i].order]=i;
        ll ans=0;
        _rep(i,1,n){
            update(tree[i],1);
            ans+=i-getsum(tree[i]);
        }
        printf("%lld\n", ans);
    }
    return 0;
}