while (1) puts(“DP神题”)
初始思路
发现按顺序加数的话,这个数只会对其他人产生影响,并不会有大于他的数出现,最多可以出现min(i-1,m)个逆序对,所以一个显然的递推式 f [ i ] [ k ] = ∑ j = 0 m i n ( k , i − 1 ) f [ i − 1 ] [ k − j ] f[i][k]=\sum\limits_{j=0}^{min(k,i-1)}f[i-1][k-j] f[i][k]=j=0∑min(k,i−1)f[i−1][k−j],就是说加入i后,一共存在k个,它产生了j个逆序对,所以之前i-1个数共k-j个逆序对,直接转移 就行了
优化
发现上面的算法是
O
(
n
m
2
)
O(nm^2)
O(nm2)的复杂度,极限是
100
0
3
1000^3
10003,拿了90分,如何优化?
发现将上面式子换个形式变成
f
[
i
]
[
k
]
=
∑
j
=
m
a
x
(
0
,
k
−
i
+
1
)
k
f
[
i
−
1
]
[
j
]
f[i][k]=\sum\limits_{j=max(0,k-i+1)}^{k}f[i-1][j]
f[i][k]=j=max(0,k−i+1)∑kf[i−1][j],为啥呢?因为我们相当于是在一个长度为k的区间框一下就有k个逆序对了,所以往右移,左端点也移动,最小在0,之后就是k-i+1了,然后发现可以求个前缀和,因为每次都是从0开始的,然后就一路前缀和过去,把k-i+1之前的方案数要去掉
代码
//By AcerMo
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int mod=10000;
int n,m;
int f[1005][1005];
signed main()
{
scanf("%d%d",&n,&m);
f[1][0]=1;
for (int i=2;i<=n;i++)
{
int s=0;
for (int k=0;k<=m;k++)
{
s=(s+f[i-1][k])%mod;
f[i][k]=s;
if (k>=i-1) s=(s-f[i-1][k-i+1]+mod)%mod;
}
}
cout<<f[n][m]%mod;
return 0;
}