loj 2038 / 洛谷 P4345 [SHOI2015] 超能粒子炮・改 题解
好玩的推式子
题目描述
曾经发明了脑洞治疗仪与超能粒子炮的发明家 SHTSC 又公开了他的新发明:超能粒子炮・改——一种可以发射威力更加强大的粒子流的神秘装置。
超能粒子炮・改相比超能粒子炮,在威力上有了本质的提升。它有两个参数 \(n\)、\(k\),它会向每个编号为 \(0\) 到 \(k\)(包含两端)的位置 \(i\) 发射威力为 \(\mathrm{C}_n^i \bmod 2333\) 的粒子流。
现在 SHTSC 给出了他的超能粒子炮・改的参数,让你求出其发射的粒子流的威力之和除以 \(2333\) 所得的余数。
输入格式:
第一行一个整数 \(t\) 表示数据组数。
之后 \(t\) 行,每行两个整数 \(n\)、\(k\),含义如题面描述。输出格式:
\(t\) 行,每行一个整数,表示其粒子流的威力之和模 \(2333\) 的值。
输入输出样例
输入样例:
3
5 5
10 7
1145 14
输出样例:
32
968
763
数据范围与约定
对于 \(10\%\) 的数据,\(t=1\),\(n,k\le 1000\);
对于 \(30\%\) 的数据,\(t=1\),\(n,k\leq 1000000\);
对于 \(50\%\) 的数据,\(t=1\),\(n\le 10^{18},k\le 1000\);
对于 \(70\%\) 的数据,\(t\le 100\),\(n,k\le 10^{18}\);
对于 \(100\%\) 的数据,\(t\le 100000\),\(n,k\le 10^{18}\)。
题解:
注:本文中所有的除法 \(/\) 都向下取整,所有的百分号 \(\%\) 都表示取模。
这个题求的是 \(\sum_{i=0}^k\mathrm{C}_n^i\bmod 2333\)。但是模数是 \(2333\) 因此可以考虑 Lucas 定理,即 \(\mathrm{C}_n^m\% p=\mathrm{C}_{n\% p}^{m\% p}\mathrm{C}_{n/p}^{m/p}\)。
我们把上面的和式推导一下,则为
\]
然后我们发现,整个过程中 \(n/2333\) 和 \(n\% 2333\) 是不变的。只需要关注 \(i/2333\) 和 \(i\%2333\) 的变化规律。
而对于连续的 \(i\in[2333k,2333k+2333)\) 它们的 \(i/2333\) 是相同的,\(i\%2333\in[0,2333)\),所以我们把需要求和的 \(k\) 个数分成 \(\left\lceil\frac{k}{2333}\right\rceil\) 段。其中前 \(\left\lfloor\frac{k}{2333}\right\rfloor\) 段一定是完整的,因此可以表示为
\]
对于加号后面的式子,\(i/2333=0\),所以是对后面一个式子求和,因此可以用杨辉三角预处理,并求出前缀和,\(O(1)\) 解决。
对于中间一个式子 \(\sum_{i=0}^{2332}\mathrm{C}_{n\%2333}^i\) ,因为 \(n\%2333<2333\) ,同理用杨辉三角。
令
\]
则原式为
\]
对于剩下的一个式子,又转化为了一个求和的子问题,因此我们递归解决。递归的边界是 \(k<2333\)。
因此可以在 \(2333^2+O(t\log k)\) 的时间复杂度内解决这个问题。
Code:
#include<cstdio>
#define ll long long
ll C[2333][2333];
ll calc(ll n,ll t)//0~t的C_n^t
{
ll ans=0,p=n%2333;
if(t/2333)
ans=C[p][p]*calc(n/2333,t/2333-1)%2333;
else
return C[p][t%2333];
ll a=n/2333,b=t/2333,tmp=1;
while(a>=2333||b>=2333)
{
if(b%2333)
tmp=tmp*(C[a%2333][b%2333]-C[a%2333][b%2333-1]+2333)%2333;
a/=2333,b/=2333;
}
if(b)
tmp=tmp*(C[a][b]-C[a][b-1]+2333)%2333;
ans=(ans+C[p][t%2333]*tmp)%2333;
return ans;
}
int main()
{
C[0][0]=1;
for(int i=1;i<=2332;++i)
{
C[i][0]=1;
for(int j=1;j<=i;++j)
C[i][j]=(C[i-1][j-1]+C[i-1][j])%2333;
}
for(int i=0;i<=2332;++i)
for(int j=1;j<=2332;++j)
C[i][j]=(C[i][j-1]+C[i][j])%2333;
int T;
ll n,k;
scanf("%d",&T);
while(T--)
{
scanf("%lld%lld",&n,&k);
printf("%lld\n",calc(n,k));
}
return 0;
}
loj 2038 / 洛谷 P4345 [SHOI2015] 超能粒子炮・改 题解的更多相关文章
- 洛谷 P4345 [SHOI2015]超能粒子炮·改 解题报告
P4345 [SHOI2015]超能粒子炮·改 题意 求\(\sum_{i=0}^k\binom{n}{i}\),\(T\)组数据 范围 \(T\le 10^5,n,j\le 10^{18}\) 设\ ...
- [洛谷P4345][SHOI2015]超能粒子炮·改
题目大意:给你$n,k$,求:$$\sum\limits_{i=0}^k\binom n i\pmod{2333}$$题解:令$p=2333,f(n,k)\equiv\sum\limits_{i=0} ...
- bzoj4591 / P4345 [SHOI2015]超能粒子炮·改
P4345 [SHOI2015]超能粒子炮·改 题意:求$\sum_{i=1}^{k}C(n,i)\%(P=2333)$ 肯定要先拆开,不然怎么做呢(大雾) 把$C(n,i)$用$lucas$分解一下 ...
- P4345 [SHOI2015]超能粒子炮·改 Lucas
\(\color{#0066ff}{ 题目描述 }\) 曾经发明了脑洞治疗仪与超能粒子炮的发明家 SHTSC 又公开了他的新发明:超能粒子炮・改--一种可以发射威力更加强大的粒子流的神秘装置. 超能粒 ...
- P4345 [SHOI2015]超能粒子炮·改
传送门 看到数据和模数大小就知道要上 lucas 了 然后开始愉快地推公式: 答案为 $\sum _{i=0}^kC_{n}^{i}\ (mod\ 2333)$ 设 $f [ i ] [ j ] = ...
- 【BZOJ4591】[SHOI2015]超能粒子炮·改 (卢卡斯定理)
[BZOJ4591][SHOI2015]超能粒子炮·改 (卢卡斯定理) 题面 BZOJ 洛谷 题解 感天动地!终于不是拓展卢卡斯了!我看到了一个模数,它是质数!!! 看着这个东西就感觉可以递归处理. ...
- Bzoj 4591: [Shoi2015]超能粒子炮·改 数论,Lucas定理,排列组合
4591: [Shoi2015]超能粒子炮·改 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 178 Solved: 70[Submit][Stat ...
- bzoj 4591: [Shoi2015]超能粒子炮·改 [lucas定理]
4591: [Shoi2015]超能粒子炮·改 题意:多组询问,求 \[ S(n, k) = \sum_{i=0}^n \binom{n}{i} \mod 2333,\ k \le n \le 10^ ...
- BZOJ_4591_[Shoi2015]超能粒子炮·改_Lucas定理
BZOJ_4591_[Shoi2015]超能粒子炮·改_Lucas定理 Description 曾经发明了脑洞治疗仪&超能粒子炮的发明家SHTSC又公开了他的新发明:超能粒子炮·改--一种可以 ...
随机推荐
- ubuntu16.04 安装openpose
安装 Anaconda3 Tensorflow-cpu python3tensorflow 1.4.1+opencv3, protobuf, python3-tk ================== ...
- 关于sleep的理解
unix是按时间片轮转调度, windows是抢占式调度 以吃蛋糕为例子,10个人吃蛋糕,如果是unix下, 假设开始时,每个人都处于就绪状态,那么操作系统调度大家排好队,按顺序吃,每个人吃1分钟, ...
- SqlServer------范式小结
说明:大多数初学者对于关系数据库中的范式很是头疼,我本人也是,所以今天又看了视频,总结了一下内容,尽量语言通俗易懂,少用专业术语以及概念. 首先要理解几个键值. 超键:在关系模式中,能唯一标识元组的属 ...
- MySql 获取服务提供的sakila数据库(Example Databases)
关于这个数据库也就是样例数据库,数据库,数据库,最可怕的就是没有数据了,对吧?没有数据你学个什么呀. 可是,没有数据,咱会自己insert,那只能适用于初学者.对于数据库的优化方面的学习,还是有大数据 ...
- AIO和NIO的理解
AIO: AIO 背后的基本思想是允许进程发起很多 I/O 操作,而不用阻塞或等待任何操作完成,可以继续做 另外的事情,等I/O操作完成,内核会通过函数回调或者信号机制通知用户进程.这样很大程度提高了 ...
- 不用EL表达式---实现product页面显示
产品页面显示 静态页面如下: <%@ page language="java" contentType="text/html; charset=UTF-8" ...
- POJ3026 Borg Maze(bfs求边+最小生成树)
Description The Borg is an immensely powerful race of enhanced humanoids from the delta quadrant of ...
- Bitmap类、BitmapFactory及BitmapFactory类中的常用方法
1.Bitmap 1.1非静态方法 public void recycle()——回收位图占用的内存空间,把位图标记为Dead public final boolean isRecycled() —— ...
- Extjs Hello extjs
<html > <head runat="server"> <title></title> <link rel="s ...
- C#NPOI.RabbitMQ.EF.Attribute.HttpRuntime.Cache.AD域.List<T>根据指定字段去重.前端JQuery.Cache.I18N(多语言).data-xx(自定义属性)
使用NPOI 操作Excel 个人使用的电脑基本默认安装Excel 操作起来 调用Excel的组件便可.如果是一台服务器.没有安装Excel,也就无法调用Excel组件. 在此推荐第三方插件.NPOI ...