种花(derangement)

题意

---

定义一个排列a，$f(a)={\sum }_{i<j,且a_i>a_j}(j-i)(a_i-a_j)$
求对于所有n的排列a，满足$a_i̸=p_i,f(a)$的和，p为排列
$n\le 5000$

---

题解

---

答案怎么算

很显然，每两个位置的贡献独立
设$f_0(x),f_1(x),f_2(x)$分别表示n的排列，前x，x-1，x-2个位置有$a_i̸=i$的限制，的个数
简单地说，$f_0$是正宗错排，$f_1$是少一个限制的错排，$f_2$是少两个限制的错排
既然这样，对于每两个位置统计对答案的贡献，对于j<i两个位置有4种摆放方式：

1. { i , j } ̸ = { p i , p j } \{i,j\}\not=\{p_i,p_j\} {i,j}̸​={pi​,pj​}，这两位的答案是$$\begin{gathered} (\sum _{1\le x\le n}\sum _{1\le y\le n}|x-y|+\sum _{1\le x\le n}|x-p_i|-\sum _{1\le y\le n}|p_j-y|+abs(p_i-p_j))(i-j)， \\ (不妨使大数放在左边小数放在右边) \end{gathered}$$化简之后就是$$(\frac{1}{6}n(n-1)(n+1)-n^2-n-p_i^2-p_j^2+(p_i+p_j)(n+1)+|p_i-p_j|)(i-j)$$其他的位置要乘上$f_2(n-2)$的方案数
2. $a_i=p_j,a_j=p_i$，答案为$\max (0,p_i-p_j)(i-j)$，要乘上$f_0(n-2)$
3. $a_i=p_j,a_j̸=p_i$或$a_j=p_i,a_i̸=p_j$，总答案为$$(\sum _{1\le x<p_i}{p}_i-x+\sum _{p_j<y\le n}y-p_j-2\max (0,p_i-p_j))$$化简之后就是$$\frac{1}{2}(p_i-1)p_i+\frac{1}{2}(n-p_j+1)(n-p_j)-2\max (0,p_i-p_j)$$应该乘上$f_1(n-2)$

f怎么求？

1. $f_0(n)$正宗·错排数，众所周知$f_0(n)=(n-1)[f_0(n-1)+f_0(n-2)]$
2. $f_1(n)$，若n放在第n为，就是n-1的错排；若n不放在第n位，故多了一个限制，是n的错排，故$f_1(n)=f_0(n)+f_0(n-1)$
3. $f_2(n)$，和1差不多，若第n-1放n-1，第n放n，就是n-2的错排；若n放n且第n-1不放n-1，或者是第n-1放n-1且第n不放n，就相当于n-1的错排；若第n不放n，且第n-1不放n-1，就相当于n的错排，故$f_2(n)=f_0(n)+2f_0(n-1)+f_0(n-2)$

---

code

#include<cstring>
#include<cstdio>
#include<algorithm>
#include<cmath>
#define fo(i,a,b) for(int i(a),_E_(b);i<=_E_;++i)
#define REP(i,a,b) for(int i(a),_E_(b);i<_E_;++i)
#define fd(i,a,b) for(int i(a),_E_(b);i>=_E_;--i)
#define pu putchar
#define ge getchar
#define mp(a,b) make_pair(a,b)
#define fs first
#define sc second
#define ll long long
#define ul unsigned long long
#define ui unsigned int
#define db double
#define ld long double
#define N 100000
#define mo 1000000009
#define ny 500000005
#define n_ 833333341

using namespace std;

void read(int &n){
	n=0;int s=1;char c;for(;(c=ge())>'9'||c<'0';)if(c=='-')s=-1;
	for(;c>='0'&&c<='9';c=ge())n=(n<<1)+(n<<3)+c-48;n*=s;
}
void write(int n){
	if(!n)pu(48);int t=0,b[20];
	for(;n;n/=10)b[++t]=n%10;while(t)pu(b[t--]+48);
}

int p[N],n;
ll f[N],ans;

ll f0(int x){return f[x];}
ll f1(int x){return ((x?f[x-1]:0)+f[x])%mo;}
ll f2(int x){return ((x>1?f[x-2]:0)+(x?2*f[x-1]:0)+f[x])%mo;}

int main(){
	scanf("%d",&n);
	if(n==1){
		printf("0");return 0;
	}
	fo(i,1,n)scanf("%d",p+i);
	f[0]=1;fo(i,2,n)f[i]=1ll*(i-1)*(f[i-1]+f[i-2])%mo;
	fo(i,1,n)REP(j,1,i){
		ll x=p[i],y=p[j];
		ans=(ans+(x>y)*(x-y)*(i-j)*f0(n-2))%mo;
		ans=(ans+((n-y+1)*(n-y)*ny+x*(x-1)*ny-2*(x>y)*(x-y)+mo)%mo*(i-j)%mo*f1(n-2))%mo;
		ans=(ans+(1ll*n*(n+1)*(n-1)%mo*n_-n*n-n-x*x-y*y+(x+y)*(n+1)+mo+1ll*abs(x-y))%mo*(i-j)%mo*f2(n-2))%mo;
	}printf("%lld",ans);
	return 0;
}