[Nowcoder 2018ACM多校第一场B] Symmetric Matrix

题目大意： 给出两个数n, m， 求构造这么一种矩阵的方案数模m： 对称矩阵； 对角元素为0； 矩阵元素由0，1，2构成； 各行的和都是2。$(n \leq 10^5, m \leq 10^9)$

题目思路： 可以做一个模型转换， 把这个矩阵看成是某个图的邻接矩阵。 那么实际上是求一个图， 没有自环， 各个点的度数都是2的方案数。 考虑dp， 用$f[n]$表示考虑前n个点时的方案数。 有两种转移， 一种是从n-1个点中拿出一个点和它连两条边， 剩下的部分方案数是$f[n-2]$， 即$(n-1)f[n-2]$； 一种是从n-1个点中拿出k个点和它一起构成一个环， 即$\sum_{k=1}^{n-3}\frac{1}{2}\frac{(n-1)!}{k!}f[k]$， 这里多除以个2是因为这个是一个项链排列， 要考虑对称性。 考虑第二种情况的式子也是可以递推的， 设其为g[n]， 比较g[n]与g[n-1]的区别容易看出只是乘一个系数再加一项， $g[n]=(n-1)g[n]+\frac{(n-1)(n-2)}{2}f[n-3]$。 就可以$O(n)$的做了。

Code:

#include <map>
#include <set>
#include <map>
#include <cmath>
#include <queue>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>

using namespace std;

#define ll long long

const int N = (int)1e5 + 10;

int n, m;
ll f[N], g[N];

int main(){

    while (scanf("%d%d", &n, &m) != EOF){
        f[0] = 1 % m; f[2] = 1 % m;

        for (int i = 3; i <= n; i ++){
            g[i] = ((i - 1) * g[i - 1] % m + 1ll * (i - 1) * (i - 2) / 2 % m * f[i - 3] % m) % m;
            f[i] = ((i - 1) * f[i - 2] % m + g[i]) % m;
        }

        printf("%lld\n", f[n]);
    }

    return 0;
}