「日常訓練」 神、上帝以及老天爺 （HDU 2048）

題意

數論中的錯排問題。記錯排爲Dn$D_n$ ，求Dnn!$\frac{D_n}{n!}$ 。

分析

顯然D1=0,D2=1$D_1=0,D_2=1$ 。當n≥3$n\ge 3$ 時，不妨設n排在了第k位，其中k≠n$k\ne n$ ，也就是1≤k≤n−1$1\le k\le n-1$ 。那麼我們現在考慮k的情況。

• 當k排在第n位時，n放k、k放n確定，因而除了n和k以外還有n-2個數，其錯排數爲Dn−2$D_{n-2}$ 。
• 當k不排在第n位時，由於n在第k位，那麼將第n位重新考慮成一個新的“第k位”，那麼注意到“n在k”是確定的。剩下的排布方式同原來一致，只是k等價於n，因此：其錯排數爲Dn−1$D_{n-1}$ 。
  我們認真思考一下這個“重新考慮”是什麼意思。 n在k，這是一個固定的條件，故這裏不作考慮；而其他的元素必定不在各自的位置，k也是：它必定不在n！因而，這裏相當於n−1$n-1$ 個元素的錯排問題。

綜上，當n$n$ 排在第k$k$ 位時共有Dn−2+Dn−1$D_{n-2}+D_{n-1}$ 種錯排方法，又k$k$ 有從1$1$ 到n−1$n-1$ 共n−1$n-1$ 種取法，我們可以得到：Dn=(n−1)(Dn−1+Dn−2)$D_n=(n-1)(D_{n-1}+D_{n-2})$

更多其他方法另請參見一位不願具名的同學的題解。

代碼

#include <iostream>
#include <cstring>
#include <cstdio>
#include <iomanip>
#include <cstdlib>

#define inf 0x3f3f3f3f
#define PB push_back
#define MP make_pair
#define fi first
#define se second
#define lowbit(x) (x & (-x))
#define rep(i, a, b) for (int i = (a); i <= (b); i++)
#define per(i, a, b) for (int i = (a); i >= (b); i--)
#define pr(x) cout << #x << " = " << x << " ";
#define prl(x) cout << #x << " = " << x << endl;
#define ZERO(X) memset((X), 0, sizeof(X))
#define ALL(X) (X).begin(), (X).end()
#define SZ(x) (int)x.size()

using namespace std;

typedef pair<int, int> PI;
typedef pair<pair<int, int>, int> PII;
typedef pair<pair<pair<int, int>, int>, int> PIII;
using ull = unsigned long long;
using ll = long long;
using ld = long double;
#define quickio                  \
    ios::sync_with_stdio(false); \
    cin.tie(0);                  \
    cout.tie(0)
#define debug(...) fprintf(stderr, __VA_ARGS__), fflush(stderr)
/*      debug("Precalc: %.3f\n", (double)(clock()) / CLOCKS_PER_SEC);
clock_t z = clock();
        solve();
        //debug("Test: %.3f\n", (double)(clock() - z) / CLOCKS_PER_SEC);
*/
template <typename T = int>
inline T read()
{
    T val = 0, sign = 1;
    char ch;
    for (ch = getchar(); ch < '0' || ch > '9'; ch = getchar())
        if (ch == '-')
            sign = -1;
    for (; ch >= '0' && ch <= '9'; ch = getchar())
        val = val * 10 + ch - '0';
    return sign * val;
}
int main()
{
    ll arr[25];
    arr[1] = 0;
    arr[2] = 1;
    ll jc[25];
    jc[1] = 1;
    jc[2] = 2;
    for (int i = 3; i <= 20; ++i)
    {
        arr[i] = (i - 1) * (arr[i - 1] + arr[i - 2]);
        jc[i] = jc[i - 1] * i;
    }

    int T;
    cin >> T;
    while (T--)
    {
        int n;
        cin >> n;
        cout << fixed << setprecision(2)
             << double(arr[n]) * 100.0 / jc[n]  << "%" << endl;
    }
    return 0;
}

修改於20180524，進一步方便自己回顧對錯排問題的理解。