HDU 3664 Permutation Counting(DP)

题目链接

弱爆啦，组合弱爆了，反正是没想出来怎么搞这个题，其实这个公式不难推啊，反正就是没推出来。今天队内赛，实在是没办法了，暴力写了个DFS，先把10以内的打出表来，发现类似杨辉三角的一个表，推不出公式，只能找规律了。也推公式，也找规律，中间还走会了神，发现borad上过的人N多了，有些着急，这样应该不难吧。。。又推了会，还是没想出来，找规律吧，估摸着应该是和上两项有关系，自己写了小程序测试一下几个数据和scf讨论了下，貌似真的是找出规律了。。。然后时间不多了，好在代码很短，马上快结束了，乱写了，最后在各种乱搞+思考之下，和暴力的写的数据对上了，中间错了2次，好在在4:53的时候终于AC了。。。。

公式的意义：p[i][j] = (i-j)*p[i-1][j-1] + (j+1)*p[i-1][j];对最后一个数进行讨论，最后一个数字，如果和a[i] > i的位置交换，最后一个数一定比那个位置大，所以之前就必须有j个a[i] > i的数，所以方案数为(j+1)*p[i-1][j] +1代表最后一个数的位置，如果最后一个数和a[i] < i的位置交换共有(i-1) - (j-1)个 ，则方案数为(i-j)*p[i-1][j-1]。

#include <cstdio>
#include <cstring>
using namespace std;
#define MOD 1000000007
__int64 p[1001][1001];
int main()
{
    int i,j,n,k;
    for(i = 1;i <= 1000;i ++)
    {
        p[i][0] = 1;
        p[i][i-1] = 1;
    }
    for(i = 3;i <= 1000;i ++)
    {
        for(j = 1;j <= (i+1)/2-1;j ++)
        {
             p[i][i-1-j] = p[i][j] = (((i-j)*p[i-1][j-1])%MOD+((j+1)*p[i-1][j])%MOD)%MOD;
        }
    }
    while(scanf("%d%d",&n,&k)!=EOF)
    {
        printf("%I64d\n",p[n][k]);
    }
    return 0;
}