牛客练习赛80 C（Lucas定理

不降数
题意：
定义 $n$ 位不降数为从左到右的组成的序列是非递减的
给定 $n$，求长度为 $n$ 的不降数有多少
思路：
对于 $n$ 比较大，但是模数不大（一般不超过 $10^6$）的计数
大概率是要用到 $Lucas$ 定理的

可以发现不降数的每一位都是 $[1,9]$ 之间的，看成序列，它至多有 $8$ 次可以上升的机会
也就是序列的差分和的范围是 $[0,8]$，枚举差分和
那么也就是 $n$ 个位置可以选择上升，要选 $i$ 个位置，每个位置可以重复选择
题目转化成了 $n$ 个位置选 $i$ 个元素，每个元素可以多次选择
答案即为 ${\sum }_{i=0}^{i=8}C(n-1+i,i)$
用 $Lucas$ 定理就可以快速得到答案了

这道题也可以 $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;
}