ccf 有趣的数

问题描述

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

分析：

参考博客：http://blog.csdn.net/liuyongjie985/article/details/44516805

首先我们还是得从题目分析

0不能再首位是吧？那首位只能是1,2,3中的一个呗是吧？你放1看看行不，肯定不行啊，所有0得在1的前面，你放1了让0情何以堪？放2可以，放3的话你让2情何以堪？

所以首位只能放2对吧？

现在来说status这个二维数组，status[i][j]到第i个位置满足第j种状态的所有可能数，总共最多有6种状态：

之后要添加元素，是加在最后一个。

0 -- 0 1 (2) 3
1 -- (0) 1 (2) 3
2 -- 0 1 (2) (3)
3 -- (0) (1) (2) 3
4 -- (0) 1 (2) (3)
5 -- (0) (1) (2) (3)

比如说status[5][4] = 70的意思就是，到第5位（首位是第1位！）为止，保持第4种状态的数一共有70个（咳咳，当然数据是我瞎掰的）

在括号里的数说明前面的k位中，只出现了括号里的数，你可以发现6个状态中2都是被括号括起来的，没办法，刚刚说了2必须是首位啊

而且这是你在填第k+1位的时候能保持的仅有的6个状态，其他的状态都是无效的

比如0（1）（2) 3就是无效的，那是不可能的啊，假如说你现在填的是第k位吧，那么你填完第k位的数后的状态是只包括1和2，那你想想，怎么可能只出现1和2呢？那你让0情何以堪？0就没地方放了。

第一个for循环不讲

第二个for循环n次だろ？很明显这个就对应n位

然后看第一个式子

status[i][0] = 1;

到第i位为止，第i位填的数字使整个数字要满足状态0，那么这样的可能有几个？1个呗，你要能填出2以外的数来我就去（qia shi la）

status[i][1] = (status[i - 1][1] * 2 + status[i - 1][0]) % mod;

到第i位为止，填第i位的数，使整个数保持状态1。要满足状态1，也就是只出现0和2，那你想前i-1位应该是个什么状态，要么是只有2，要么是0,2对吧？就两种状态

然后如果前i-1位是状态0的话（只有2），那你这一位只能填0了嘛，所以status[i - 1][0]就是这样来的，如果前i-1位是状态1（有0,2），那么你现在有两种选择，填0或者是2两种情况，status[i - 1][1] * 2就是这样来的，然后把两种情况相加，就表示到第i位为止，状态为1的数共有status[i][1]这么多个

然后后面的都差不多啦，因为每一位可能会进入6种不同的状态，所以都要算出来

最后输出System.out.println(status[number][5])，因为题目要求0,1,2,3必须至少出现一次，所有我们要输出状态5的数的个数

每次都求余，这个很好理解吧？因为题目只要求求余数就可以了呀，中间结果不断取余不会有影响的

采用数位dp求解，到某一个阶段的各个状态递推。

#include <iostream>
#include <algorithm>
#include <string>
#include <list>
#include <vector> 
#include <iterator>
#include<string.h>
using namespace std;
#define mod 1000000007 
int main(void)  
{  
  int n;
  cin>>n;
  long long dp[n+1][6];
  memset(dp,0,sizeof(dp));
  for (int i=1;i<n+1;i++)
  {
  	dp[i][0]=1%mod;//2
  	dp[i][1]=(dp[i-1][0]+dp[i-1][1]*2)%mod;//0 2
  	dp[i][2]=(dp[i-1][0]+dp[i-1][2])%mod;//2 3
  	dp[i][3]=(dp[i-1][1]+dp[i-1][2]+dp[i-1][3]*2)%mod;//0 2 3
  	dp[i][4]=(dp[i-1][1]+dp[i-1][4]*2)%mod;//0 1 2
  	dp[i][5]=(dp[i-1][3]+dp[i-1][4]+dp[i-1][5]*2)%mod;//0 1 2 3
  }
  cout<<dp[n][5];
   return 0;  
}