51nod 1556 计算（默慈金数）

默慈金数

部分摘自百度文库
默慈金数是在数学中，一个给定的数n的默慈金数是“在一个圆上的n个点间，画出彼此不相交的弦的全部方法的总数”。
也可以想象成在一个“网格”上，若限定每步只能向右移动一格(可以向右上、向右下、向右)，并禁止移动到 y = 0 以下的地方”，则以这种走法用 n 步从 (0,0) 移动至 (n,0) 的可能形成的路径的总数为 n 的默慈金数。

前几个默慈金数

1, 2, 4, 9, 21, 51, 127, 323, 835, 2188, 5798, 15511, 41835, 113634, 310572, 853467, 2356779, 6536382, 18199284, 50852019, 142547559, 400763223, 1129760415, 3192727797, 9043402501, 25669818476, 73007772802, 208023278209, 593742784829…
第四个：9 种

递推公式

代码实现

	m[1] = 1;
    m[2] = 2;
    for (int i = 3; i < n; i++) {
        m[i] = (1l * (2 * i + 1) * m[i - 1] % mod + 1l * 3 * (i - 1) * m[i - 2] % mod) 
        		% mod * QuickPow(i + 2, mod - 2) % mod;
    }

1556 计算

有一个 1*n 的矩阵 固定第一个数为 1 其他填正整数 且相邻数的差不能超过 1 求方案数 %1e9+7 的结果

输入

一个数 n 表示 1 * n 的矩阵（n<=10^6）

输出

一个数 表示方案数 %1e9+7 的结果

输入样例

3

输出样例

5

代码

如果在x轴上就只有两种方案，其余有3种方案，所以每种方案 * 3 后减去，第一种方案即对应的默慈金数。

#include <bits/stdc++.h>

#define ll long long
using namespace std;
const int mod = 1e9 + 7;
const int MAXN = 1e6 + 10;
int n;
// f 答案，m 默慈金数
ll f[MAXN], m[MAXN];

ll QuickPow(ll x, ll N) {
    ll res = x % mod;
    ll ans = 1;
    while (N) {
        if (N & 1) {
            ans = ans * res % mod;
        }
        res = res * res % mod;
        N = N >> 1;
    }
    return ans;
}

int main() {
    cin >> n;
    m[1] = 1;
    m[2] = 2;
    f[1] = 2;
    for (int i = 3; i < n; i++) {
        m[i] = (1l * (2 * i + 1) * m[i - 1] % mod
                + 1l * 3 * (i - 1) * m[i - 2] % mod) % mod * QuickPow(i + 2, mod - 2) % mod;
    }
    for (int i = 2; i < n; i++) {
        f[i] = (1l * f[i - 1] * 3 - m[i - 1] + mod) % mod;
    }
    printf("%lld\n", f[n - 1]);
    return 0;
}
/*
n = 12
1        1        2        4        9       21       51      127      323      835     2188     5798
0        1        2        5       12       30       76      196      512     1353     3610     9713
0        0        1        3        9       25       69      189      518     1422     3915    10813
0        0        0        1        4       14       44      133      392     1140     3288     9438
0        0        0        0        1        5       20       70      230      726     2235     6765
0        0        0        0        0        1        6       27      104      369     1242     4037
0        0        0        0        0        0        1        7       35      147      560     2002
0        0        0        0        0        0        0        1        8       44      200      814
0        0        0        0        0        0        0        0        1        9       54      264
0        0        0        0        0        0        0        0        0        1       10       65
0        0        0        0        0        0        0        0        0        0        1       11
0        0        0        0        0        0        0        0        0        0        0        1
*/