快速乘
Facico
生命是闪耀的此刻,不是过程,就像芳香不需要道路一样。
展开
-
WYF的盒子
Description 给出n,m,k,p,求 ∑i=nmikmodp\sum_{i=n}^mi^k\mod p Solution 首先,自然数幂求和,好多方法啊。解决自然数幂和的各种方法 由于第一类斯特林数不用中国剩余定理来对于mod操作特殊处理,所以这个好。原创 2016-03-15 22:05:52 · 1648 阅读 · 0 评论 -
【NOI2012】随机数生成器
Description给你一个式子 Xn+1=(a∗Xn+c)modmX_{n+1}=(a*X_n+c) \mod m 求第n项Solution用什么发现n十分的大101810^{18},又只给1000ms,只能用O(logn)O(\log n)的方法 有什么可以做呢? 很明显是矩阵乘法!!!构造一个转移矩阵。先要弄初始矩阵。 先设初始矩阵有一位XnX_n,首先需要乘,并不用多开一位。要从原创 2016-03-10 18:41:09 · 1298 阅读 · 0 评论 -
处理两数相乘的黑科技
黑科技来源我们发现如果两个数都接近long long,乘起来就爆掉了,那该怎么办。黑科技快速乘简介这个应该是白科技。 虽然叫快速乘,但是十分的慢。 速度比快速幂多一个log 。 比如说是a*b,我们设f(i)=a*i 那么把b拆一下,f(b)=f(b/2) * 2+a *(b mod 2) 递归的做Codell qsc(ll x,ll y){ ll z=0; if(y==原创 2016-03-16 14:21:14 · 1621 阅读 · 0 评论 -
【51Nod 1691】比大小
Description有两个数列A和B有两个数列A和B 已知A0,a,b,N已知A_0,a,b,N An=A(n−1)∗a+b(n>=1)A_n=A_(n-1)*a+b (n>=1) B数列满足B数列满足 Bn=2∗B(n/2)+1(n为偶数)B_n=2*B_(n/2) + 1 (n为偶数) Bn=2∗B((n−1)/2)+(n+1)/2(n为奇数)B_n=2*B_((n-1)/2) +原创 2016-07-06 10:34:45 · 741 阅读 · 0 评论 -
【NOIP模拟】快速荷叶叶变换
DescriptionSolution以前做过一道加强版的题目,叫做模积和。 所以第一题,秒切。Code#include<iostream> #include<cstdio> #include<cstring> #include<cmath> #include<algorithm> #define fo(i,a,b) for(i=a;i<=b;i++) using namespace std; c原创 2016-09-04 12:27:59 · 628 阅读 · 0 评论