HDU 3664（DP）

解题思路：用dp[i][j]表示第n个排列E值为j。a[k]表示k位置处的值。对于n+1这个数，放在排列中有三种情况，（1）放在末尾E值不变。（2）放在a[k] > k处，E值也不变，有j种情况。（3）放在a[k] <= k处，E值减少1，有i - j种情况。

故dp[i][j] = dp[i - 1] + j * dp[i - 1][j] + (i - j) * dp[i - 1][ j - 1]。 

#include<cstdio>const int mod = 1e9 + 7;long long dp[1010][1010];void init(){    dp[1][0] = 1;    for(int i = 1; i <= 1000; i++){        dp[i][0] = 1;        for(int j = 1; j <= i; j++)            dp[i][j] = (dp[i - 1][j]  + j * dp[i - 1][j]  + (i - j) * dp[i - 1][j - 1]) % mod;    }}int main(){    init();   int n, k;   while(~scanf("%d%d", &n, &k)){      printf("%I64d\n", dp[n][k]);   }}