有趣的数

• CCF
• 组合数

---

题干描述

201312-4

我们把一个数称为有趣的，当且仅当：

1.它的数字只包含 0,1,2,3，且这四个数字都出现过至少一次。

2.所有的 0都出现在所有的 1之前，而所有的2都出现在所有的 3 之前。

3.最高位数字不为0。

因此，符合我们定义的最小的有趣的数是 2013。

除此以外 ，4位的有趣的数还有两个：2031 和 2301。

请计算恰好有n位的有趣的数的个数。

由于答案可能非常大，只需要输出答案除以10^9+7 的余数。

输入格式

输入只有一行，包括恰好一个正整数n。

输出格式

输出只有一行，包括恰好 n 位的整数中有趣的数的个数除以 10^9+7 的余数。

数据范围

4≤n≤1000

输入样例：

4

输出样例：

3

思路分析

根据题干我们需要满足三个要求：

①0，1，2，3至少出现一次。 ②0在1之前，2在3之前。 ③0不为首位。

据此我们可以把问题分为若干类。例如：01出现2次，23出现n-2次；01出现3次，23出现n-3次……01出现n-2次，23出现2次。

01 （次数）	2	3	4	…………	n-2
23 (次数）	n-2	n-3	n-4	……………	2

也即是我们算出每一个组合的方案数相加即为总的方案数。

①对于k位01而言，假设前t位放0的话，后k-t位放1。相应的对于每一个t值，都分别对应一种组合。由于0至少出现一次所以对应又k-1种放法。

②通里得到23的不同的序列选法为n-k-1

综合上述，最终，我们可以得到方案数之和公式。

注：第二种方法

下面开始探索这样的状态有哪些，以及状态之间的相关性。

按照题意，

可以出现的数字类型有0、1、2、3且要求的是n位数字的个数，那么对于一个n位数字，

可以将其状态定义为：n位中有特定类型数字，且满足0在1前、2在3前的个数。

具体来说：f(n, 0123)表示前n位中同时出现0、1、2、3且满足0在1前、2在3前的数字个数；

f(n, 012)表示前n位中同时出现0、1、2且满足0在1前、2在3前的数字个数（当然这里没有3，所有2在3前可以视为妥妥能达到）。

定义了上述状态后，其实可以比较顺畅地想到状态机dp了。在进行状态转移之前，我们可以剔除掉一些无效状态，即肯定不会出现的状态：f(n, 0) 当n=1时，该状态可以出现，但n>1时，数字不能全是0且n位数大于1，由于本题中满足要求的数字至少4位，所以该状态可以剔除；f(n, 1) f(n, 3)由于1前面必须出现0，3前面必须出现2，所以不可能单独有1和3出现，该状态可剔除；……按照上述类似方法，最终可能出现的状态只能有f(n, 2), f(n, 02), f(n, 01), f(n, 23), f(n, 012), f(n, 023), f(n, 0123)，共6种

注：指向状态0123圈圈的箭头表示，0123状态只能由本身及012和023状态转移过来`

#include <cstring>
#include <algorithm>
#include <iostream>
using namespace std;

typedef long long LL;

const int N = 1010, mod = 1e9 + 7;
LL a2, a02, a23, a012, a023, a0123;
LL t2, t02, t23, t012, t023, t0123;//默认设为0
int n;

int main() {
    cin >> n;
    a2 = 1;
    for (int i = 2; i <= n; i++) {
        // 备份
        t2 = a2, t02 = a02, t23 = a23, t012 = a012, t023 = a023, t0123 = a0123;
        a2 = t2 % mod;
        // 系数02表示在末尾可以加上0或2两种，下同
        a02 = (t2 + 2 * t02) % mod;
        a23 = (t2 + t23) % mod;
        a012 = (t02 + 2 * t012) % mod;
        a023 = (t02 + t23 + 2 * t023) % mod
        a0123 = (t012 + t023 + 2 * t0123) % mod;
    }

    cout << a0123 << endl;

    return 0;
}