题目链接:https://vjudge.net/problem/UVA-11827
前几天再离散书上看到了欧几里得算法的证明,写一下我的理解:
常见的gcd模板:
long long gcd(long long a,long long b){
long long r=a%b;
while(r){
a=b;
b=r;
r=a%b;
}
return b;
}
递归版本:
long long gcd(long long a,long long b){
return a%b==0?b:gcd(b,a%b);
}
递归版本其实更加直接地指出了欧几里得算法的原理:
如果a%b==0 ,显然gcd(a,b)=b;
否则gcd(a,b)=gcd(b,a%b)
证明:
不妨认为a>b,设a=kb+r
其中r为a除以b的余项,即a%b
所以我们要证明的就是 gcd(a,b)=gcd(b,r)
①令d为a,b任一公因子
②r=a-kb
则依据定理:若d是a和b的因子,则d必然也是xa+yb的因子,x和y为任意整数
所以,d必然也是r的因子
③d是b的因子,又是r的因子,故a,b任一公因子都是b和r的公因子
④同理,可以推出b,r任一公因子都是a和b的公因子
综上,a,b与b,r具有完全相同的公因子,所以它们的最大公因数gcd自然相等
回到这道很恶心的题目:
这道题数据量非常小,直接暴力O()求最大gcd值
坑点在输入,因为没给你输入数的个数,只能自己判断了
个人用的方法:
应该算是比较简洁的了
char s;
long long a[105];
long long t=10,num=0,cnt=0;
while(~scanf("%c",&s)){
if(s>='0'&&s<='9')
num=t*num+s-'0';
else if(num!=0)
a[cnt++]=num,num=0;
else
num=0;
if(s=='\n') break;
}
完整ac代码:
#include <iostream>
#include <cstdio>
#include <algorithm>
#define eps 1e-8
using namespace std;
const long long mod=1e9+7;
long long gcd(long long a,long long b){
return a%b==0?b:gcd(b,a%b);
}
int main()
{
int N;
scanf("%d ",&N);
while(N--){
char s;
long long a[105];
long long t=10,num=0,cnt=0;
while(~scanf("%c",&s)){
if(s>='0'&&s<='9')
num=t*num+s-'0';
else if(num!=0)
a[cnt++]=num,num=0;
else
num=0;
if(s=='\n') break;
}
long long ans=0;
for(int i=1;i<cnt;i++)
for(int j=0;j<i;j++)
ans=max(ans,gcd(a[i],a[j]));
cout<<ans<<endl;
}
return 0;
}