hdoj 3664 Permutation Counting

类型：DP

题目：http://acm.hdu.edu.cn/showproblem.php?pid=3664

来源：2010 Asia Regional Harbin

思路：状态dp(n, k)表示n个数有k个ai > i的排列数。考虑第n个数，设当前排列为a1 a2 ...ai...an -1

（1）放在末尾，不影响当前排列的结果

 （2） 和某一个ai交换，如果ai <= i，则结果+1，如果ai > i，则结果不变

这样构造可以保证不重不漏

可以得到：dp(n, k) = (k + 1)dp(n - 1, k) + ((n - 1)- (k - 1))dp(n - 1, k - 1)

// hdoj 3664 Permutation Counting
#include <iostream>
#include <sstream>
#include <string>
#include <queue>
#include <stack>
#include <vector>
#include <list>
#include <map>
#include <algorithm>
#include <numeric>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <ctime>
#include <set>
#include <deque>
#include <bitset>
#include <functional>
#include <utility>
#include <iomanip>
#include <cctype>
using namespace std;

#define FOR(i,a,b) for(i = (a); i < (b); ++i)
#define FORE(i,a,b) for(i = (a); i <= (b); ++i)
#define FORD(i,a,b) for(i = (a); i > (b); --i)
#define FORDE(i,a,b) for(i = (a); i >= (b); --i)
#define max(a,b) ((a) > (b)) ? (a) : (b)
#define min(a,b) ((a) < (b)) ? (a) : (b)
#define CLR(a,b) memset(a,b,sizeof(a))
#define PB(x) push_back(x)

typedef long long LL;
typedef vector<int> VI;

const int MAXN = 1010;
const int MAXM = 0;
const int hash_size = 25000002;
const int inf = 1000000007;

long long dp[MAXN][MAXN];

int n, k;

void creat() {
    int i, j;

    dp[0][0] = 1, dp[1][0] = 1;
    FORE(i, 2, 1000)
        FORE(j, 0, i) {
            dp[i][j] = ((j + 1) * dp[i - 1][j] % inf + (i - j) * dp[i - 1][j - 1] % inf) % inf;
        }
}

int main()
{
    creat();
    while(scanf("%d %d", &n, &k) != EOF) {
        cout<<dp[n][k]<<endl;
    }
    return 0;
}