【BZOJ-4591】超能粒子炮·改 数论 + 组合数 + Lucas定理
2016-05-18 21:23
537 查看
4591: [Shoi2015]超能粒子炮·改
Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 95 Solved: 33
[Submit][Status][Discuss]
Description
曾经发明了脑洞治疗仪&超能粒子炮的发明家SHTSC又公开了他的新发明:超能粒子炮·改--一种可以发射威力更加强大的粒子流的神秘装置。超能粒子炮·改相比超能粒子炮,在威力上有了本质的提升。它有三个参数n,k。它会
向编号为0到k的位置发射威力为C(n,k) mod 2333的粒子流。现在SHTSC给出了他的超能粒子炮·改的参数,让你求
其发射的粒子流的威力之和模2333。
Input
第一行一个整数t。表示数据组数。之后t行,每行二个整数n,k。含义如题面描述。
k<=n<=10^18,t<=10^5
Output
t行每行一个整数,表示其粒子流的威力之和模2333的值。Sample Input
15 5
Sample Output
32HINT
Source
By 佚名上传Solution
Lucas定理算是裸题?大概就是预处理出组合数和前缀和,然后Lucas搞搞...
Code
#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> #include<cmath> using namespace std; long long read() { long long x=0,f=1; char ch=getchar(); while (ch<'0' || ch>'9') {if (ch=='-') f=-1; ch=getchar();} while (ch>='0' && ch<='9') {x=x*10+ch-'0'; ch=getchar();} return x*f; } #define mod 2333 #define maxn 2500 int T;long long N,K; int C[maxn][maxn],Sum[maxn][maxn]; void GetC(int n) { C[0][0]=1; for (int i=1; i<=n; i++) { C[i][0]=1; 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<=n; i++) { Sum[i][0]=C[i][0]; for (int j=1; j<=n; j++) Sum[i][j]=(Sum[i][j-1]+C[i][j])%mod; } } int Lucas(long long n,long long m) { if (!m) return 1; return C[n%mod][m%mod]*Lucas(n/mod,m/mod)%mod; } int Calc(long long n,long long k) { if (k<0) return 0; return ((Calc(n/mod,k/mod-1)*Sum[n%mod][mod-1])%mod+(Lucas(n/mod,k/mod)*Sum[n%mod][k%mod])%mod)%mod; } int main() { T=read(); GetC(2333); while (T--) { N=read(),K=read(); printf("%d\n",Calc(N,K)); } return 0; }
相关文章推荐
- 第十二周技术博客
- PHP AJAX返回 "TEXT"
- 解码H264文件的一些基础知识
- 极角排序
- Java文件流 字节流和字符流
- 深入浅出: Java回调机制(异步)
- 我的Unity3D学习日记-01
- JAVA获取本机IP地址
- 从头认识多线程-2.23 静态内部类的同步机制跟普通类相同
- Try Git 译文
- java中重载和重写
- 矩阵快速幂+map 雷神之路 scu dp练习A题
- POJ1426 建议打表做该题/用数组模拟队列(比STL队列速度要快)
- Yii2.0 自定义日志类
- Libgdx之Label Image
- 面向项目(十一)—— 库的使用
- JAVA-007前面所学总结加深
- 第十三周的学习进度表
- IOI2000 邮局
- iOS蓝牙开发(一)蓝牙相关基础知识