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

链接

http://115.28.138.223/view.page?gpid=T2

问题描述

　　我们把一个数称为有趣的，当且仅当：
　　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

分析

这题难点在于要想到使用动态规划，并找出合适的状态、转移方程。

状态如下：

状态[i][j]表示在前i位中，根据j的不同，分为如下几种情况
0－－用了2，剩0，1，3
1－－用了0，2，剩1，3
2－－用了2，3，剩0，1
3－－用了0，1，2，剩3
4－－用了0，2，3，剩1
5－－全部用了

为什么是这5种情况，也可自己分析出来:分析要写出长度为i位的满足要求的数时，
前i-1位可能包含哪些情况，在前i-1位可能的情况的最后添加1位即可写出i位的。
注意，前i-1位可以不合法(只出现4个数字中的一部分)，但要保证在前i-1位后面添加1位后合法！

状态转移方程可参考代码

已ac

//ccf 201312-4 有趣的数
//动态规划题
#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <algorithm>
#include <math.h>
#include <queue>
#include <stack>
#include <map>
#define INF 123123123
#define INS 0.000005
#define MAX_INDEX 1005
using namespace std;


long long result[MAX_INDEX][6];
/*
状态[i][j]表示在前i位中，是如下的情况
为什么是这5种情况，也可自己分析出来:分析要写出长度为i位的满足要求的数时，
前i-1位可能包含哪些情况，在前i-1位可能的情况的最后添加1位即可写出i位的。
注意，前i-1位可以不合法(只出现4个数字中的一部分)
0－－用了2，剩0，1，3
1－－用了0，2，剩1，3
2－－用了2，3，剩0，1
3－－用了0，1，2，剩3
4－－用了0，2，3，剩1
5－－全部用了*/
long long mods = 1000000007;
int main()
{
    int n;
    cin >> n;
    for(int j = 0; j<6; j++)
        result[0][j] = 0;
    /*result[4][0] = 1;//2222
    result[4][1] = 7;//2000,2002,2020,2200,2022,2202,2220
    result[4][2] = 3;//2333,2233,2223
    result[4][3] = 5;//2012,2021,2201,2001,2011
    result[4][4] = 9;//2003,2030,2300,2330,2303,2033,2023,2203,2230
    result[4][5] = 3;//2013,2301,2031*/
    for(int i=1; i<=n; i++)
    {
        result[i][0] = 1;
        result[i][1] = (result[i-1][0]+result[i-1][1]*2) % mods;
        result[i][2] = (result[i-1][2]+result[i-1][0]) % mods;
        //2只能在3前，前i-1位已有3，所以只能在最后添加3
        result[i][3] = (result[i-1][3]*2+result[i-1][1]) % mods;
        result[i][4] = (result[i-1][4]*2+result[i-1][1]+result[i-1][2]) % mods;
        result[i][5] = (result[i-1][3] + result[i-1][4] + result[i-1][5]*2) % mods;
    }
    cout << result[n][5] << endl;
    return 0;
}