Description
Sheng bill有着惊人的心算能力,甚至能用大脑计算出两个巨大的数的GCD(最大公约
数)!因此他经常和别人比赛计算GCD。有一天Sheng bill很嚣张地找到了你,并要求和你比
赛,但是输给Sheng bill岂不是很丢脸!所以你决定写一个程序来教训他。
Input
共两行:
第一行:一个数A。
第二行:一个数B。
Output
一行,表示A和B的最大公约数。
Sample Input
12
54
Sample Output
6
Data Constraint
Hint
数据范围:
对于20%的数据,1<= A,B<= 10^18。
对于100%的数据,1 <= A,B <= 10^10000。
The Solution
题目大意:高精度的gcd
普通的辗转相除法本来就是一个类似的递归形式,
如果把它打成高精度的话,打起来比较麻烦,且时间复杂度有点悬
其实是我懒得(bu hui)打。。。囧
于是我们可以采用中华民族伟大的智慧结晶——更相减损术
顾名思义,这是一个主要只用除和减的方法求最大公约数
所以它打高精度来,会比较顺手。
具体内容如下:
1、若a为偶数,b为奇数: gcd(a,b)=gcd(a2,b)
2、若a为奇数,b为偶数: gcd(a,b)=gcd(a,b2)
3、若a,b都是偶数: gcd(a,b)=2∗gcd(a2,b2)
3、若a,b都是奇数: gcd(a,b)=gcd(a−b,b)(a>b)
代码实现起来比较恶心
这里提供参考
CODE
#include <cstdio>
#include <iostream>
#include <cmath>
#include <algorithm>
#include <cstring>
#define fo(i,a,b) for (int i=a;i<=b;i++)
#define fd(i,a,b) for (int i=a;i>=b;i--)
#define N 10005
#define yw 1000000000
using namespace std;
char s1[N],s2[N];
int len1,len2,cnt;
struct node
{
int x[1205],l;
}a,b;
bool pd()
{
if (a.l > b.l) return true;
if (a.l < b.l) return false;
fd(i,a.l,1)
if (a.x[i] > b.x[i]) return true;
else if (a.x[i] < b.x[i]) return false;
return true;
}
node Pressure(node a,node b)
{
int k;
node c;
fo(i,1,1200)
{
if (i <= b.l) c.x[i] = a.x[i] - b.x[i];
else if (i <= a.l) c.x[i] = a.x[i];
else c.x[i] = 0;
if (c.x[i] < 0) c.x[i] += yw, a.x[i+1] --;
}
c.l = a.l;
while (c.x[c.l] == 0 && c.l) c.l --;
return c;
}
void Diva()
{
fo(i,1,a.l)
{
if (a.x[i] & 1) a.x[i-1] += yw / 2;
a.x[i] >>= 1;
}
if (! a.x[a.l]) a.l --;
}
void Divb()
{
fo(i,1,b.l)
{
if (b.x[i]&1) b.x[i - 1] += yw / 2;
b.x[i] >>= 1;
}
if (! b.x[b.l]) b.l--;
}
void Mul()
{
fd(i,a.l,1)
{
a.x[i] <<= 1;
a.x[i+1] += a.x[i] / yw;
a.x[i] %= yw;
}
while (a.x[a.l] > 0) a.l ++;
fd(i,b.l,1)
{
b.x[i] <<= 1;
b.x[i+1] += b.x[i] / yw;
b.x[i] %= yw;
}
while (b.x[b.l] > 0) b.l ++;
}
int main()
{
scanf("%s%s",s1 + 1,s2 + 1);
len1 = strlen(s1 + 1);
len2 = strlen(s2 + 1);
if (len1 % 9) a.l = len1 / 9 + 1; else a.l = len1 / 9;
if (len2 % 9) b.l = len2 / 9 + 1; else b.l = len2 / 9;
fo(i,1,a.l)
{
int k1 = max(1,len1 - i * 9 + 1);
int k2 = len1 - (i - 1) * 9;
fo(j,k1,k2) a.x[i] = a.x[i] * 10 + s1[j] - '0';
}
fo(i,1,b.l)
{
int k1 = max(1,len2 - i * 9 + 1);
int k2 = len2 - (i - 1) * 9;
fo(j,k1,k2) b.x[i] = b.x[i] * 10 + s2[j] - '0';
}
while (true)
{
if ((a.x[1] % 2 == 0) && (b.x[1] % 2 == 0)) Diva(),Divb(),cnt ++;
else if (a.x[1] % 2 == 0) Diva();
else if (b.x[1] % 2 == 0) Divb();
if (pd())
{
a = Pressure(a,b);
if (! a.l)
{
while (cnt --) Mul();
while (b.x[b.l] == 0) b.l --;
fd(i,b.l,1)
if (i == b.l) printf("%d",b.x[i]);
else printf("%09d",b.x[i]);
break;
}
}
else
{
b = Pressure(b,a);
if (! b.l)
{
while (cnt --) Mul();
while (a.x[a.l] == 0) a.l --;
fd(i,a.l,1)
if (i == a.l) printf("%d",a.x[i]);
else printf("%09d",a.x[i]);
break;
}
}
}
return 0;
}