【CCF-CSP-2013-12-4有趣的数】 DP

2013-12-4-有趣的数

题目链接
http://118.190.20.162/view.page?gpid=T2
题解
$本题考察点在于观察满足性质的序列种类，发现种类并不多$
$一共有13个合法状态，于是我们只要可以推理出13个合法状态之间的转移方程就可以。$
$首先我们给出13个合法状态：$
$1由1转移----1状态$
$2由2转移----2状态$
$3由3转移----3状态$
$02由24转移----4状态$
$03由35转移----5状态$
$12由126转移----6状态$
$13由137转移----7状态$
$23由28转移----8状态$
$012由49转移----9状态$
$013由510转移----10状态$
$023由4811转移----11状态$
$123由612转移----12状态$
$0123由91113转移----13状态$
$我们定义dp[i][j]为数字长度为i状态为j的数字种类数$
$很明显初始状态的定义为：$
$dp[1][1]=1,dp[1][2]=1,dp[1][3]=1;$
$我们很容易得到状态转移方程：$
$dp[i][1]=dp[i-1][1];$
$dp[i][2]=dp[i-1][2];$
$dp[i][3]=dp[i-1][3];$
$dp[i][4]=dp[i-1][2]+2\ast dp[i-1][4];$
$dp[i][5]=dp[i-1][3]+2\ast dp[i-1][5];$
$dp[i][6]=dp[i-1][1]+dp[i-1][2]+2\ast dp[i-1][6];$
$dp[i][7]=dp[i-1][1]+dp[i-1][3]+2\ast dp[i-1][7];$
$dp[i][8]=dp[i-1][2]+dp[i-1][8];$
$dp[i][9]=dp[i-1][4]+2\ast dp[i-1][9];$
$dp[i][10]=dp[i-1][5]+2\ast dp[i-1][10];$
$dp[i][11]=dp[i-1][4]+dp[i-1][8]+2\ast dp[i-1][11];$
$dp[i][12]=dp[i-1][6]+2\ast dp[i-1][12];$
$dp[i][13]=dp[i-1][9]+dp[i-1][11]+2\ast dp[i-1][13];$
$最终dp[n][13]即为最终答案。$
$本题坑点：$
$①:同状态之间的转移一定要注意不能使用本状态不存在的元素$
$②:不同状态之间的转移注意先后关系$
$③:注意不能存在前导0的问题$
代码

/*
***********************************************
Author        :biubiubiu
Created Time  :2018/10/17/03:46
File Name     :F:\CSP\2013-12\4.cpp
************************************************
*/
#include<stdio.h>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn = 1005;
const int Mod =1000000007;
ll dp[maxn][14];//dp[i][j]表示到达i长度满足j状态的种类数
int main()
{
    /*--------------------
    一个
    1      1       ----1状态
    2      2       ----2状态
    3      3       ----3状态
    两个
    02     2 4     ----4状态
    03     3 5     ----5状态
    12     1 2 6   ----6状态
    13     1 3 7   ----7状态
    23     2 8     ----8状态
    三个
    012    4 9      ----9状态
    013    5 10     ----10状态
    023    4 8 11   ----11状态
    123    6 12     ----12状态
    四个
    0123   9 11 13  ----13状态
    */
    dp[1][1]=1,dp[1][2]=1,dp[1][3]=1;
    int n;
    scanf("%d",&n);
    for(int i=2;i<=n;i++)
    {
        dp[i][1]=dp[i-1][1];
        dp[i][2]=dp[i-1][2];
        dp[i][3]=dp[i-1][3];
        dp[i][4]=dp[i-1][2]+2*dp[i-1][4];
        dp[i][5]=dp[i-1][3]+2*dp[i-1][5];
        dp[i][6]=dp[i-1][1]+dp[i-1][2]+2*dp[i-1][6];
        dp[i][7]=dp[i-1][1]+dp[i-1][3]+2*dp[i-1][7];
        dp[i][8]=dp[i-1][2]+dp[i-1][8];
        dp[i][9]=dp[i-1][4]+2*dp[i-1][9];
        dp[i][10]=dp[i-1][5]+2*dp[i-1][10];
        dp[i][11]=dp[i-1][4]+dp[i-1][8]+2*dp[i-1][11];
        dp[i][12]=dp[i-1][6]+2*dp[i-1][12];
        dp[i][13]=dp[i-1][9]+dp[i-1][11]+2*dp[i-1][13];
        for(int j=1;j<=13;j++) dp[i][j]%=Mod;
    }
    printf("%lld",dp[n][13]);
    return 0;
}

---