题目:
问题描述
我们把一个数称为有趣的,当且仅当:
1. 它的数字只包含0, 1, 2, 3,且这四个数字都出现过至少一次。
2. 所有的0都出现在所有的1之前,而所有的2都出现在所有的3之前。
3. 最高位数字不为0。
因此,符合我们定义的最小的有趣的数是2013。除此以外,4位的有趣的数还有两个:2031和2301。
请计算恰好有n位的有趣的数的个数。由于答案可能非常大,只需要输出答案除以1000000007的余数。
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
题解:对于给定的n位,我们将0,1,2,3四个数字填入n个空位中,我们将0,1组合为一组,2,3组合为一组,考虑(0,1)组,我们假设共有a个数字,由于0不能在第一位,所以总有n-1个空位可选,也就是C(n-1,a)中情况,剩余空位自然填(2,3)组,然后考虑每个组中的情况,由于每个数字至少出现一次,所以(0,1)组内排列共有(a-1)中情况(所有情况n+1中减去全0和全1的情况),(2,3)组组内排列共有(n-a-1)中。所以答案为
然后就是求解的问题,对于组合数,由于我们要连续求解,所以每次都可以利用上次的结果。比如,
所以只要将C(n-1,1)作为初始值,每次乘一个数,再除以一个数,即可得到下一个组合数,对于取模问题,我们可以用逆元来求解,这里有一个线性递推求逆元的方法:
//p表示取的模数,MAX表示要求的最大值
inv[1] = 1;
for (int i = 2; i<MAX; i++)
inv[i] = inv[P%i]*(P-P/i)%P;
这样问题就解决了,初始化逆元O(n),循环一次求解O(n),妥妥的。
代码:
#include <bits/stdc++.h>
using namespace std;
const int MOD=1000000007;
const int MAX=1000+10;
typedef long long LL;
LL inv[MAX];
void init(int P)
{
//p表示取的模数,MAX表示要求的最大值
inv[1] = 1;
for (int i = 2; i<MAX; i++)
inv[i] = inv[P%i]*(P-P/i)%P;
}
int main()
{
init(MOD);
int n;
while(scanf("%d",&n)!=-1)
{
LL ans=0;
LL col=n-1;
for(int i=2;i<=n-2;i++)
{
col=(col*(n-i)%MOD)*inv[i]%MOD;
ans=(ans+(col*(i-1)*(n-i-1)%MOD))%MOD;
}
cout<<ans<<endl;
}
return 0;
}
不是很懂用DP的大佬们。。。