不降数
题意:
定义
n
n
n 位不降数为从左到右的组成的序列是非递减的
给定
n
n
n,求长度为
n
n
n 的不降数有多少
思路:
对于
n
n
n 比较大,但是模数不大(一般不超过
1
0
6
10^6
106)的计数
大概率是要用到
L
u
c
a
s
Lucas
Lucas 定理的
可以发现不降数的每一位都是
[
1
,
9
]
[1,9]
[1,9] 之间的,看成序列,它至多有
8
8
8 次可以上升的机会
也就是序列的差分和的范围是
[
0
,
8
]
[0,8]
[0,8],枚举差分和
那么也就是
n
n
n 个位置可以选择上升,要选
i
i
i 个位置,每个位置可以重复选择
题目转化成了
n
n
n 个位置选
i
i
i 个元素,每个元素可以多次选择
答案即为
∑
i
=
0
i
=
8
C
(
n
−
1
+
i
,
i
)
\sum_{i=0}^{i=8} C(n-1+i,i)
∑i=0i=8C(n−1+i,i)
用
L
u
c
a
s
Lucas
Lucas 定理就可以快速得到答案了
这道题也可以
d
p
dp
dp + 优化取模卡过去
code:
#include<bits/stdc++.h>
#define endl '\n'
#define ll long long
#define ull unsigned long long
#define ld long double
#define all(x) x.begin(), x.end()
#define mem(x, d) memset(x, d, sizeof(x))
#define eps 1e-6
using namespace std;
const int maxn = 1e6 + 9;
const int mod = 100019;
const int inf = 0x3f3f3f3f;
const ll INF = 0x3f3f3f3f3f3f3f3f;
ll n, m;
ll fac[maxn], inv[maxn];
void init(){
fac[0] = fac[1] = 1;
inv[0] = inv[1] = 1;
for(int i = 2; i <= maxn - 9; ++i)
fac[i] = fac[i-1] * i % mod,
inv[i] = (mod - mod / i) * inv[mod % i] % mod;
for(int i = 2; i <= maxn - 9; ++i)
inv[i] = inv[i] * inv[i-1] % mod;
}
ll C(ll n, ll m){
if(n < m) return 0ll;
return fac[n] * inv[m] % mod * inv[n-m] % mod;
}
ll Lucas(ll n, ll m){
if(!m) return 1ll;
return Lucas(n / mod, m / mod) * C(n % mod, m % mod) % mod;
}
void work()
{
init();
cin >> n;
ll ans = 0;
for(int i = 0; i < 9; ++i){
(ans += Lucas(n - 1 + i, i)) %= mod;
}
cout << ans;
}
int main()
{
ios::sync_with_stdio(0);
// int TT;cin>>TT;while(TT--)
work();
return 0;
}