雅礼集训2019 day7

Inverse

$1\le n\le 500,0\le k\le 50$, P是一个1到n的排列.
概率$dp$好题。
思路：
定义$f_{i,j,k}$表示$k$轮变换之后$a_{p_i}>a_{p_j}$的概率。
然后考虑对当前轮翻转的区间$[l,r]$分类转移。

1. $r<iorl>jori<l<r<j$，则$f_{i,j,k-1}\to f_{i,j,k}$
2. $l\le i\le r<j$，则$f_{l+r-i,j,k-1}\to f_{i,j,k}$
3. $i<l\le j\le r$，则$f_{i,l+r-j,k-1}\to f_{i,j,k}$
4. $l\le i<j\le r$，则$1-f_{l+r-j,l+r-i}\to f_{i,j,k}$

利用二阶前缀和来优化转移可以做到$O(n^{2}k)$
代码：

#include<bits/stdc++.h>
#define ri register int
using namespace std;
const int rlen=1<<18|1;
inline char gc(){
   
	static char buf[rlen],*ib,*ob;
	(ib==ob)&&(ob=(ib=buf)+fread(buf,1,rlen,stdin));
	return ib==ob?-1:*ib++;
}
inline int read(){
   
	int ans=0;
	char ch=gc();
	while(!isdigit(ch))ch=gc();
	while(isdigit(ch))ans=((ans<<2)+ans<<1)+(ch^48),ch=gc();
	return ans;
}
typedef long long ll;
const int mod=1e9+7;
inline int add(const int&a,const int&b){
   return a+b>=mod?a+b-mod:a+b;}
inline int dec(const int&a,const int&b){
   return a>=b?a-b:a-b+mod;}
inline int mul(