![](https://img-blog.csdnimg.cn/20201014180756930.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数论、数学
文章平均质量分 72
bubbleoooooo
这个作者很懒,什么都没留下…
展开
-
同余问题之中国剩余定理 _(:з」∠)_
中国剩余定理:中国剩余定理又指孙子定理。叙述如下:令m1,…,mr两两互素, b1,…,br为整数, x≡b1(mod m1),x≡b2(mod m2)……x≡br(mod mr);有唯一正整数解x,其形式为:x=∑biMi’Mi(mod M)(1M=∏mi(1符号定义完了。下面开始证明。我们需要证明两个问题。x是方程组的解(负数+M取模为正就好)。(2)x是唯一解。原创 2014-08-04 21:41:30 · 1022 阅读 · 0 评论 -
hdu 4818 RP problem
//高斯消元法解线性方程组//每一个方程都是 进==出//Ai1 * x1 + Ai2 *x2 + Ai3 * x3 +...+ An-1 * xn-1 = xi;//所以别的点流进==本点流出xi//移项发现每个A[i][i]=-1;//A[i][j]=1/d[j];d[j]为j点出度。//n个方程解n个未知数。要么是唯一解,要么是无穷解。//若是唯一解,则枚举所有A->j(j为与原创 2015-07-20 12:05:28 · 978 阅读 · 0 评论 -
CF571A - Lengthening Sticks(容斥)
//571A - Lengthening Sticks//数学题,容斥原理//总的选择数C(l+3,3) la+lb+lc+unused=l;//确定最长边,酱紫不构成三角形的充分必要条件就为 最长边>=其他两边和//枚举最长边增量i,其他两边为b+lb,c+lc//a+la>=b+lb+c+lc lb+lc<=l-la;//x=min(a+la-b-c,l-la);//则不构原创 2015-08-25 09:23:41 · 512 阅读 · 0 评论 -
codeforces 559C Gerald and Giant Chess(组合数学)
//CF559C//1.从(0,0)到(x,y)的方法数为C(x+y,x);//2.转换成总的种数-k个禁走点的方法数*k的禁走点到终点的方法数。//这里的k个禁走点的方法数不经过之前的禁走点,所以相互独立//3.大组合数取模公式 C(x+y,x)=(x+y)!/(x!*y!)//4 (a/b)%mod = a*b^(mod-2)%mod;#include#include#incl原创 2015-08-09 18:25:12 · 527 阅读 · 0 评论 -
poj 1845 Sumdiv
知识点:1推论5 (除数和函数)设a是正整数,σ(a)表示a的所有正除数之和,那么σ(1)=1,当a有标准素因数分解式时,σ(a)=∏((pj^αj+1)-1/(pj – 1))(12(a + b) % p = (a % p + b % p) % p(a * b) % p = (a % p * b % p) % p3取模运算中不能出现除法,等比数列求和公式不成立原创 2015-03-11 16:08:23 · 580 阅读 · 0 评论 -
算数基本定理
算数基本定理 符号:a|b a整除b(a,b) a,b最大公约数定理1:设p是素数, p|a1a2,那么p|a1或p|a2至少有一个成立。一般的,若p|a1…ak, 则p|a1,…,p|ak至少有一个成立。 定理2:设a>1, 那么必有a=p1p2…ps;其中pj(1将相同的素数合并,即得a=(p1^a1)***(ps^as)(p1标准素因数分解式)(与原创 2015-03-11 12:00:10 · 1083 阅读 · 0 评论 -
容斥原理之素数部分模板
二进制:void getfactor(__int64 num){ fac.clear(); for(int i=2;i*i<=num;i++) { if(num%i==0) { fac.push_back(i); while(num%i==0)num/=i; } }原创 2014-11-01 14:25:57 · 836 阅读 · 0 评论 -
欧拉phi函数与欧拉定理
欧拉phi函数欧拉phi(x)函数等于不超过x且和x互素的整数个数。欧拉phi函数证明:1. p为质数时。2. 设n =,则n前有n-1,即个整数。因为n的因子只有1,p,n。只需要考虑由p生成的小于n的合数个数,即{p,2*p, … ,()*p}共有 个。所以.3. 设n=p*q. 则n前有n-1,即pq-1个整数。因为n的因子只有1,p , q , pq原创 2014-08-23 10:52:46 · 1696 阅读 · 0 评论 -
最大公约数gcd算法及其扩展
gcd(greatest common divisor)模板如下:简单写下证明:求证gcd(a,b)=gcd(b,a%b)令a=k*b+mod;x=gcd(a,b),y=gcd(b,a%b)=gcd(b,mod);x是a,b的最大公约数,显然可以整除a,b,由于mod=a-k*b,x可以整除mod。x是b,mod公约数。y是b,mod的最大公约数,显然可以整原创 2014-08-23 13:03:17 · 1215 阅读 · 0 评论 -
java 大数用法 处理高精度
在用C或者C++处理大数时感觉非常麻烦,但是在JAVA中有两个类BigInteger和BigDecimal分别表示大整数类和大浮点数类,至于两个类的对象能表示最大范围不清楚,理论上能够表示无线大的数,只要计算机内存足够大。这两个类都在java.math.*包中,因此每次必须在开头处引用该包。Ⅰ基本函数:1.valueOf(parament); 将参数转换为制定的类型转载 2014-08-16 18:57:32 · 697 阅读 · 0 评论 -
排列组合之生成组合_(:з」∠)_
1.生成集合的所有子集原创 2014-08-19 12:39:20 · 884 阅读 · 1 评论 -
同余问题之模方程模板
线性同余方程原创 2014-08-23 17:16:53 · 629 阅读 · 0 评论 -
排列组合之生成排列_(:з」∠)_
1.生成全排列第一类:按字典序shengcheng原创 2014-08-19 11:42:59 · 985 阅读 · 0 评论 -
高斯消元模板整理
//高斯消元法解异或方程组,返回方程解得个数。const int N = 30;int A[N][N];//关系矩阵int Gauss(int equ,int var){//返回解得个数。 int row,col; for(row=0,col=0;row<equ&&col<var;col++,row++){ int max_r=row;//默认最大为本行原创 2015-07-20 12:18:47 · 2396 阅读 · 2 评论