好玩的推式子

题目描述

曾经发明了脑洞治疗仪与超能粒子炮的发明家 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}\)。

我们把上面的和式推导一下,则为

\[\sum_{i=0}^k\mathrm{C}_{n/2333}^{i/2333}\mathrm{C}_{n\%2333}^{i\%2333}
\]

然后我们发现,整个过程中 \(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\) 段一定是完整的,因此可以表示为

\[\sum_{t=0}^{\left\lfloor\frac{k}{2333}\right\rfloor}\mathrm{C}_{n/2333}^t\sum_{i=0}^{2332}\mathrm{C}_{n\%2333}^i+\sum_{i=k-k\%2333}^k\mathrm{C}_{n/2333}^{i/2333}\mathrm{C}_{n\%2333}^{i\%2333}
\]

对于加号后面的式子,\(i/2333=0\),所以是对后面一个式子求和,因此可以用杨辉三角预处理,并求出前缀和,\(O(1)\) 解决。

对于中间一个式子 \(\sum_{i=0}^{2332}\mathrm{C}_{n\%2333}^i\) ,因为 \(n\%2333<2333\) ,同理用杨辉三角。

\[\sum_{i=0}^{2332}\mathrm{C}_{n\%2333}^i=S,\\\ \sum_{i=k-k\%2333}^k\mathrm{C}_{n/2333}^{i/2333}\mathrm{C}_{n\%2333}^{i\%2333}=A
\]

则原式为

\[S\sum_{t=0}^{\left\lfloor\frac{k}{2333}\right\rfloor}\mathrm{C}_{n/2333}^t+A
\]

对于剩下的一个式子,又转化为了一个求和的子问题,因此我们递归解决。递归的边界是 \(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] 超能粒子炮・改 题解的更多相关文章

  1. 洛谷 P4345 [SHOI2015]超能粒子炮·改 解题报告

    P4345 [SHOI2015]超能粒子炮·改 题意 求\(\sum_{i=0}^k\binom{n}{i}\),\(T\)组数据 范围 \(T\le 10^5,n,j\le 10^{18}\) 设\ ...

  2. [洛谷P4345][SHOI2015]超能粒子炮·改

    题目大意:给你$n,k$,求:$$\sum\limits_{i=0}^k\binom n i\pmod{2333}$$题解:令$p=2333,f(n,k)\equiv\sum\limits_{i=0} ...

  3. bzoj4591 / P4345 [SHOI2015]超能粒子炮·改

    P4345 [SHOI2015]超能粒子炮·改 题意:求$\sum_{i=1}^{k}C(n,i)\%(P=2333)$ 肯定要先拆开,不然怎么做呢(大雾) 把$C(n,i)$用$lucas$分解一下 ...

  4. P4345 [SHOI2015]超能粒子炮·改 Lucas

    \(\color{#0066ff}{ 题目描述 }\) 曾经发明了脑洞治疗仪与超能粒子炮的发明家 SHTSC 又公开了他的新发明:超能粒子炮・改--一种可以发射威力更加强大的粒子流的神秘装置. 超能粒 ...

  5. P4345 [SHOI2015]超能粒子炮·改

    传送门 看到数据和模数大小就知道要上 lucas 了 然后开始愉快地推公式: 答案为 $\sum _{i=0}^kC_{n}^{i}\ (mod\ 2333)$ 设 $f [ i ] [ j ] = ...

  6. 【BZOJ4591】[SHOI2015]超能粒子炮·改 (卢卡斯定理)

    [BZOJ4591][SHOI2015]超能粒子炮·改 (卢卡斯定理) 题面 BZOJ 洛谷 题解 感天动地!终于不是拓展卢卡斯了!我看到了一个模数,它是质数!!! 看着这个东西就感觉可以递归处理. ...

  7. Bzoj 4591: [Shoi2015]超能粒子炮·改 数论,Lucas定理,排列组合

    4591: [Shoi2015]超能粒子炮·改 Time Limit: 10 Sec  Memory Limit: 256 MBSubmit: 178  Solved: 70[Submit][Stat ...

  8. bzoj 4591: [Shoi2015]超能粒子炮·改 [lucas定理]

    4591: [Shoi2015]超能粒子炮·改 题意:多组询问,求 \[ S(n, k) = \sum_{i=0}^n \binom{n}{i} \mod 2333,\ k \le n \le 10^ ...

  9. BZOJ_4591_[Shoi2015]超能粒子炮·改_Lucas定理

    BZOJ_4591_[Shoi2015]超能粒子炮·改_Lucas定理 Description 曾经发明了脑洞治疗仪&超能粒子炮的发明家SHTSC又公开了他的新发明:超能粒子炮·改--一种可以 ...

随机推荐

  1. ubuntu16.04 安装openpose

    安装 Anaconda3 Tensorflow-cpu python3tensorflow 1.4.1+opencv3, protobuf, python3-tk ================== ...

  2. 关于sleep的理解

    unix是按时间片轮转调度, windows是抢占式调度 以吃蛋糕为例子,10个人吃蛋糕,如果是unix下, 假设开始时,每个人都处于就绪状态,那么操作系统调度大家排好队,按顺序吃,每个人吃1分钟, ...

  3. SqlServer------范式小结

    说明:大多数初学者对于关系数据库中的范式很是头疼,我本人也是,所以今天又看了视频,总结了一下内容,尽量语言通俗易懂,少用专业术语以及概念. 首先要理解几个键值. 超键:在关系模式中,能唯一标识元组的属 ...

  4. MySql 获取服务提供的sakila数据库(Example Databases)

    关于这个数据库也就是样例数据库,数据库,数据库,最可怕的就是没有数据了,对吧?没有数据你学个什么呀. 可是,没有数据,咱会自己insert,那只能适用于初学者.对于数据库的优化方面的学习,还是有大数据 ...

  5. AIO和NIO的理解

    AIO: AIO 背后的基本思想是允许进程发起很多 I/O 操作,而不用阻塞或等待任何操作完成,可以继续做 另外的事情,等I/O操作完成,内核会通过函数回调或者信号机制通知用户进程.这样很大程度提高了 ...

  6. 不用EL表达式---实现product页面显示

    产品页面显示 静态页面如下: <%@ page language="java" contentType="text/html; charset=UTF-8" ...

  7. POJ3026 Borg Maze(bfs求边+最小生成树)

    Description The Borg is an immensely powerful race of enhanced humanoids from the delta quadrant of ...

  8. Bitmap类、BitmapFactory及BitmapFactory类中的常用方法

    1.Bitmap 1.1非静态方法 public void recycle()——回收位图占用的内存空间,把位图标记为Dead public final boolean isRecycled() —— ...

  9. Extjs Hello extjs

    <html > <head runat="server"> <title></title> <link rel="s ...

  10. C#NPOI.RabbitMQ.EF.Attribute.HttpRuntime.Cache.AD域.List<T>根据指定字段去重.前端JQuery.Cache.I18N(多语言).data-xx(自定义属性)

    使用NPOI 操作Excel 个人使用的电脑基本默认安装Excel 操作起来 调用Excel的组件便可.如果是一台服务器.没有安装Excel,也就无法调用Excel组件. 在此推荐第三方插件.NPOI ...