俄罗斯套娃【DP】【滚动数组】

>Description

---

>Input

>Output

---

>Sample Input
10 1000

>Sample Output
3628800

---

>解题思路
感谢lyf巨佬帮我改码风然后从全wa到a，说明码风真的很重要TT
DP题。

设$f[i][j]$为前$i$个数（长度为n的序列为1~n），有$j$对逆序对的方案数，每次到一个位置i就枚举第i个数放的位置，从而添加逆序对，第i个数可以添加的逆序对为0 ~ i-1
因此得出状态转移方程（l枚举的是当前数添加的逆序对）：

然而这样只能拿到60分。
从状态转移方程可以看出，$f[i][j]$的转移只考虑$f[i-1][j...j-i+1]$这一段的和，所以我们可以累积前缀和，转移时只用O(1)就可以完成，然后再加一个滚动数组就OK了。

---

>代码

#include <iostream>
#include <cstdio>
#include <cstring>
#define ll long long
using namespace std;

const ll p = 10000000007;
int n, k, t;
ll f[5][5005], sum[5][5005];

int main()
{
	scanf ("%d%d", &n, &k);
	f[0][0] = 1; //dp
	for(int i = 0; i <= k; i++)
	  sum[0][i] = 1; //前缀和
	for (int i = 2; i <= n; i++)
	{
		t ^= 1; //滚动
		for (int j = 0; j <= k; j++)
		{
	  		if (j - i >= 0) f[t][j] = (p + sum[t ^ 1][j] - sum[t ^ 1][j - i])% p;
	  		                       //对于差的模法，结果可能为负数，所以要加一个p在模防止出错TT
	  		else f[t][j] = sum[t ^ 1][j];
	  		if (j > 0) sum[t][j] = (sum[t][j - 1] + f[t][j]) % p;
	  		else sum[t][j] = f[t][j];
	  		//要判断边界TT
		}
	}
	printf ("%lld", (sum[t][k] + p) % p);
	return 0;
}