![](https://img-blog.csdnimg.cn/20201014180756724.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数论
King yu
KING FOR ME
展开
-
R - Harmonic Number LightOJ - 1234 (分块打表)
了解了分块打表的思想:#include<stdio.h>#include<string.h>#include<math.h>#define ll long long const int maxn=2e6+10;const int maxx=1e8;double a[maxn]={0};int main(){ memset(a,0,sizeof...原创 2019-09-23 20:40:16 · 127 阅读 · 0 评论 -
大数取模
char b[10000];int len=strlen(b); for(int i=0;i<len;i++) { p=(p*10+b[i]-'0')%k; } printf("%lld\n",p);原创 2019-08-26 21:20:54 · 117 阅读 · 0 评论 -
解同余方程组
已知: *N=a1x+b1;N=a2y+b2;。。。。。N=anz+bn;求N;解:a1x+b1=a2y+b2; ==> a1x - a2y=b2-b1; (因为求x 所以a2前正负不影响结果);令c=b2-b1;通过扩展欧几里得算法对此式求解 求出x, 并得到gcd(a1,a2)=d;如果c%d!=0,则此方程组无解;如果有解则:带入到 N=a1...原创 2019-08-13 10:19:03 · 2114 阅读 · 0 评论 -
欧拉函数及其线性筛法 O(n)
欧拉函数:定义phi (n) 为 1 到 n-1 范围内的与互素的数的个数;公式:phi(n)= n*(1- 1/p1) * (1 - 1/p2) * (1 - 1/p3) 。。。。。。(1 - 1/pn);pi为n的素因子;性质:1。 n为1时, phi(1)==1;2。 n为素数时, phi (n)==n-1;3。 若m , n互质 phi(n...原创 2019-08-16 10:56:22 · 254 阅读 · 0 评论 -
求逆元
a*x = 1(mod p); 求 x;当p为素数时 : x==pow( a ,p-2) % p;一般情况:扩展欧几里得:a * x + p * y== 1;// a*x==1(mod b); ==> a*x+b*y==1; 求x; #include<bits/stdc++.h>using namespace std;i...原创 2019-08-06 19:09:54 · 111 阅读 · 0 评论 -
素数表示成两个数平方的和
素数p满足 p恒等于1(mod 4) 则可以表示成两个数的平方和原创 2019-08-06 11:05:00 · 2455 阅读 · 0 评论 -
fmod函数 和 modf函数
fmod函数用来算小数取模:fmod()与求模运算符%的比较:1).%只用于整型的计算,后一个数不能为0;fmod()可以对浮点型数据进行取模运算,后一个数可以为0,返回NaN(NaN,是Not a Number的缩写,用于处理计算中出现的错误情况,比如 0.0 除以 0.0 或者求负数的平方根)。2).%是运算符,fmod()是函数double fmod(double x,double...原创 2019-10-09 20:17:54 · 1070 阅读 · 0 评论