试题编号: | 201312-4 |
试题名称: | 有趣的数 |
时间限制: | 1.0s |
内存限制: | 256.0MB |
问题描述: | 问题描述 我们把一个数称为有趣的,当且仅当: 输入格式 输入只有一行,包括恰好一个正整数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;
}