qeroLai的博客_ACM之路--数论一

大整数分解pollard（注意2^59次方分解不了）

void Pollard(long long n); void Factor(long long n) { long long d = 2; while (true) { if (n % d == 0) { Pollard(d); Pollard(n / d); return; } d++; } } void Pollard(long long n) { if (n <= 0) printf

return

原创 2014-02-13 17:43:39 1382 阅读

poj3518筛法

/*题意：给出一个数，若为合数则求出其最靠近的两边的质数之差（距离），若是质数，则输出0即可解题：阿拉斯托散筛法的应用，试模板的题目 63msACc++代码*/ #include <iostream> #include <cstdio> #include <cstring> #include <cmath> #include<algorith

include

visit

原创 2014-02-13 17:43:03 588 阅读

poj3641大素数判断（拉宾米勒）和快速模平方（蒙哥马利）

/*题意：给出p和a，如果p不是素数且a^p === a（mod p）的话，输出yes，否者输出no 解题：蒙哥马利模平方计算+拉宾米勒素数判定方法 0msACc++代码*/ #include <iostream> #include <cstdio> #include <cstring> #include <cmath> #include<al

include

蒙哥马利

原创 2014-02-13 17:42:29 1216 阅读

poj1730Pollard分解

/*题意：给出一个数，注意有可能是负数，这个数能表示为x=b^p,求最大的p 解题：这道题比较卡精度，所以直接用__int64，负数的话就输入后处理一下，用Pollard分解一下这个数，获得这个数得标准分解式，即这个数的素数乘积表达式，问题就是求组成这个数的这些素数的幂数的最大公约数，注意负数的答案只能是奇数，若得到的答案是偶数，需要不断除以2直至获得一个奇数答案为止。 16msACc++代码

a

原创 2014-02-13 17:41:53 619 阅读

拉宾米勒素数判断

int_type Montgomery(int_type n,int_type p,int_type m) { //蒙哥马利法快速计算(n^p)%m的值 int_type k = 1; n%=m; while(p!=1) { if(0!=(p&1))k=(k*n)%m; n=(n*n)%m; p>>=1; } return(n*k)%m; } bool Miller_rabin

return

false

蒙哥马利

原创 2014-02-13 17:41:11 1040 阅读

hdu1573非互素的中国剩余定理（疑似有错）

/*大牛的思路：典型的中国剩余定理，但是这里是非互质情况下的中国剩余定理。解题思路： 1.因为(a1,a2,a3,a4,….,ak)不一定互质，所以不能够直接用中国剩余定理。 2.x=r1+a1*k1,x=r2+a2*k2,所以有r1+a1*k1=r2+a2*k2，化简后得到 a1*k1=(r2-r1) mod(a2); 用扩展欧几里得可以得到最小的k1,所以x=r1+a1*k1+a1*a2/

中国

原创 2014-02-13 17:40:29 919 阅读

hdu1452整数因子和与模平方计算方法

#include <iostream> #include <cstdio> #include <cstring> #include <cmath> using namespace std; int GetRemainder(int baseNum, int power, int modelNum) { //判断各个值正确性 baseNum %= mo

include

power

计算方法

原创 2014-02-13 17:39:49 698 阅读

hdu1212大数求余

/*解题思路：分解为十进制的一位，每位求余一次。即设str[0--n]为大数(0为个位，1为十位，2为百位，...) 有 ans[0] = str[0] % b； ans[1] = str[1] * 10 % b; . . . ans[n] = str[n]*10^n % b; 这里要解决对10的n次方求余的问题，其实很简单，设temp[k]=10^k%b,必有 temp[k] = temp[k]

十进制

原创 2014-02-13 17:38:55 633 阅读

poj3517数学方法解约瑟环

#include <iostream> #include <cstdio> #include <cstring> #include <cmath> using namespace std; int n,k,m; int main() { while(scanf("%d %d %d",&n,&k,&m),n||m||k) { i

include

数学

return

原创 2014-02-13 17:38:11 673 阅读

poj1781找规律约瑟环

/*这题一开始用约瑟环的递推公式做，发现无限的超时，果断打表一看，前五十如下： 1 1 2 1 3 3 4 1 5 3 6 5 7 7 8 1 9 3 10 5 11 7 12 9 13 11 14 13 15 15 16 1 17 3 18 5 19 7 ... 规律出来了，先用等比求和算出在第几个循环里面，再用等差数列递推式算出是该循环的第几个值*/ #include <iostream

a

原创 2014-02-13 17:37:42 841 阅读

官方博客	全部文章	热门标签	班级博客
了解我们	网站地图	意见反馈

鸿蒙开发者社区	51CTO学堂
51CTO	软考资讯

51CTO博客

Qero

大整数分解pollard（注意2^59次方分解不了）

poj3518筛法

poj3641大素数判断（拉宾米勒）和快速模平方（蒙哥马利）

poj1730Pollard分解

拉宾米勒素数判断

hdu1573非互素的中国剩余定理（疑似有错）

hdu1452整数因子和与模平方计算方法

hdu1212大数求余

poj3517数学方法解约瑟环

poj1781找规律约瑟环