ccf试题有趣的数，动态规划

关键：将二维数组看作二维函数，由前面的递推到后面
思路：之前尝试用深度搜索，一个个位置试探，可是本人只能求出一个最小的有趣数，而无法统计数量，后看到一些博客，大多乱七八糟，叙述不明，有一篇倒是解释的非常清楚，博主连接：https://blog.csdn.net/qq_37969433/article/details/86551488

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

贴代码：

#include<iostream>
using namespace std;
//动态分配，前面一步影响后面的分配。

/*
	将二维数组看作由两个因变量决定的函数
	任何一个n位有趣数，都可以看做是n-1位数，再在其后面追加0-3其中一个数而组成的

	要满足第一位为2,2在3前，0在1前
	数字只会出现六种情况，
	0.只有数字2，缺少0、1、3
	1.只有数字2、0，缺少1、3
	2.只有2、3，缺少0、1
	3.只有2、1、0，缺少3
	4.只有2、3、0，缺少1
	5.所有数字都有，已为有趣数,(末尾可以加1或3)

	最后需要得到的是状态5的n位数的个数，由已知逻辑关系可以发现一种递推关系
	记所求值为f（n,5），
	f(n,5)=f(n-1,5)*2+f(n-1,4)+f（n-1，3）  等于n-1位的状态5末尾加1或者3 加上 n-1位的状态4末尾加1 状态3加3
	f(n,4)=f(n-1,4)*2 +f(n-1,2)+f(n-1,1) 等于n-1位状态4的末尾加0或3 加上 状态2加0 或者状态1加3
	f(n,3)=f(n-1,1)+f(n-1,3)*2
	f(n,2)=f(n-1,2)+f(n-1,0)
	f(n,1)=f(n-1,1)*2+f(n-1,0)
	f(n,0)=f(n-1,0)
 */

int main()
{
	int n;
	long long mod = 1000000007;
	while (cin >> n)//输入位数
	{//动态创建二维数组,用二维数组来代替两个变量：位数和状态
		long long** f = new long long* [n + 1];//多申请一行防止溢出
		for (int i = 0; i < n + 1; i++)
		{
			f[i] = new long long[6];
			if (i == 0)
				for (int j = 0; j < 6; j++)
				{
					f[i][j] = 0;//0位数的任何状态个数都为0
				}
		}
		/*
		0.只有数字2，缺少0、1、3
		1.只有数字2、0，缺少1、3
		2.只有2、3，缺少0、1
		3.只有2、1、0，缺少3
		4.只有2、3、0，缺少1
		5.所有数字都有，已为有趣数,(末尾可以加1或3)*/
		for (int i = 1; i <= n; i++)
		{//从一位数开始,递推到n位数
		//不取余的话最后的数可能会超过界限，导致错误，由同余定理只要最后再取一次模所得还是相同的

			f[i][0] = 1;//i位数满足状态为0的只有一个数
			f[i][1] = (f[i - 1][0] + f[i - 1][1] * 2) % mod;
			f[i][2] = (f[i - 1][0] + f[i - 1][2]) % mod;
			f[i][3] = (f[i - 1][1] + f[i - 1][3] * 2) % mod;
			f[i][4] = (f[i - 1][1] + f[i - 1][2] + f[i - 1][4] * 2) % mod;
			f[i][5] = (f[i - 1][3] + f[i - 1][4] + f[i - 1][5] * 2) % mod;
			

		/*	
			f[i][0] = 1;//i位数满足状态为0的只有一个数
			f[i][1] = (f[i - 1][0] + f[i - 1][1] * 2) ;
			f[i][2] = (f[i - 1][0] + f[i - 1][2]);
			f[i][3] = (f[i - 1][1] + f[i - 1][3] * 2);
			f[i][4] = (f[i - 1][1] + f[i - 1][2] + f[i - 1][4] * 2);
			f[i][5] = (f[i - 1][3] + f[i - 1][4] + f[i - 1][5] * 2);
		*/}
		cout <<f[n][5]% mod<< endl;//
	}
	return 0;
}