【BZOJ4591】[SHOI2015]超能粒子炮·改（卢卡斯定理）

2024-04-12 18:03:14

【BZOJ4591】[SHOI2015]超能粒子炮·改（卢卡斯定理）

题面

BZOJ
洛谷

题解

感天动地！终于不是拓展卢卡斯了！我看到了一个模数，它是质数！！！
看着这个东西就感觉可以递归处理。
令\(f(n,k)\)表示答案。
\[\begin{aligned} f(n,k)&=\sum_{i=0}^k {n\choose i}\\ &=\sum_{i=0}^k {n/p\choose i/p}*{n\%p\choose i\%p}\\ &=\sum_{x=0}^{p-1}{n\%p\choose x}*\sum_{i=0}^k[i\%p=x]{n/p\choose i/p}\\ &=\sum_{x=0}^{p-1}{n\%p\choose x}*\sum_{i=0}^{(k-x)/p}{n/p\choose i}\\ &=\sum_{x=0}^{p-1}{n\%p\choose x}*f(n/p,(k-x)/p) \end{aligned}\]
前面那个东西可以提前预处理好前缀和，而后面那个东西最多递归两次。而递归层数也就最多\(6\)层。所以单次复杂度\(O(2^6)\)。卡卡常就洛谷rk1了。。。

#include<iostream>
#include<cstdio>
using namespace std;
#define ll long long
#define MOD 2333
#define MAX 2350
inline ll read()
{ll x=0;bool t=false;char ch=getchar();while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();if(ch=='-')t=true,ch=getchar();while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();return t?-x:x;
}
int C[MAX][MAX];
int f(ll n,ll k)
{if(n<MOD)return C[n][min(n,k)];if(!k)return 1;int ret=0,x=k%MOD,y=n%MOD;ret=C[y][min(y,x)]*f(n/MOD,k/MOD);if(k-x)ret=(ret+(C[y][y]-C[y][min(y,x)]+MOD)*f(n/MOD,(k-x-1)/MOD))%MOD;return ret;
}
ll n,k;
int main()
{for(int i=0;i<MAX;++i)C[i][0]=1;for(int i=1;i<MAX;++i)for(int j=1;j<=i;++j)C[i][j]=(C[i-1][j]+C[i-1][j-1])%MOD;for(int i=0;i<MAX;++i)for(int j=1;j<=i;++j)C[i][j]=(C[i][j]+C[i][j-1])%MOD;int T=read();while(T--){n=read();k=read();printf("%d\n",f(n,k));}return 0;
}

转载于:https://www.cnblogs.com/cjyyb/p/10180457.html

【BZOJ4591】[SHOI2015]超能粒子炮·改（卢卡斯定理）相关推荐

[BZOJ4591][SHOI2015]超能粒子炮·改(Lucas定理+数位DP)
大组合数取模可以想到Lucas,考虑Lucas的意义,实际上是把数看成P进制计算. 于是问题变成求1~k的所有2333进制数上每一位数的组合数之积. 数位DP,f[i][0/1]表示从高到低第i位,这 ...
BZOJ4591 SHOI2015超能粒子炮·改（卢卡斯定理+数位dp）
注意到模数很小,容易想到使用卢卡斯定理,即变成一个2333进制数各位组合数的乘积.对于k的限制容易想到数位dp.可以预处理一发2333以内的组合数及组合数前缀和,然后设f[i][0/1]为前i位是否卡 ...
[2019.2.24]BZOJ4591 [Shoi2015]超能粒子炮·改
以下除法一律为整除. 求\(\sum_{i=0}^kC_n^i\ mod\ p,p=2333\) 设\(f(i,j)=\sum_{k=0}^jC_i^k\) \(f(n,k)=\sum_{i=0}^k ...
P4345 [SHOI2015]超能粒子炮·改
P4345 [SHOI2015]超能粒子炮·改题意: 求解式子∑i=0kCni%p\sum_{i=0}^{k}C_{n}^{i} \% p∑i=0kCni%p n,k<=1e18 题解: ...
[SHOI2015] 超能粒子炮·改
题目链接:戳我设\(f[n][k]=\sum_{i=0}^kC_n^i\pmod {2333}\) 那么根据卢卡斯定理我们知道\(f[n][k]=\sum_{i=0}^kC_{n/p}^{i/p}\ ...
【[SHOI2015]超能粒子炮·改】
就是运用\(Lucas\)推一个柿子首先是前置芝士\(Lucas\)定理 \[C_{n}^{m}\%p=C_{n/p}^{m/p}*C_{n\%p}^{m\%p}\%p\] 至于证明我建议去问一下 ...
loj2038「SHOI2015」超能粒子炮・改（卢卡斯的理解和延伸+数位dp）
我们发现暴力+裸的卢卡斯可以骗到50分所以就应该朝卢卡斯的方向去想. 卢卡斯是把n,k拆成两个p进制数进行统计的. 注意到如果在某一位ai>nia_i>n_iai>ni,a表示 ...
P4345-[SHOI2015]超能粒子炮·改【Lucas定理,类欧】
正题题目链接:https://www.luogu.com.cn/problem/P4345 题目大意 TTT组询问,给出n,kn,kn,k求 ∑i=0k(ni)\sum_{i=0}^{k}\bino ...
超能粒子炮·改 HYSBZ - 4591
题目链接: http://www.lydsy.com/JudgeOnline/problem.php?id=4591 就是求: C(n,0)+C(n,1)+... ...

最新文章

热门文章