一道比较不错的组合思想题

$方法一(O(n^3))$

$将原数组进行排序，二重循环枚举两个数的值，但是发现很难计算方案数$
$所以就再加上一个这个值所在的位置，记为k$
$可以将原数组进行排序，较小的数为x，较大的数为y，比x小的数为sm，比x大的数为bg，在两个中间的为mid$
$很显然的sm的数一定是在k之前，bg的数一定是在k之后，mid是在x之后，y之前$
$$\begin{gathered} 由于已经知道k的位置了，那么a前有几个sm，b后有几个bg都可以直接解决了，这个方案数即为 \\ {C}_{sm}^{k-1}\ast C_{bg}^{n-k}\ast |a[i]-a[j]|\ast C_{mid}^{2\ast k-2-sm}\ast 2(有点繁琐，自己推) \end{gathered}$$
这样就可以直接了

#include<bits/stdc++.h>
using namespace std;
const int N=5e3+10;
#define Mod 998244353
#define ll long long
ll power(ll x,ll y){
	if(!y)return 1;
	ll sum=power(x,y>>1);
	if(y&1)return sum*sum%Mod*x%Mod;
	return sum*sum%Mod;
}
int n,a[N];ll fac[N],inv[N],ans;
ll C(int n,int k){
    if(n<k)return 0;
	return fac[n]*inv[k]%Mod*inv[n-k]%Mod;
}
ll P(int n,int k){
	return fac[n]*inv[n-k]%Mod;
}
int main()
{
	scanf("%d",&n);fac[0]=1;
	for(int i=1;i<=2*n;i++)
	  scanf("%d",&a[i]);
	for(int i=1;i<=2*n;i++)
	  fac[i]=fac[i-1]*i%Mod;
	inv[2*n]=power(fac[2*n],Mod-2);
	for(int i=2*n-1;i+1;i--)
	  inv[i]=inv[i+1]*(i+1)%Mod;
	sort(a+1,a+2*n+1);
	for(int i=1;i<=2*n;i++)
	  for(int j=i+1;j<=2*n;j++){
	  	int sm=i-1,bg=2*n-j,mid=j-i-1;
	  	for(int k=n-bg;k<=sm+1;k++){
	  	  ans=(ans+C(sm,k-1)*C(bg,n-k)%Mod*abs(a[i]-a[j])%Mod
			*2%Mod*C(mid,2*k-2-sm)%Mod)%Mod;
		}
	  }
	printf("%lld",ans);
	return 0;
}

$方法二(O(n^2))$

$可以发现这个题目又有一个性质：如果a数组确定了，那么b数组自然而然地也确定了$
$可以将|ai-bi|分成(ai-bi)和(bi-ai)，只需要考虑这个式子其中一个值应该是正的还是负的$
$所以可以枚举i和j(即为所在的位置)的值$
$示意图如下$
$可以发现当j\leq i-j时，即2*j<=i时，显然a[i]为正数，否则a[i]为负数\$
$方案数求解也非常好写$

#include<bits/stdc++.h>
using namespace std;
const int N=5e3+10;
#define Mod 998244353
#define ll long long
ll power(ll x,ll y){
	if(!y)return 1;
	ll sum=power(x,y>>1);
	if(y&1)return sum*sum%Mod*x%Mod;
	return sum*sum%Mod;
}
int n;ll fac[N],inv[N],ans,a[N];
ll C(int n,int k){
	if(k>n)return 0; 
	return fac[n]*inv[k]%Mod*inv[n-k]%Mod;
}
signed main()
{
	scanf("%d",&n);fac[0]=1;
	for(int i=1;i<=2*n;i++)
	  scanf("%d",&a[i]);
	for(int i=1;i<=2*n;i++)
	  fac[i]=fac[i-1]*i%Mod;
	inv[2*n]=power(fac[2*n],Mod-2);
	for(int i=2*n-1;i+1;i--)
	  inv[i]=inv[i+1]*(i+1)%Mod;
	sort(a+1,a+2*n+1);
	for(int i=1,p;i<=2*n;i++)
	  for(int j=1;j<=min(i,n);j++){
	  	if(i>=2*j)p=1;else p=-1;
	  	ans=((ans+2ll*p*a[i]%Mod*C(i-1,j-1)
		%Mod*C(2*n-i,n-j)%Mod)%Mod+Mod)%Mod;
	  }
	printf("%lld",ans);
	return 0;
}

$方法三(O(nlo{g}^{2}n))$

不会…