接着昨天的继续看了一些数论的有关知识。
1.扩展欧几里得算法(已知(a,b),求一组(p,q),使得p*a+q*b=gcd(a,b).)
代码如下:
#include<bits/stdc++.h>
using namespace std;
int extended_gcd(int a,int b,int &x,int &y)
{
int ret,tmp;
if(b==0) //用到了递归当b=0时,q=1,q=0,然后返回去就得到了所求结果。
{
x=1;
y=0;
return a;
}
ret=extended_gcd(b,a%b,x,y);
tmp=x;
x=y;
y=tmp-a/b*y;
return ret;
}
int main()
{
int a,b,x,y,z;
scanf("%d%d",&a,&b);
z=extended_gcd(a,b,x,y);
cout<<z<<" "<<x<<" "<<y<<endl;
}
2.求解线性同余方程
定理1:对于方程a*x+b*y=c该方程等价于a*x三c(mod b),有整数解的充分必要条件是:c%gcd(a,b)=0.
定理2:若gcd(a,b)=1,且x0,y0为a*x+b*y=c的一组解,则该方程的任意解可表示为:x=x0+bt皆成立*t,y=y0-a*t,任意整数t皆成立。
核心代码如下:代码基本都能理解,疑问就是函数的返回值为什么定义成bool类型而不是直接输出x和y;33.
int extended_euclid(int a,int b,int &x,int &y)
{
if(b==0)
{
x=1;
y=0;
return a;
}
int d=extend_euclid(b,a%b,x,y);
int temp=x;
x=y;
y=temp-a/b*y;
return d;
}
bool linearequation(int a,int b,int c,int &x,int &y)
{
int d=extend_euclid(a,b,x,y);
if(c%d)return 0;
int k=c/d;
x*=k;
y*=k;
return 1;
}
3.逆元,
定义:若a*x三1(mod b),a,b互质,则称x为a的逆元。
求你元的代码有点不太理解,
4.中国剩余定理讲解还算大体能懂,但是代码就不怎么能看懂了,最近会仔细研究代码,努力把它弄懂。
5.斐波那契额数,基本的还可以利用递归写出来,但是延伸的东西,就不怎么能看懂,等碰到有关的例题还会回来去看这里。
明天要赶进度了,因为晚来了两天,会尽快赶上去的,感觉数论学起来有点吃力,要多看一看博客上的讲解。