(^-^) --------数论 --------(^-^)
SY_Pistachio
低调的前行,越努力越幸运!
展开
-
B - Number Sequence
A number sequence is defined as follows: f(1) = 1, f(2) = 1, f(n) = (A * f(n - 1) + B * f(n - 2)) mod 7. Given A, B, and n, you are to calculate the value of f(n). Input The i...原创 2018-04-14 16:57:29 · 274 阅读 · 0 评论 -
51Nod 1136 欧拉函数(欧拉函数模板题)
1136 欧拉函数 基准时间限制:1 秒 空间限制:131072 KB 分值: 0 难度:基础题 收藏 关注对正整数n,欧拉函数是少于或等于n的数中与n互质的数的数目。此函数以其首名研究者欧拉命名,它又称为Euler's totient function、φ函数、欧拉商数等。例如:φ(8) = 4(Phi(8) = 4),因为1,3,5,7均和8互质。Input输入一个数...原创 2018-10-11 20:08:39 · 314 阅读 · 0 评论 -
数论 之 筛法总结(艾托拉斯特尼筛法+欧拉筛法)
1.筛法: 2.埃拉托斯特尼筛法(素数/质数筛选法):2.1 步骤: 给出要筛数值的范围n,找出以内的素数。先用2去筛,即把2留下,把2的倍数剔除掉;再用下一个素数,也就是3筛,把3留下,把3的倍数剔除掉;接下去用下一个素数5筛,把5留下,把5的倍数剔除掉;重复进行......详细列出算法如下:将2的倍数(用红色标出),序列变成: 2 3 4 5 6 7 8...原创 2018-10-12 16:11:06 · 1086 阅读 · 0 评论 -
51 Nod 1046 A^B Mod C (数论+矩阵快速幂)
1046 A^B Mod C 基准时间限制:1 秒 空间限制:131072 KB 分值: 0 难度:基础题 收藏 关注给出3个正整数A B C,求A^B Mod C。例如,3 5 8,3^5 Mod 8 = 3。Input3个正整数A B C,中间用空格分隔。(1 <= A,B,C <= 10^9)Output输出计算结果Input示例...原创 2018-10-16 16:55:44 · 220 阅读 · 0 评论 -
51 Nod 1014 X^2 Mod P (数论+二次剩余)
1014 X^2 Mod P 基准时间限制:1 秒 空间限制:131072 KB 分值: 10 难度:2级算法题 收藏 关注X*X mod P = A,其中P为质数。给出P和A,求<=P的所有X。Input两个数P A,中间用空格隔开。(1 <= A < P <= 1000000, P为质数)Output输出符合条件的X,且0 <=...原创 2018-10-16 17:10:11 · 279 阅读 · 0 评论 -
数论 之 中国剩余定理(孙子定理)
1.中国剩余定理理解: 剩余定理,顾名思义就是和余数有关的操作,比如中国剩余定理解决的经典问题: 在《孙子算经》中有这样一个问题:“今有物不知其数,三三数之剩二(除以3 余2),五五数之剩三(除以5 余3),七七数之剩二(除以7 余2),问物几何?” 类似求解这样的问题,便属于中国剩余定理的范畴。但是这样的问题如何来解决呢?在解决本问题之前,我们先来学习两个相关的定理:...原创 2018-10-10 20:26:47 · 1719 阅读 · 0 评论 -
数论 之 乘法逆元介绍(扩展欧几里得算法求解乘法逆元)
1.扩展欧几里得: 已知整数a、b,扩展欧几里得算法可以在求得a、b的最大公约数的同时,能找到整数x、y(其中一个很可能是负数),使它们满贝祖等式:。 2.乘法逆元:ax≡1 (mod p) 这个等式可以描述为: a乘一个数x并模p等于1,即 a%p*x%p=res,res%p=1;看上去就是同余定理的一个简单等式。 逆元x可以描...原创 2018-10-10 21:17:12 · 1653 阅读 · 0 评论 -
中国剩余定理应用 (poj 1006 Biorhythms)
Biorhythms1.中国剩余定理回顾: 这个问题源自于我国数学古书《孙子算经》中的一道问题:“今有物,不知其数,三三数之,剩二;五五数之,剩三;七七数之,剩二。问物几何?”意思是一个整数除以三余二,除以五余三,除以七余二,求这个整数(满足条件且最小)。做法是: n%3==2,n%5==3,n%7==2,并且3,5,7折三个数是互质的。 我们找到使(5*7*a)%...原创 2018-10-23 19:32:23 · 343 阅读 · 0 评论 -
数论——线性筛素数(欧拉筛) 洛谷 3383
P3383 【模板】线性筛素数#include<iostream>using namespace std;const int maxn=1e7+6;int n,m;int flag[maxn];int p[maxn];void getPrime(int n){ int r=0; for(int i=2;i<=n;i++) { if(...原创 2018-10-30 19:27:53 · 282 阅读 · 0 评论 -
洛谷——P1075 质因数分解
P1075 质因数分解质数一个定理:一个数能且只能分解为一组质数的乘积#include<bits/stdc++.h>using namespace std;int main(){ int n; cin>>n; for(int i=2;i<=n;i++) { if(n%i==0) { cout<<n/i<<...原创 2018-10-30 20:11:06 · 1457 阅读 · 0 评论 -
数论学习 之 欧拉函数详解
1.积性函数——欧拉函数——欧拉函数性质——欧拉函数模板:#include<bits/stdc++.h>using namespace std;int Euler(int n){ int ans=n; for(int i=2;i<=sqrt(n);i++) { if(n%i==0) { ans=ans/i*(i-1);//欧拉函数性质三 w...原创 2018-10-11 19:49:27 · 478 阅读 · 0 评论 -
c++中关于M_PI,M_E的解释和使用,以及参考值
1.M_PI 是一个宏定义,圆周率的定义 #define M_PI 3.14159265358979323846 (参考) 此宏定义和编译器有关,TC中M_PI宏就定义在<math.h>里面。 但vc的<math.h>中没有了M_PI的宏定义。2.M_E 是自然对数的一个宏定义 #define M_E 2.7182818284590452353...原创 2018-10-11 18:44:05 · 31187 阅读 · 0 评论 -
51Nod 1130 N的阶乘的长度 V2(斯特林近似) M_PI使用
点击转到1130 N的阶乘的长度 V2(斯特林近似) 基准时间限制:1 秒 空间限制:131072 KB 分值: 0 难度:基础题 收藏 关注输入N求N的阶乘的10进制表示的长度。例如6! = 720,长度为3。Input第1行:一个数T,表示后面用作输入测试的数的数量。(1 <= T <= 1000)第2 - T + 1行:每行1个数N。(1 <...原创 2018-10-11 17:28:47 · 161 阅读 · 0 评论 -
C - u Calculate e
A simple mathematical formula for e is where n is allowed to go to infinity. This can actually yield very accurate approximations of e using relatively small values of n. Output ...原创 2018-04-14 17:25:15 · 307 阅读 · 0 评论 -
欧几里得
什么是欧几里得算法:欧几里得算法就是辗转相除法,用来求得最大公约数。1.1模板一: int gcd(int a,int b){ //return b==0?a:gcd(b,a%b); return b?gcd(b,a%b):a;}1.2完整代码: #include<algorithm> #include<iostream>#include<cstdio&...原创 2018-04-15 09:54:45 · 879 阅读 · 0 评论 -
数论 快速幂
1.快速幂的目的就是做到快速求幂,所以以后遇到要求多次幂指数的时候,可以考虑快速幂算法。2.快速幂是基于二进制的: 假设我们要求a^b,按照朴素算法就是把a连乘b次,这样一来时间复杂度是O(b)也即是O(n)级别,快速幂能做到O(logn),快了好多好多。它的原理如下: 假设我们要求a^b,那么其实b是可以拆成二进制的,该二进制数第i位的权为2^(i-1),例如当b==11时...原创 2018-04-13 20:27:39 · 367 阅读 · 1 评论 -
Parity check
Parity checkTime Limit: 2000 ms Memory Limit: 524288 KiB Submit Statistic Discuss Problem DescriptionFascinated with the computer games, Gabriel even forgets to study. Now sh...原创 2018-04-13 21:09:47 · 211 阅读 · 0 评论 -
Goldbach 米勒 罗宾算法
1.来源:https://blog.csdn.net/qq_38749759/article/details/800425272.算数检测法介绍:https://blog.csdn.net/zengaming/article/details/518672403.米勒 罗宾算法原理介绍:https://blog.csdn.net/fisher_jiang/article/details/98...原创 2018-04-23 21:16:35 · 1354 阅读 · 2 评论 -
数论 斯特林公式 HDU 1018 Big Number
Big NumberTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 41597 Accepted Submission(s): 20377Problem Description In many applications very la...原创 2018-04-29 09:39:44 · 206 阅读 · 0 评论 -
数论 斯特林公式
斯特林公式是一条用来取n阶乘近似值的数学公式。一般来说,当n很大的时候,n阶乘的计算量十分大,所以斯特灵公式十分好用。从图中可以看出,即使在n很小的时候,斯特灵公式的取值已经十分准确。 公式为: 从图中看出,对于足够大的整数n,这两个数互为近似值。更加精确地: ...原创 2018-04-27 17:46:53 · 3314 阅读 · 0 评论 -
51Nod 1256 乘法逆元模板题
点击转到1256 乘法逆元 基准时间限制:1 秒 空间限制:131072 KB 分值: 0 难度:基础题 收藏 关注给出2个数M和N(M < N),且M与N互质,找出一个数K满足0 < K < N且K * M % N = 1,如果有多个满足条件的,输出最小的。Input输入2个数M, N中间用空格分隔(1 <= M < N <= 1...原创 2018-10-11 15:21:36 · 257 阅读 · 0 评论 -
51Nod 1079 中国剩余定理(中国剩余定理模板题)
点击转到1079 中国剩余定理 基准时间限制:1 秒 空间限制:131072 KB 分值: 0 难度:基础题 收藏 关注一个正整数K,给出K Mod 一些质数的结果,求符合条件的最小的K。例如,K % 2 = 1, K % 3 = 2, K % 5 = 3。符合条件的最小的K = 23。Input第1行:1个数N表示后面输入的质数及模的数量。(2 <= N &l...原创 2018-10-11 16:02:13 · 218 阅读 · 0 评论 -
洛谷——P3912 素数个数 Meissel-Lehmer算法 找素数快速算法
P3912 素数个数 线性筛交了一次,一个测试点不过得90分。#include<bits/stdc++.h> using namespace std;const int MAXN=1e8+6; int primes[MAXN],tot=0;bool isPrime[MAXN];int getPrime(int n){ memset(isPrime,true...原创 2018-10-30 20:26:01 · 1892 阅读 · 0 评论