裴蜀定理 对于整系数方程ax+by=m,设d =(a,b) 方程有整数解当且仅当d|m 这个定理实际上在之前学习拓展欧几里得解不定方程的时候就已经运用到 拓展到多元的方程一样适用BZOJ1441 给出n个数(A1...An)现求一组整数序列(X1...Xn)使得S=A1*X1+...An*X...
分类:
其他好文 时间:
2015-04-15 21:01:18
阅读次数:
149
poj 2115 (解单变元模线性方程)
http://poj.org/problem?id=2115
题意:
给出a,b,c,k,求x,使得(a+c*x)%(2^k)=b
限制:
0
思路:
拓展欧几里得单变元模线性方程
令 A=c;C=((b-a)%(2^k)+2^k)%(2^k);B=2^k
则这道题就化为Ax%n=B
对于Ax%B=C
-> Ax...
分类:
其他好文 时间:
2015-01-31 16:30:04
阅读次数:
242
求关于 x 的同余方程 ax ≡ 1 (mod b)的最小正整数解。
最大公约数的求法中最过著名的莫过于欧几里得辗展相除法,它有两种形式(递归与非递归,其实是一样的,任何递归都可以写成非递归)...
分类:
编程语言 时间:
2014-10-31 15:49:21
阅读次数:
373
数论拓展欧几里得,计算mx+ny=d的一组解(m,n为已知)
int xx=x+n/d*i;
int yy=y-m/d*i;//xx,yy分别为其他通解
void extend_gcd(int a,int b,int &x,int &y,int &d)
{
if(!b)
{
x=1;
y=0;
d=a; ...
分类:
其他好文 时间:
2014-10-29 13:04:53
阅读次数:
155
A/BTime Limit: 1000/1000 MS (Java/Others)Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 2390Accepted Submission(s): 1731Problem Descrip...
分类:
编程语言 时间:
2014-10-28 21:23:56
阅读次数:
156
做这道题之前需要先做一下POJ2417,我的题解:http://blog.csdn.net/wyfcyx_forever/article/details/40538515
现在来看这个问题:Ax≡B(mod C)
已知A,B,C
注意这里的A,B,C没有任何限制!
那么考虑我们的传统的GSBS算法为何不能解决这个问题:假设枚举的某个i,我们要利用拓展欧几里得求出存不存在某个A^j(0
那...
分类:
编程语言 时间:
2014-10-28 15:38:14
阅读次数:
260
题目链接:hdu 4828 Grids题目大意:略。解题思路:将上一行看成是入栈,下一行看成是出栈,那么执着的方案就是卡特兰数,用递推的方式求解。#include #include typedef long long ll;const int N = 1000005;const ll MOD = 1...
分类:
其他好文 时间:
2014-10-09 01:40:27
阅读次数:
378
1 void gcd(LL a,LL b,LL &d,LL &x,LL &y) 2 { 3 if (!b) 4 { 5 x=1; 6 y=0; 7 d=a; 8 } 9 else10 {11 g...
分类:
其他好文 时间:
2014-10-07 14:53:03
阅读次数:
149
很容易看出来一个同余式,说到底是解一个线性同余方程,计算机解通常有拓展欧几里得和欧拉定理两种算法,参照去年的NOIP水题,问题是这题数据范围是2^32所以要int64 TAT#include#include#include#includeusing namespace std;__int64 exg...
分类:
其他好文 时间:
2014-10-05 16:15:28
阅读次数:
150
https://vijos.org/p/1781学习了下拓欧。。求exgcd时,因为a*x1+b*y1=a*x2+b*y2=b*x2+(a-b*[a/b])*y2然后移项得a*x1+b*y1=b*x2+a*y2-(b*[a/b])*y2a*(x1-y2)+b*y1-b*x2+(b*[a/b]*y2)...
分类:
其他好文 时间:
2014-10-03 14:03:54
阅读次数:
235