P4345 [SHOI2015]超能粒子炮·改（卢卡斯定理推式子）

题目链接：点击这里

题目大意：
$t$ 组样例，每组样例给定 $n,k$ 求：
$$\sum _{i=0}^k{C}_n^i\mathrm{m}\mathrm{o}\mathrm{d}2333$$

题目分析：
设 $F(n,k)=\sum _{i=0}^k{C}_n^i,p=2333$ 有：
$$F(n,k)=\sum _{i=0}^k{C}_n^i$$
$$=\sum _{i=0}^k{C}_{n/p}^{i/p}{C}_{n\%p}^{i\%p}$$
对于这个式子我们观察 $C_{n/p}^{i/p}$ 发现这东西跟整除分块很像，我们可以对其进行分块对 $C_{n/p}^{i/p}$ 的每一个 $i/p$ 进行合并：
$$=C_{n/p}^0\sum _{i=0}^{p-1}{C}_{n\%p}^i+C_{n/p}^1\sum _{i=0}^{p-1}{C}_{n\%p}^i+...+C_{n/p}^{k-1}\sum _{i=0}^{p-1}{C}_{n\%p}^i+C_{n/p}^{k/p}\sum _{i=0}^{k\%p}{C}_{n\%p}^i$$
前面 $k-1$ 个是整块，最后一个不一定是完整的将其称之为余项，我们先考虑整块的处理，先将 $\sum _{i=0}^k{C}_{n\%p}^i$ 提出了可得：
$$=\sum _{i=0}^{p-1}{C}_{n\%p}^i\sum _{j=0}^{k/p-1}{C}_{n/p}^j$$
$$=F(n\%p,p-1)F(n/p,k/p-1)$$
接下来我们考虑处理余项：
$$C_{n/p}^{k/p}\sum _{i=0}^{k\%p}{C}_{n\%p}^i$$
$$=C_{n/p}^{k/p}F(n\%p,k\%p)$$
综上有：
$$F(n,k)=F(n\%p,p-1)F(n/p,k/p-1)+C_{n/p}^{k/p}F(n\%p,k\%p)$$
对于 $n<p$ 且 $k<p$ 的情况由 $F(n,k)$ 的定义有：
$$F(n,k)=F(n,k-1)+C_n^k$$
边界为 $F(n,0)=C_n^0=1$
因此对于 $F(n\%p,p-1),F(n\%p,k\%p)$ 我们可以打表求出，对于 $C_{n/p}^{k/p}$ 我们可以直接通过卢卡斯定理求出

时间复杂度为 $O(tlo{g}^{2}n)$

具体细节见代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<set>
#include<map>
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
ll read()
{
	ll res = 0,flag = 1;
	char ch = getchar();
	while(ch<'0' || ch>'9')
	{
		if(ch == '-') flag = -1;
		ch = getchar();
	}
	while(ch>='0' && ch<='9')
	{
		res = (res<<3)+(res<<1)+(ch^48);//res*10+ch-'0';
		ch = getchar();
	}
	return res*flag;
}
const int maxn = 2333+5;
const int mod = 2333;
const double pi = acos(-1);
const double eps = 1e-8;
ll c[maxn][maxn],f[maxn][maxn];

ll lucas(ll n,ll m)
{
	if(!m) return 1;
	return lucas(n/mod,m/mod)*c[n%mod][m%mod]%mod;
}
ll F(ll n,ll k)
{
	if(k < 0) return 0;
	if(!n || !k) return 1;
	if(n < mod && k < mod) return f[n][k];
	return (F(n/mod,k/mod-1)*f[n%mod][mod-1]%mod+lucas(n/mod,k/mod)*f[n%mod][k%mod]%mod)%mod;
}
int main()
{
	ll t = read();
	c[0][0] = f[0][0] = 1;
	for(int i = 1;i < maxn;i++)
		c[i][i] = c[i][0] = f[i][0] = 1;
	for(int i = 1;i < maxn;i++)
		for(int j = 1;j < maxn;j++)
			c[i][j] = (c[i-1][j]+c[i-1][j-1])%mod;
	for(int i = 0;i < maxn;i++)
		for(int j = 1;j < maxn;j++)
			f[i][j] = (c[i][j]+f[i][j-1])%mod; 
	while(t--)
	{
		ll n = read(),k = read();
		printf("%lld\n",F(n,k));
	}
	return 0;
}