51nod 1020 1到n的全排列逆序对为k的方案数

题目链接：

http://www.51nod.com/Challenge/Problem.html#problemId=1020

题意：

在一个排列中，如果一对数的前后位置与大小顺序相反，即前面的数大于后面的数，那么它们就称为一个逆序。一个排列中逆序的总数就称为这个排列的逆序数。

如2 4 3 1中，2 1，4 3，4 1，3 1是逆序，逆序数是4。

1-n的全排列中，逆序数最小为0（正序），最大为n*(n-1) / 2（倒序）

给出2个数n和k，求1-n的全排列中，逆序数为k的排列有多少种？

例如：n = 4 k = 3。

1 2 3 4的排列中逆序为3的共有6个，分别是：

1 4 3 2

2 3 4 1

2 4 1 3

3 1 4 2

3 2 1 4

4 1 2 3

由于逆序排列的数量非常大，因此只需计算并输出该数 Mod 10^9 + 7的结果就可以了。

输入

第1行：一个数T，表示后面用作输入测试的数的数量。（1 <= T <= 10000)
第2 - T + 1行：每行2个数n，k。中间用空格分隔。（2 <= n <= 1000, 0 <= k <= 20000)

输出

共T行，对应逆序排列的数量 Mod (10^9 + 7)

输入样例

1
4 3

输出样例

6

题解：

用dp[i][j]表示i个不同的数字全排列，逆序对为j的方案数

dp[i][j]=sigma dp[i-1][k]     max(0,j-(i-1))<=k<=j

状态转移是根据最末尾的数字的贡献度进行枚举的，如果最末尾数字是s，那么s贡献的逆序对数是s-1.

由于最末尾数字是1的时候贡献度最大，是i-1.所以 边界条件是 max(0,j-(i-1)).

根据状态转移方程可以看出，每个状态等于上一行的某个区间和，区间和可以看成一个滑动窗口，可以用一个变量来维护。

预处理的复杂度是O(nk) ,每次查询复杂度O(1)

题目给了128M的内存，如果数组用long long 大概需要160M的内存，int大概是80M.所以用int数组。

代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=1e9+7;
int f[1001][20001];
void init(){
    f[1][0]=1;
    for(int i=2;i<=1000;i++){
        ll sum=0;
        for(int j=0;j<=min(20000,i*(i-1)/2);j++){
             sum=(sum+f[i-1][j])%mod;
             if(j>=i) sum=(sum-f[i-1][j-i]+mod)%mod;
             f[i][j]=sum;
        }
    }
}
int main(){
    init();
    int t;
    scanf("%d",&t);
    while(t--){
        int n,k;
        scanf("%d%d",&n,&k);
        printf("%d\n",f[n][k]);
    }
    return 0;
}