2016年09月_不归路___

原创素数判定

一个恒成立的等式a^(p-1)%p≡1(p为素数 a 但如果反过来（倘若a^(p-1)%p=1,a 但是，它的否逆定理成立。即若a^(p-1)%p≠1，则，p为合数。所以，对于需要判断素数的p，我们随机枚举若干个a（a 显然a^(p-1)需要用快速幂，虽然可以边乘边模，但是a*a显然也有可能超。所以定一个函数f(a,b,p)来计算a*b%p; 如果b为偶数 f(a,b,p)=f((

2016-09-14 22:12:38 1215

原创斐波那契数列矩阵快速幂

#include #include using namespace std; const int init[2][2]={{0,1},{1,1}}; struct node{ int x[2][2]; node(){ memset(x,0,sizeof(x)); } }; node ci(node a){ node turn; int i,j,k; for(i=0;i<2;i++)

2016-09-12 21:56:07 987

原创树状数组

树状数组满足树形结构，只需要n的数组。上图更直观：上面的树状数组是c，下面队列数组是a。 c[2]是c[1]和a[2]的父亲节点（也可以理解为c[1]=a[1],c[2]是a[1]和a[2]的父亲节点，注意不是所有的i都有a[i]=c[i],a[2]≠c[2]，显然只有a[2*k+1]=c[2*k+1]）c[6]是a[5](c[5])和a[6]的父亲节点，c[4]是c[2](

2016-09-05 21:51:55 736

原创欧拉函数

欧拉函数求的是小于等于n的正整数中与n互质的数的个数。例如φ(8)=4，因为1,3,5,7均和8互质。通式：其中p1, p2……pn为x的所有质因数，x是不为0的整数。 φ(1)=1（唯一和1互质的数(小于等于1)就是1本身）。注意：每种质因数只一个。比如12=2*2*3那么φ（12）=12*（1-1/2）*(1-1/3)=4 若n是质数p的k次幂，，因为

2016-09-03 16:43:35 876

转载中国剩余定理

中国剩余定理（CRT）的表述如下设正整数两两互素，则同余方程组有整数解。并且在模下的解是唯一的，解为其中，而为模的逆元。代码： int CRT(int a[],int m[],int n) {

2016-09-03 15:37:47 856 1

转载求逆元的方法汇总

逆元它是一个可以取消另一给定元素运算的元素。对于正整数a和m，如果有a⋅x≡1(mod m)，那么把这个同余方程中x的最小正整数解叫做a模m的逆元。前提求a(modm)意义下的逆元,要求a与m互质,否则不存在乘法逆元定理欧拉定理(费马小定理) 欧拉定理：若a与m互质，那么有aφ(m)≡1(mod m) 方法 1.欧拉定理根据欧拉定理 aϕ

2016-09-03 15:06:38 2240

原创拓扑排序

拓扑序列是将一个有向无环图的点排序后的序列，排序规则是，对任意一条由v指向u的边都满足在序列中v在u左边，显然，对于一个有向无环图，拓扑序列可能不唯一。为什么要有向？因为拓扑排序的规则就是方向。为什么要无环？如果对于u和v两个点在同一个环内，就意味着u可以到v，v可以到u，那么在拓扑序列中，u是在v的左边还是右边呢？说一下拓扑排序的思路。（看了很多代码都没看懂只看懂了拓扑序列下面是

2016-09-03 14:46:32 866

My博客