CCF201312-4 有趣的数(动态规划)

试题编号：	201312-4
试题名称：	有趣的数
时间限制：	1.0s
内存限制：	256.0MB
问题描述：	问题描述 　　我们把一个数称为有趣的，当且仅当： 　　1. 它的数字只包含0, 1, 2, 3，且这四个数字都出现过至少一次。 　　2. 所有的0都出现在所有的1之前，而所有的2都出现在所有的3之前。 　　3. 最高位数字不为0。 　　因此，符合我们定义的最小的有趣的数是2013。除此以外，4位的有趣的数还有两个：2031和2301。 　　请计算恰好有n位的有趣的数的个数。由于答案可能非常大，只需要输出答案除以1000000007的余数。 输入格式 　　输入只有一行，包括恰好一个正整数n (4 ≤ n ≤ 1000)。 输出格式 　　输出只有一行，包括恰好n 位的整数中有趣的数的个数除以1000000007的余数。 样例输入 4 样例输出 3

看到“答案可能非常大”就知道暴力是不可能暴力的了，然而身为一枚蒟蒻，借鉴来了一种动规的做法（蒟蒻的事，能算偷么？），即第n位数的解来源于第n-1位的最优解。

对第n位而言，经过归纳，有六种状态：

状态	已用	剩余
0	2	0 1 3
1	2 0	1 3
2	2 3	0 1
3	2 0 1	3
4	2 0 3	1
5	2 0 1 3

接着，我们需要分析，每个状态可以通过哪些状态转移得到：

状态0：只有一种可能，即状态0本身；

状态1：①第n-1位是状态0的情况下，第n位填0（以下略写）；②状态1填2或0（2与0之间并无冲突，故在第n-1位已达到状态1的情况下，在第n位填2或0都是合法的，由于有两种可能，所以需要将在状态1情况下的第n-1位的可能数乘以2）；

状态2：①状态0填3；②状态2填3（由于2与3有冲突，且3必须出现在2的后面，所以在2、3已经出现的情况下，只能在第n位填3，否则违反了规则）；

状态3：①状态1填1；②状态3填2或1，故乘以2；

状态4：①状态1填3；②状态2填0；③状态4填0或3，故乘以2；

状态5：①状态3填3；②状态4填1；③状态5填1或3，故乘以2；

 

由此，所有状态转移的公式就总结完毕了，下面是代码：

#include <bits/stdc++.h>
#define ll long long
#define INF 0x3f3f3f3f
using namespace std;
const int maxn = 1e4+10;
const ll mod = 1e9+7;
int n;
ll dp[maxn][6];

/*
state used   rest
0     2      013
1     20     13
2     23     01
3     201    3
4     203    1
5     0123
*/

void solve()
{
    cin >> n;
    for(int i = 1; i <= n; ++i)
    {
        int j = i-1;
        dp[i][0] = 1;
        dp[i][1] = (dp[j][0] + dp[j][1]*2) % mod; //由于时时有可能会溢出
        dp[i][2] = (dp[j][0] + dp[j][2]) % mod; //故在递推的同时就进行取模
        dp[i][3] = (dp[j][1] + dp[j][3]*2) % mod;
        dp[i][4] = (dp[j][1] + dp[j][2] + dp[j][4]*2) % mod;
        dp[i][5] = (dp[j][3] + dp[j][4] + dp[j][5]*2) % mod;
    }
    cout << dp[n][5]; //解为第n位的状态5的可能数
}

int main()
{
	solve();
	return 0;
}