HDU5792 World is Exploding(树状数组)

#include<cstdio>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<string>
#include<algorithm>
#include<map>
#include<queue>
#include<vector>
#include<cmath>
#include<utility>
using namespace std;
typedef long long LL;
const int N = 60008, INF = 0x3F3F3F3F;
#define MS(a, num) memset(a, num, sizeof(a))
#define PB(A) push_back(A)
#define FOR(i, n) for(int i = 0; i < n; i++)
int C[N];
int n;
int le[N], ri[N], le1[N], ri1[N];
int val[N];
int b[N], tp[N];
inline int lowbit(int x){
    return x&-x;
}
inline void add(int x, int val){
    for(int i=x;i<=n;i+=lowbit(i)){
        C[i] += val;
    }
}
inline int sum(int x){
    int ret = 0;
    for(int i=x;i>0;i-=lowbit(i)){
        ret+=C[i];
    }
    return ret;
}

int main(){
    while(~scanf("%d", &n)){
        for(int i = 1; i <= n; i++){
            scanf("%d", &val[i]);
            tp[i] = val[i];
        }

        sort(val + 1, val + n + 1);
        for(int i = 1; i <= n; i++){
            b[i]  = lower_bound(val + 1, val + n + 1, tp[i]) - val;
        }
        MS(C , 0);
        for(int i = 1; i <= n; i++){
            le[i] = sum(b[i] - 1);
            le1[i] = sum(n) - sum(b[i]);
            add(b[i], 1);
        }

        MS(C, 0);
        for(int i = n; i >= 1; i--){
            ri[i] = sum(b[i] - 1);
            ri1[i] = sum(n) - sum(b[i]);
            add(b[i], 1);
        }
        LL ans= 0;
        LL sum = 0;
        for(int i = n; i >= 1; i--){
            sum += ri1[i];
        }
        ans = sum;
        sum = 0;
        for(int i = n; i >= 1; i--){
            sum += ri[i];
        }
        ans *= sum;
        sum = 0;
        for(int i = 1; i <= n; i++){
            sum += (LL)ri[i] * (LL)ri1[i];
            sum += (LL)ri1[i] * (LL)le1[i];
            sum += (LL)le[i] * (LL)le1[i];
            sum += (LL)le[i] * (LL)ri[i];
        }
        ans -= sum;
        printf("%I64d\n", ans);
    }
    return 0;
}