一、.gcd原理及模板
1.gcd利用欧几里得算法(辗转相除法);
2.实现代码:(一般写作一个函数)
非递归
int gcd(int a,int b)
{
int r=a%b;
while®
{ a=b;b=r;r=a%b;}
return b;
}
递归 (常用)
int gcd(int a,int b)
{
return b?gcd(b,a%b):a;
}
二、LCM
lcm 最小公倍数=两数的成积除以两者的最大公约数;
lcm(a,b)=a*b/gcd(a,b)
题目
problem A:最大公约数和最小公倍数
nefu 1077
#include <iostream>
#include <bits/stdc++.h>
using namespace std;
int gcd(int a,int b)
{
return b?gcd(b,a%b):a;
}
int main()
{
int a,b,c,d;
while(scanf("%d%d",&a,&b)!=EOF)
{c=gcd(a,b);
printf("%d ",c);
d=a/gcd(a,b)*b;/*防止数据过大,先除以再乘*/
printf("%d\n",d);
}
return 0;
}
problem B:又见GCD
nefu 992
#include <iostream>
#include <bits/stdc++.h>
using namespace std;
int gcd(int x,int y)
{
return y?gcd(y,x%y):x;
}
int main()
{
int a,b,c;
while(scanf("%d%d",&a,&b)!=EOF)
{
for(c=b*2;gcd(a,c)!=b;c=c+b);/*c肯定是b的倍数*/
printf("%d\n",c);
}
return 0;
}
problem C:多个数的最大公约数
nefu 764
#include <iostream>
#include <bits/stdc++.h>
using namespace std;
int gcd(int x,int y)
{
return y?gcd(y,x%y):x;
}
int main()
{
int n,a[10],s,i;
while(scanf("%d%d",&n,&a[0])!=EOF)
{
s=a[0];
for(i=1;i<n;i++)
{scanf("%d",&a[i]);
s=gcd(s,a[i]);/*求多个数,一个一个接着求*/
}
printf("%d\n",s);
}
return 0;
}
problem D:多个数的最小公倍数
nefu 765
#include <iostream>
#include <bits/stdc++.h>
using namespace std;
long long int gcd(long long int x,long long int y)
{
return y?gcd(y,x%y):x;
}
int main()
{
long long int n,a[10],s,z;
int i;
while(scanf("%lld%lld",&n,&a[0])!=EOF)
{s=z=a[0];
for(i=1;i<n;i++)
{
scanf("%lld",&a[i]);
s=gcd(z,a[i]);/*前两个数的最小公倍数和后一个数找最大公约数再找最小公倍数*/
z=z*a[i]/s;
}
printf("%lld\n",z);
}
return 0;
}
problem E:人见人爱gcd
nefu 1221
由题意的两个条件可得x2+y2=aa-2gcd(a,b)*b;
#include <iostream>
#include <bits/stdc++.h>
using namespace std;
int gcd( int x, int y)
{
return y?gcd(y,x%y):x;
}
int main()
{
int n;
int a,b,t;
while(scanf("%d",&n)!=EOF)
{while(n--)
{
scanf("%d%d",&a,&b);
printf("%d\n",a*a-2*b*gcd(a,b));
}
}
return 0;
}
problem F:LCM&GCD
nefu 1411
#include <iostream>
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;/*表示long long 用ll表示*/
ll gcd(ll x,ll y)
{
return y?gcd(y,x%y):x;
}
ll x,y,t;
int main()
{
ios::sync_with_stdio(false);
cin>>t;
while(t--)
{
cin>>x>>y;
ll flag=0;
for(ll i=x;i<=y;i+=x)
{
if(x*y%i==0)
{
ll j=x*y/i;
/* 不用再判断最小公倍数是否为y,因为j=x*y/i是根据x,y分别为最大公约数
和最小公倍数推导出的,这三个条件满足其中两个,第三个也满足*/
if(gcd(i,j)==x&&j<=y&&j>=x)
flag++;
}
}
printf("%lld\n",flag);
}
return 0;
}
problem G:高木同学的因子
nefu1669
#include <bits/stdc++.h>
using namespace std;
long long gcd(long long x,long long y)
{
if(y==0)return x;
else return gcd(y,x%y);
}
int main()
{
long long x,y,t;
int flag=0,i;
scanf("%lld%lld",&x,&y);
t=gcd(x,y);
for(i=1;i*i<t;i++)/*是最大公约数的因子肯定是两数的公因子*/
{
if(t%i==0)
flag+=2;
}
if(i*i==t)
flag++;
printf("%d\n",flag);
return 0;
}