/* 快速乘法 */ #include #include void read (long long &now) { register char word = getchar (); for (now = 0; !isdigit (word); word = getchar ()); for (; i... ...
分类：
其他好文 时间：
2017-08-08 23:00:05
阅读次数：
382
#include using namespace std; const int N = 1e6+5; char a[N]; int main() { scanf("%s",a+1); int n = strlen(a+1); int l = 1,r = n; while(l = 0 && a[r] ... ...
分类：
其他好文 时间：
2017-07-26 01:58:40
阅读次数：
159
仓鼠大神博客：http://www.cnblogs.com/linyujun/p/5194184.html http://www.cnblogs.com/linyujun/p/5194170.html ①gcd 一行版： 位运算版： ②快速幂和快速乘 ③扩展欧几里得（扩展gcd） ④乘法逆元 扩展欧 ...
分类：
编程语言 时间：
2017-07-14 23:58:51
阅读次数：
389
说到素数不得不说素数判定算法。 其中极为经典的为Rabin Miller测试。 通过二次探测的方法，我们可以将其正确率上升到一个很高的高度。 二次探测的原理我还是不太懂，所以NOI前我暂时只是梳理一下这个算法的流程。 首先，我来介绍一些小Trick。 $O(1)$的快速乘。 在一些卡常数而且爆lon ...
分类：
其他好文 时间：
2017-07-10 16:28:08
阅读次数：
228
题目背景 矩阵快速幂 题目描述 给定n*n的矩阵A，求A^k 输入输出格式 输入格式： 第一行，n,k 第2至n+1行，每行n个数，第i+1行第j个数表示矩阵第i行第j列的元素 输出格式： 输出A^k 共n行，每行n个数，第i行第j个数表示矩阵第i行第j列的元素，每个元素模10^9+7 输入输出样例 ...
分类：
其他好文 时间：
2017-07-07 22:33:27
阅读次数：
137
时间限制:20000ms 单点时限:1000ms 内存限制:256MB 描述 在写代码时，我们经常要用到类似 x × a 这样的语句( a 是常数)。众所周知，计算机进行乘法运算是非常慢的，所以我们需要用一些加法、减法和左移的组合来实现乘一个常数这个操作。具体来讲， 我们要把 x × a 替换成：( ...
分类：
其他好文 时间：
2017-06-27 12:56:29
阅读次数：
153
题意：c( n, m)%M M = P1 * P2 * ......* Pk （其中Pk是素数) 思路：Lucas定理中C(n,m)%M，M必须是素数，当M不是素数时，我们可以把它拆成素数的乘积 如果x=C(n,m)%M ,M=p1*p2*..*pk; a[i]=Lucas(n,m)%pi； xΞa ...
分类：
其他好文 时间：
2017-06-13 12:44:44
阅读次数：
231
Given the number, you are to answer the question: "Is it prime?"Solutions to this problem can be submitted in C, C++, Pascal, Perl, Python, Ruby, Lisp ...
分类：
其他好文 时间：
2017-06-03 17:19:39
阅读次数：
275
