codeforces 285E Positions in Permutations（DP+容斥，hard）

题意：
定义1～N构成的排列 p(1)p(2)…p(N)
满足 abs(p(i) - i) = 1 的称之为 good position
统计有K个good position的排列数

思路：
每个位置要么是 good position，要么不是。。
我们可以先考虑占坑，有N个坑，我们先占r个(0<=r<=N)作为good position，其余的数任意排列，这样我们就知道 good position 数大于等于 r 的排列数了。
假设a[i]是 good position值大于等于i的排列数。
这样，a[i] - a[i+1] 是否就是值恰好为i的呢？
当然不是。。
因为a[i]多数了，a[i+1]被计算了C(i+1, i)次。。
所以需要容斥一下。。
题解代码中直接两重循环减重复的部分。。
然后我自己yy了一下容斥。。原理不是很明白。。
//==========//
好了，dp部分怎么来呢。。
当考虑 i 是否为 good position 这两种状态时， 需要知道 i-1 和 i+1 是否已经被使用过了
如果定义 前i个位置恰好有j个good position 的计数
那么 i 位置只需要知道 i-1是否已被用了。。
这个dp出来是占坑的方法数，然后a[k] = dp[k] * (n-k)!

LL dp[Maxn+5][Maxn+5][2][2], n, k;
LL C[Maxn+5][Maxn+5], ans[Maxn+5];

inline void mod_add(LL &x, LL y) {
    x += y;
    if (x > Mod) x -= Mod;
}

int main() {
#ifndef ONLINE_JUDGE
    freopen("input.in", "r", stdin);
#endif
    cin >> n >> k;

    C[0][0] = 1;
    rep(i, 1, n) {
        C[i][0] = 1;
        rep(j, 1, i)
            mod_add( C[i][j], C[i-1][j] ), mod_add( C[i][j], C[i-1][j-1] );
    }

    // 前i个位置，有j个good position，i, i+1 是否被使用
    dp[0][0][1][0] = 1;
    rep(i, 1, n) rep(j, 0, i-1) rep(now, 0, 1) rep(nxt, 0, 1) {
        if (!now)  mod_add(dp[i][j+1][nxt][0], dp[i-1][j][now][nxt]);
        if (i < n) mod_add(dp[i][j+1][nxt][1], dp[i-1][j][now][nxt]);
        mod_add(dp[i][j][nxt][0], dp[i-1][j][now][nxt]);
    }

    // 对good position之外的位置全排列
    // 方程第二维语义变成 有大于等于j个good position
    rep(i, 0, n) {
        rep(now, 0, 1) rep(nxt, 0, 1) mod_add(ans[i], dp[n][i][now][nxt]);
        rep(j, 2, n-i) ans[i] = ans[i]*j%Mod;
    }

    // 容斥
    int flag = -1;
    rep(i, k+1, n) {
        ans[k] = (ans[k] + flag * C[i][k] * ans[i]) % Mod;
        if (ans[k] < 0) ans[k] += Mod;
        flag *= -1;
    }

    cout << ans[k];

    return 0;
}