[SDOI2009]SuperGCD 高精度GCD

Description

Sheng bill有着惊人的心算能力，甚至能用大脑计算出两个巨大的数的GCD（最大公约 数）！因此他经常和别人比

赛计算GCD。有一天Sheng bill很嚣张地找到了你，并要求和你比 赛，但是输给Sheng bill岂不是很丢脸！所以你

决定写一个程序来教训他。

Input

共两行： 第一行：一个数A。 第二行：一个数B。

0 < A , B ≤ 10 ^ 10000。

Output

一行，表示A和B的最大公约数。

Sample Input

12
54

Sample Output

6

solution

---

用更相减损术来计算GCD

优化：

1.

对于lena>lenb的情况

a=a-b*(lena-lenb-1)  (具体实现 只需要标记一下，见代码)

2.

把2都提出来  (但是好像并没有快)

3.

如果重载运算符，取地址可以变快...

---

这题最坑的是 phthon ...

 

 

#include<cstdio>
#include<cstring>
#include<iostream>
#define mem(a,b) memset(a,b,sizeof(a))
#define ll long long
using namespace std;
const int Maxn=10066;

int comp;

struct bign
{
    int s[Maxn],len;
    bign(){len=1;mem(s,0);}
    __attribute__((optimize("-O3"))) inline bign operator * (const int &c) const
    {
        bign now=*this,temp;
        temp.len=now.len+15;
        int x;
        for(int i=1;i<=temp.len;++i)
        {
            temp.s[i]+=now.s[i]*c;
            temp.s[i+1]+=temp.s[i]/10;
            temp.s[i]%=10;
        }
        while(temp.s[temp.len]==0&&temp.len>1)
          --temp.len;
        return temp;
    }
    __attribute__((optimize("-O3"))) inline bign operator - (const bign &c) const//取地址会变快... 
    {
        bign now=*this;
        int q1=(now.len>c.len?now.len:c.len);
        for(int i=1+comp;i<=q1;++i)
        {
            now.s[i]-=c.s[i-comp];
            if(now.s[i]<0)
            {
                --now.s[i+1];
                now.s[i]+=10;
            }
        }
        while(now.len>1&&now.s[now.len]==0)
          --now.len;
        return now;
    }
    __attribute__((optimize("-O3"))) inline bign operator / (const int &c) const
    {
        bign now=*this,temp;
        int x=0;
        for(int i=now.len;i>=1;--i)
        {
            temp.s[i]=(x*10+now.s[i])/c;
            x=(x*10+now.s[i])%c;
        }
        temp.len=now.len;
        while(temp.len>1&&temp.s[temp.len]==0)
          --temp.len;
        return temp;
    }
};

bign a,b;
char s[10066];
int ji;

__attribute__((optimize("-O3"))) inline void read()
{
    scanf("%s",s);
    int len=strlen(s);
    a.len=len;
    for(int i=0;i<len;++i)
      a.s[len-i]=s[i]-'0';
    scanf("%s",s);
    len=strlen(s);
    b.len=len;
    for(int i=0;i<len;++i)
      b.s[len-i]=s[i]-'0';
}

__attribute__((optimize("-O3"))) inline void out11(bign s)
{
    for(int i=s.len;i>=1;--i)
      printf("%d",s.s[i]);
    printf("\n");
}

__attribute__((optimize("-O3"))) inline int com()
{
    if(a.len>b.len)
      return 3;
    if(b.len>a.len)
      return 1;
    int flag;
    for(int i=(a.len>b.len?a.len:b.len);i>=1;--i)
    {
        if(a.s[i]==b.s[i])
          continue;
        if(a.s[i]>b.s[i])
          return 3;
        if(a.s[i]<b.s[i])
          return 1;
    }
    return 2;
}

__attribute__((optimize("-O3"))) inline void Gcd()
{
    int flag;
    while(1)
    {
        flag=com();
        if(flag==2)
          break;
        while((a.s[1]&1)==0&&(b.s[1]&1)==0)
        {
            a=a/2;
            b=b/2;
            ++ji;
        }
        if(flag==1)
        {
            comp=((b.len-a.len-1)>0?(b.len-a.len-1):0);
            b=b-a;
        }
        else
        {
            comp=((a.len-b.len-1)>0?(a.len-b.len-1):0);
            a=a-b;
        }
    }
}

__attribute__((optimize("-O3"))) int main(){
    
    //freopen("1.txt","r",stdin);
    
    read();
    Gcd();
    for(int i=1;i<=ji;++i)
      a=a*2;
    out11(a);
}