【概率】CF749E Inversions After Shuffle

【题目】
原题地址
给定一个$n$的排列，现在等概率选取其中一个区间，将这个区间中的数字随机重排，求现在整个区间的期望逆序对个数。$n\le 10^5$

【解题思路】
我们可以从期望的定义入手，计算出选取每个区间的期望逆序对个数，再除以区间个数。进一步转化，我们需要求出每一个数对选取每个选取区间的贡献次数。这样子想已经可以做了，不过我们可以再做一步转化：计算选取每个区间时，会变化多少个逆序对的贡献和（记为$x$），这样我们用原来的逆序对个数加上$\frac{x}{n(n+1)}$就是答案了。

这样的话我们枚举选取区间的右端点$r$，那么我们要计算的就是当前数$a_r$与之前比它大的数的贡献，假设一个数$a_i>a_r$，那么他们的贡献就是$i$，于是这个东西我们可以用一个$\text{BIT}$来维护，这样子可以统计出下标$<r$所有大于$a_r$的数的贡献，这个贡献是负的。而我们只需要用$r(r-1)$与之相减就可以计算出正的贡献。（当然反过来算也可以）

这样子就做完了，复杂度是$O(n\log n)$的。
但是精度上需要注意一下，所以贡献的除以二我们放到最后，然后最后除以区间数时要写两个除法，不能先乘起来再除掉，不然就挂了啊。
然而这是个zz错误。

【参考代码】

#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
typedef double db;
const int N=1e5+10;
int n,a[N];
db ans;
ll sum,res;

int read()
{
	int ret=0;char c=getchar();
	while(!isdigit(c)) c=getchar();
	while(isdigit(c)) ret=ret*10+(c^48),c=getchar();
	return ret;
}

namespace BIT
{
	#define lowbit(x) (x&(-x))
	ll c[N];
	void clear(){memset(c,0,sizeof(c));}
	void update(int x,int d){for(;x<N;x+=lowbit(x))c[x]+=d;}
	ll query(int x){ll ret=0;for(;x;x-=lowbit(x))ret+=c[x];return ret;}
};
using namespace BIT;

int main()
{
#ifndef ONLINE_JUDGE
	freopen("CF749E.in","r",stdin);
	freopen("CF749E.out","w",stdout);
#endif
	n=read();
	for(int i=1;i<=n;++i) a[i]=read();
	for(int i=n;i;--i) ans+=query(a[i]),update(a[i],1);
	clear();
	for(int i=1;i<=n;++i)
	{
		ll t1=query(a[i]),t2=(ll)i*(i-1)/2-t1;
		sum+=(res+=t1-t2);update(a[i],i);
	}
	ans+=(db)sum/n/(n+1);
	printf("%0.15lf\n",ans);

	return 0;
}

【总结】
啊，这个转化的思路还是可以学习的，当然直接暴力去算贡献也是可以的233。