1, 定义:线性不定方程是指ax + by = c。 2, 给出整数a, b, c,先求出一组解,然后讨论如何表示通解(所有解)。首先设gcd(a; b) = d。如果c不是d的倍数,那么方程无解,因为左边是d的倍数,而右边不是。在有解的情况,设c = c0 £ d,则可以先求出ax + by = ...
分类:
其他好文 时间:
2016-07-19 23:23:18
阅读次数:
194
8-1 buychicken1.c 8-2 buychicken2.c 8-3 bank.c 8-4 ladder.c ...
分类:
其他好文 时间:
2016-07-19 20:51:33
阅读次数:
223
算法竞赛中围绕数论三个重要问题 素数运算 求解不定方程和同余房程 积性函数的应用 1.素数运算实验范例 1).计算[2,n]区间中所有素数 2) 大整数素数测试 1.1 使用筛法生成[2,n]区间中所有素数 a. 最简单之埃拉托斯特尼筛法 思想:辅助数组做筛子,递增搜索筛子中最小数,把被搜索数的倍数... ...
分类:
其他好文 时间:
2016-05-31 15:52:13
阅读次数:
200
扩展欧几里德算法 要是 a*x + b*y = gcd(a,b) ,则一定存在整数对(x0 , y0)使其有解 则 x = x0 + (b/gcd)*t , y = y0 - (a/gcd)*t 是不定方程的通解 因为 a*x + b*y = gcd(a , b); 又因为 gcd(a , b) = ...
分类:
其他好文 时间:
2016-05-04 22:26:10
阅读次数:
143
形如a*x+b*y=c 为不定方程,a,b>0其实无所谓,因为gcd(a,b)=gcd(|a|,|b|) //gcd为最大公约数 由数论的定理所知,当c%gcd==0,不定方程有解,现在我们来求这个解. gcd=gcd(a,b);a*b=gcd*lcm; //lcm为最小公倍数 a'=a/gcd;b ...
分类:
其他好文 时间:
2016-03-29 23:45:18
阅读次数:
278
刚学习的扩展欧几里得算法,刷个水题 求解 线性不定方程 和 模线性方程 求方程 ax+by=c 或 ax≡c (mod b) 的整数解 1、ax+by=gcd(a,b)的一个整数解: <span style="font-size:14px;">void ex_gcd(int a,int b,int ...
分类:
其他好文 时间:
2016-03-29 19:22:43
阅读次数:
178
http://poj.org/problem?id=2115 题解:一个变量从A开始加到B,每次加C并mod2^k,问加多少次。转化为不定方程:C*x+2^K*Y=B-A //poj2115 #include<cstdio> #include<cstdlib> #include<cstring> #
分类:
其他好文 时间:
2016-02-02 21:35:54
阅读次数:
219
http://poj.org/problem?id=1061 题意:裸题。注意负数。 //poj1061 #include<cstdio> #include<cstdlib> #include<cstring> #include<iostream> using namespace std; type
分类:
其他好文 时间:
2016-02-02 21:27:44
阅读次数:
156
今天接到scy的压缩包,开始做数论专题。那今天就总结一下拓展欧几里得求解不定方程和同余方程组。 首先我们复习一下欧几里得算法: 1 int gcd(int a,int b){ 2 if(b==0) return a; 3 return gcd(b,a%b);4 } 拓展欧几里得算法: 推导过程: 给
分类:
其他好文 时间:
2016-02-01 23:51:53
阅读次数:
281