解题思路:用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]);
}
}