首先为了使状态更方便转移,我们显然要从大到小或从小到大往当前序列插入进行转移。
我们尝试在一个已经有j个逆序对和前i-1个数字组成的数列中插入,我们会发现将数字i插入可以很轻松的算出新的逆序对贡献。
作为这个数列中最大的数字,显然呆在它后面的都是逆序对咯。
于是得到动归方程:
f[i][j]=Sigma(f[i-1][j-i+k]) k∈[1,i]
首先发现一个性质,用的总是i-1状态,那么我们滚动一下就轻松剩下了一大堆空间。(直接去掉了一维)
然后重点来了!原来这个DP是O(n^3)的,但是这个递推式居然是用连续区间和,那么加个前缀和维护不就秒杀了。
------------------------------------------------------------------------------------------------------------------
Code:
#include<iostream>
#include<cstdio>
#include<cstring>
#define For(i,l,r) for(int (i)=(l);(i)<=(r);(i)++)
using namespace std;
const int N=1010,MOD=10000;
int n,k,f[N],sum[N];
int main(){
scanf("%d%d",&n,&k);
For(i,0,k) sum[i]=1;
For(i,2,n){
f[0]=1;
For(j,1,k){
if(j-i>=0){
f[j]=(sum[j]-sum[j-i])%MOD;
while(f[j]<0) f[j]+=MOD;
}else f[j]=sum[j]%MOD;
}
For(j,1,k){
sum[j]=sum[j-1]+f[j];
sum[j]%=MOD;
}
}printf("%d\n",f[k]);
return 0;
}
要注意一下算术溢出,我也是对拍才发现的QAQ