【SDOI2009】Super Gcd （高精度gcd）

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(\dfrac{a}{2},b)$

2、若a为奇数，b为偶数：$gcd(a,b)=gcd(a,\dfrac{b}{2})$

3、若a,b都是偶数：$gcd(a,b)=2*gcd(\dfrac{a}{2},\dfrac{b}{2})$

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;
}