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]);
   }
}