luogu P1774 最接近神的人_NOI导刊2010提高（02）

题目描述

破解了符文之语，小FF开启了通往地下的道路。当他走到最底层时，发现正前方有一扇巨石门，门上雕刻着一幅古代人进行某种活动的图案。而石门上方用古代文写着“神的殿堂”。小FF猜想里面应该就有王室的遗产了。但现在的问题是如何打开这扇门……

仔细研究后，他发现门上的图案大概是说：古代人认为只有智者才是最容易接近神明的。而最聪明的人往往通过一种仪式选拔出来。仪式大概是指，即将隐退的智者为他的候选人写下一串无序的数字，并让他们进行一种操作，即交换序列中相邻的两个元素。而用最少的交换次数使原序列变成不下降序列的人即是下一任智者。

小FF发现门上同样有着n个数字。于是他认为打开这扇门的秘诀就是找到让这个序列变成不下降序列所需要的最小次数。但小FF不会……只好又找到了你，并答应事成之后与你三七分……

输入输出格式

输入格式：

第一行为一个整数n，表示序列长度

第二行为n个整数，表示序列中每个元素。

输出格式：

一个整数ans，即最少操作次数。

输入输出样例

输入样例#1：

4
2 8 0 3

输出样例#1：

3
   样例说明：开始序列为2 8 0 3，目标序列为0 2 3 8，可进行三次操作的目标序列：
    1．Swap (8,0):2  0  8  3
    2．Swap (2,0):0  2  8  3
    3．Swap (8,3):0  2  3  8

说明

对于30%的数据1≤n≤10^4。

对于100%的数据1≤n≤5*10^5；

-maxlongint≤A[i]≤maxlongint。

归并排序求逆序对：

#include<iostream>
#include<cstdio>
#include<algorithm>
#define ll long long int

using namespace std;
const ll N=500010;

ll a[N];
ll ans[N];
ll answer;
ll n;
ll now;

ll read()
{
    ll x=0;char c=getchar();ll f=1;
    while(c<‘0‘||c>‘9‘){if(c==‘-‘)f=-1;c=getchar();}
    while(c>=‘0‘&&c<=‘9‘)x=x*10+c-‘0‘,c=getchar();
    return x*f;
}

void T_S(ll start,ll endd)
{
    if(start==endd)
        return ;
    int mid=(start+endd)/2;
    T_S(start,mid);
    T_S(mid+1,endd);
    ll i=start,j=mid+1;
    now=start;
    while(i<=mid&&j<=endd)
        if(a[i]<=a[j])
            ans[now++]=a[i++];
        else
            answer+=mid-i+1,
            ans[now++]=a[j++];
    while(i<=mid)
        ans[now++]=a[i++];
    while(j<=endd)
        ans[now++]=a[j++];
    for(int i=start;i<=endd;i++)
        a[i]=ans[i];
}

int main()
{
    n=read();
    for(int i=1;i<=n;i++)
        a[i]=read();
    T_S(1,n);
    printf("%lld",answer);
    return 0;
}

树状数组求逆序对（较慢）：

#include<bits/stdc++.h>

using namespace std;
const int maxn=1e6+7;
int n;
int sz[maxn],c[maxn];
struct node
{
    int dat,pos;
}nd[maxn];

bool cmp(node a,node b)
{
    return a.dat<b.dat;
}
inline int lowbit(int x)
{
    return x&(-x);
}
int getsum(int x)
{
    int sum=0;
    for(;x>0;x-=lowbit(x))sum+=c[x];
    return sum;
}
void updat(int x)
{
    for(;x<=n;x+=lowbit(x))c[x]++;
}
int main()
{
    cin>>n;
    for(int i=1;i<=n;i++)
    {
        scanf("%d",sz+i);
        nd[i].dat=sz[i];nd[i].pos=i;
    }
    sort(nd+1,nd+1+n,cmp);
    int js=1;
    sz[nd[1].pos]=1;
    for(int i=2;i<=n;i++)
    {
        if(nd[i].dat==nd[i-1].dat)sz[nd[i].pos]=js;
        else sz[nd[i].pos]=++js;
    }
    long long ans=0;
    for(int i=1;i<=n;i++)
    {
        ans+=i-getsum(sz[i])-1;//<=sz[i]
        updat(sz[i]);
    }
    cout<<ans;
    return 0;
}