![](https://img-blog.csdnimg.cn/20201014180756927.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数学----数论
pxlsdz
冲就完了!
展开
-
Codeforces 450B Jzzhu and Sequences(矩阵快速幂)
A -Jzzhu and SequencesTime Limit:1000MSMemory Limit:262144KB64bit IO Format:%I64d & %I64uSubmitStatusPracticeCodeForces 450BAppoint description:DescriptionJzzhu has inve...原创 2019-03-28 15:18:17 · 158 阅读 · 0 评论 -
最小质因子
最小质因子欧拉线性筛打表实现原理:欧拉线性筛素数法,每个数仅使用其最小素因数筛去。#include <iostream> using namespace std; inline int mini(int a, int b) {return a < b ? a : b;} // 使用欧拉线性筛打最小质因子表const int LIM = 1e6 + 10...原创 2019-10-25 19:36:30 · 2092 阅读 · 0 评论 -
数学一本通 教堂 数论
教堂 题目描述ROMA城中有一些古典的印度式建筑,这些建筑和周围的欧洲建筑风格格格不入。这些伪装成教堂的建筑其实是某国特工的基地。Tomas接受了一项任务,就是从某个教堂出发,逐个访问这些教堂,搞清楚每一个教堂的内部结构,并回到出发的地方。这些教堂很有规律地构成了一个m*n的矩形,每个教堂和它的八个方向的教堂有直接的路径相连。水平或垂直方向相邻的教堂之间的路程均为1。请问Tomas至...原创 2018-07-30 09:36:17 · 755 阅读 · 0 评论 -
扩展扩展欧几里得算法求逆元
若a*x≡1(mod b) ,a,b互质,则称x为a的逆元。根据逆元的定义,则可以转化为a*x+b*y=1。这样就可以用扩展欧几里得算法求x了。注意:在gcd不为1说明逆元不存在(因为c=1,c%gcd==0为有整数解的充分必要条件),若为1,调整x0到0~m-1的范围中即可int ex_gcd(int a,int b,int &x,int &y) //...原创 2018-07-30 20:30:54 · 8027 阅读 · 0 评论 -
青蛙的约会 扩展的欧几里德算法
青蛙的约会 Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 127383 Accepted: 27617 Description两只青蛙在网上相识了,它们聊得很开心,于是觉得很有必要见一面。它们很高兴地发现它们住在同一条...原创 2018-07-31 11:12:31 · 211 阅读 · 0 评论 -
POJ 2689 Prime Distance 大区间素数筛选
Prime Distance Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 22308 Accepted: 5928 DescriptionThe branch of mathematics called number...原创 2018-08-09 11:28:24 · 174 阅读 · 0 评论 -
POJ 2891 Strange Way to Express Integers 中国剩余定理 不互质
Strange Way to Express Integers Time Limit: 1000MS Memory Limit: 131072K Total Submissions: 19891 Accepted: 67...原创 2018-08-09 11:51:13 · 171 阅读 · 0 评论 -
POJ 2115 C Looooops 扩展欧几里德
C Looooops Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 31879 Accepted: 9269 Descrip...原创 2018-08-09 14:14:29 · 178 阅读 · 0 评论 -
中国剩余定理(CRT) 证明 互质与不互质
中国剩余定理(CRT) 图来自博客:https://blog.csdn.net/u012717411/article/details/43168405我对此做一些解释。首先,(Mi,mi)=1表示它两互质下面这个方程更有利于理解扩展欧几里德如何求解是Mi模mi的逆元,另一种好理解证明:代码模板:typedef long long ll;void ...原创 2018-07-31 22:09:41 · 1279 阅读 · 0 评论 -
POJ Biorhythms 1006 中国剩余定理
Biorhythms Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 145780 Accepted: 46994 DescriptionSome people believe that there are three ...原创 2018-07-31 22:22:43 · 284 阅读 · 0 评论 -
POJ 2407 Relatives 欧拉函数(不打表)
Relatives Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 16563 Accepted: 8410 DescriptionGiven n, a positive integer, how many positi...原创 2018-08-11 08:09:03 · 208 阅读 · 0 评论 -
POJ2487 Farey Sequence 欧拉函数模板题
Farey Sequence Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 18838 Accepted: 7581 Des...原创 2018-08-10 21:46:17 · 208 阅读 · 0 评论 -
线性筛素数 证明详解
线性筛素数代码实现:typedef long long ll; const ll N = 200000; ll prime[N] = {0},num_prime = 0; //prime存放着小于N的素数 int isNotPrime[N] = {1, 1}; // isNotPrime[i]如果i不是素数,则为1 int Prime() ...原创 2018-08-01 08:01:25 · 6166 阅读 · 2 评论 -
十进制小数转换成二进制小数
来自:http://whudongyang.iteye.com/blog/1208120十进制小数转换成二进制小数采用"乘2取整,顺序排列"法。具体做法是:用2乘十进制小数,可以得到积,将积的整数部分取出,再用2乘余下的小数部分,又得到一个积,再将积的整数部分取出,如此进行,直到积中的整数部分为零,或者整数部分为1,此时0或1为二进制的最后一位。或者达到所要求的精度为止。 然后把取...原创 2018-08-10 08:23:39 · 493 阅读 · 0 评论 -
素数 (性质,费马小定理 miller_rabin_素性测试)
转载自Matrix大牛的博客 把代码翻译成C++http://www.matrix67.com/blog/archives/234一个数是素数(也叫质数),当且仅当它的约数只有两个——1和它本身。规定这两个约数不能相同,因此1不是素数。对素数的研究属于数论范畴,你可以看到许多数学家没事就想出一些符合某种性质的素数并称它为某某某素数。整个数论几乎就围绕着整除和素数之类的词转过去转过来。对于写...原创 2018-08-01 11:22:39 · 2036 阅读 · 0 评论 -
计算 n! 中末尾0的个数 求n!中p的重数
题意:输入正整数n,计算 n! 中末尾0的个数 输入:输入一个正整数n (1≤n≤1 000 000 000) 输出:输出 n! 末尾0的个数 样例输入:3 100 1024 样例输出:0 ...原创 2018-08-06 08:49:50 · 414 阅读 · 0 评论 -
NEFU 119 组合素数 (n!素因数p的幂的求法)
组合素数 Problem : 119 Time Limit : 1000ms Memory Limit : 65536K description 小明的爸爸从外面旅游回来给她带来了一个礼物,小明高兴地跑...原创 2018-08-06 09:05:49 · 777 阅读 · 0 评论 -
POJ 3358 Period of an Infinite Binary Expansion 经典数论:欧拉函数+欧拉定理+快速幂取模
Period of an Infinite Binary Expansion Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 3707 ...原创 2018-08-10 10:07:40 · 277 阅读 · 0 评论 -
POJ 2773 Happy 2006 欧拉函数+区间
Happy 2006 Time Limit: 3000MS Memory Limit: 65536K Total Submissions: 14003 ...原创 2018-08-10 11:01:11 · 314 阅读 · 0 评论 -
POJ 1811 Prime Test Pollard_rho算法+Miller_rabin算法
Prime Test Time Limit: 6000MS Memory Limit: 65536K Total Submissions: 37712 Accepted: 10140 Case Time Limit: 4000MS Description Given a big integ...原创 2018-08-01 15:37:57 · 199 阅读 · 0 评论 -
Pollard_rho算法+Miller_rabin算法 大整数的分解
原理证明这个博客写得能看懂:https://www.cnblogs.com/fzl194/p/9047710.html简单例题:POJ 1811 Prime Test这里贴代码很详细的解释,方便套用#include<stdio.h>#include<string.h>#include<stdlib.h>#include<time.h...原创 2018-08-01 15:46:47 · 537 阅读 · 0 评论 -
HDU 5391 Zball in Tina Town 威尔逊定理+快速素数判断法(218ms)
Zball in Tina TownTime Limit: 3000/1500 MS (Java/Others) Memory Limit: 262144/262144 K (Java/Others)Total Submission(s): 2766 Accepted Submission(s): 1295 Problem DescriptionTina Town is...原创 2018-08-06 10:33:43 · 465 阅读 · 0 评论 -
比较高效的素数判定(适用于判断一个数是否为素数)
来自:https://blog.csdn.net/songyunli1111/article/details/78690447标准版:大部分人都知道的比较快的方法:判断从2到sqrt(n)是否存在其约数,时间复杂度O(sqrt(n))高配版:判断2之后,就可以判断从3到sqrt(n)之间的奇数了,无需再判断之间的偶数,时间复杂度O(sqrt(n)/2) 尊享版:首先看一个关于质...原创 2018-08-06 10:38:55 · 348 阅读 · 0 评论 -
数论当中定理
数论当中定理【算术基本定理】定理:每一个大于1的正整数n都可以唯一地写成素数的乘积,在乘积中的素因子按照非降序排列,正整数n的分解式n=(p1^α1)*(p2^α2)*(p3^α3)* ....... *(pk^αk)称为n的标准分解式,其中p1,p2,p3......pk是素数,p1<p2<p3.....,且α1,α2,α3.......是正整数。性质:(1)设d(n)为n的...原创 2018-08-06 10:44:00 · 545 阅读 · 0 评论 -
因子个数 和因子和
【因子个数】实现代码:int count(int n)///求因子的个数{ int s=1;///记录总共的素因子的个数 for(int i=2;i*i<=n;i++) if(n%i==0) { int a=0;///记录的是每个素因子的个数 while(n%i==0) ...原创 2018-08-06 11:45:30 · 1649 阅读 · 0 评论 -
HDU 1215 七夕节 因子和打表
七夕节Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 52255 Accepted Submission(s): 16824 Problem Description七夕节那天,月老来到数字王国,他在城门上贴了一张告示,...原创 2018-08-06 11:48:34 · 269 阅读 · 0 评论 -
POJ 2773 Happy 2006 二分+容斥原理(二进制枚举或dfs)
Happy 2006 Time Limit: 3000MS Memory Limit: 65536K Total Submissions: 14003 ...原创 2018-08-10 16:13:43 · 321 阅读 · 0 评论 -
欧拉函数的定义与计算
欧拉函数 【定义】数论,对正整数n,欧拉函数是小于或等于n的数中与n互质的数的数目。此函数以其首名研究者欧拉命名,它又称为Euler's totient function、φ函数、欧拉商数等。 例如φ(8)=4,因为1,3,5,7均和8互质。 从欧拉函数引伸出来在环论方面的事实和拉格朗日定理构成了欧拉定理的证明...原创 2018-08-01 20:58:56 · 10441 阅读 · 0 评论 -
POJ 2034 Anti-prime Sequences dfs
Anti-prime Sequences Time Limit: 3000MS Memory Limit: 30000K Total Submissions: 4117 Accepted: 1853 ...原创 2018-08-10 20:18:53 · 313 阅读 · 0 评论 -
HDU 3579 Hello Kiki 中国剩余定理(不互质情况)
Hello KikiTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 5099 Accepted Submission(s): 2009 Problem DescriptionOne day I was shopping...原创 2018-08-06 19:38:37 · 287 阅读 · 0 评论 -
逆元方法总结
逆元(inv)1.逆元若a*x≡1(mod b) ,a,b互质,则称x为a的逆元。A/B %p=A*inv(B) %p inv(B)为B的逆元2.应用(a/b)%m 时,因b可能会过大,会出现爆精度的情况,且除法不满足同余拆分定理。3.求逆元的方法(1).费马小定理在m是素...原创 2018-08-02 00:08:46 · 2646 阅读 · 0 评论 -
最大公约数(欧几里得算法和二进制算法)
欧几里得算法原理欧几里德算法又称辗转相除法,用于计算两个整数a,b的最大公约数。 gcd(a,b)=gcd(b,a mod b);证明过程证明一:证明:a可以表示成a = kb + r,则r = a mod b假设d是a,b的一个公约数,则有d|a, d|b,而r = a - kb,因此d|r因此d是(b,a mod b)的公约数假设d 是(b,a mod b)...原创 2018-08-02 00:23:46 · 2007 阅读 · 0 评论 -
数论重要内容总结
数论重要内容快速幂与快速乘 二进制算法:点这里矩阵快速幂之矩阵加速中矩阵构造:点这里最大公约数和最小公倍数(欧几里德算法和二进制算法):点这里线性同余方程(扩展欧几里得算法):点这里同余方程组(中国剩余定理):点这里高斯消元:点这里数论一些容易忽略的定理:点这里素数的判定1)高效小范围素数判定(<1e9):点这里2) Miller—Rabin素数测定(...原创 2018-08-02 00:30:38 · 835 阅读 · 0 评论 -
POJ 1595 Prime Cuts 素数打表+找中间数规律
Prime Cuts Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 12015 Accepted: 4569 DescriptionA prime number is a counting number (1, 2, ...原创 2018-08-11 10:05:49 · 216 阅读 · 0 评论 -
POJ 3641 Pseudoprime numbers
Pseudoprime numbers Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 12461 Accepted: 5367 DescriptionFermat's theorem states that for a...原创 2018-08-11 10:50:01 · 205 阅读 · 0 评论 -
poj 2191Mersenne Composite Numbers 大整数因式分解+素数判断
Mersenne Composite Numbers Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 2471 Accepted: 1147 DescriptionOne of the world-wide cooper...原创 2018-08-11 15:23:49 · 465 阅读 · 0 评论 -
POJ 2429 GCD & LCM Inverse(素数判定Miller-Rabin+素因子分解Pollard-rho+二进制枚举)
GCD & LCM Inverse Time Limit: 2000MS Memory Limit: 65536K Total Submissions: 18825 Accepted: 3486 DescriptionGiven two positive integers a an...原创 2018-08-07 09:29:19 · 282 阅读 · 0 评论 -
POJ poj 2142 The Balance 扩展欧几里得 |x|+|y|最小
The Balance Time Limit: 5000MS Memory Limit: 65536K Total Submissions: 8784 Accepted: 3817 DescriptionMs. Iyo Kiffa-Australis has a balance and o...原创 2018-08-11 16:28:09 · 288 阅读 · 1 评论 -
poj 2262 Goldbach's Conjecture 素数打表
Goldbach's Conjecture Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 48891 Accepted: 18614 DescriptionIn 1742, Christian Goldbach, a ...原创 2018-08-11 17:12:07 · 162 阅读 · 0 评论 -
POJ 2909 Goldbach's Conjecture
Goldbach's Conjecture Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 11574 Accepted: 6807 DescriptionFor any even number n greater th...原创 2018-08-11 17:26:26 · 163 阅读 · 0 评论