[SDOI2009][BZOJ 1876]SuperGCD

Description

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

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

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

Input

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

0 < A , B ≤ 10 ^ 10000。

Output

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

题解：

　　高精度

代码：

#include<bits/stdc++.h>
using namespace std;
inline int read(){
    int x=,f=;char ch=getchar();
    while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
    while(ch>=''&&ch<=''){x=(x<<)+(x<<)+ch-'';ch=getchar();}
    return x*f;
}
#define ll long long
#define MN 1255
#define Base 100000000
#define Width 8
char readin[MN*];
struct HPint{
    ll s[MN],len;
    void init(){memset(s,,sizeof s);len=;}
    HPint operator =(ll num){
        init();
        for(;num>;){len++;s[len]=num%Base;num/=Base;}
        return *this;
    }
    HPint operator =(const char *str){
        init();
        ll numlen=strlen(str);len=(numlen+Width-)/Width;
        for(int i=numlen-,t=,w;i>=;i--,w*=){
            if((numlen-i-)%Width==){w=;t++;}
            s[t]+=w*(str[i]-);
        }
        return *this;
    }
    bool operator <(const HPint &a)const{
        if(a.len!=len) return len<a.len;
        for(int i=a.len;i>=;i--)if(a.s[i]!=s[i])return s[i]<a.s[i];
        return ;
    }
    bool operator ==(const HPint &a)const{
        return !(a<*this)&&!(*this<a);
    }
    HPint operator -(const HPint &a)const{
        HPint ans;ans.init();ans.len=len;
        for(int i=;i<=len;i++){
            ans.s[i]+=s[i]-a.s[i];
            if(ans.s[i]<){ans.s[i]+=Base;ans.s[i+]--;}
        }
        for(;ans.s[ans.len]==&&ans.len>;ans.len--);
        return ans;
    }
    void MUL2(){
        for(int i=;i<=len;i++)s[i]*=;len+=;
        for(int i=;i<=len;i++){
            s[i+]+=s[i]/Base;
            s[i]%=Base;
        }
        for(;s[len]==&&len>;len--);
    }
    void DIV2(){
        len+=;
        for(int i=len;i>=;i--){
            if(s[i]&)s[i-]+=Base;
            s[i]/=;
        }
        for(;s[len]==&&len>;len--);
    }
    bool check(){
        if(len==)return ;
        else return s[]%==;
    }
    void read(){memset(readin,,sizeof(readin));scanf("%s",readin);*this=readin;}
    void print(){
        if(len==)printf("");
        else{printf("%lld",s[len]);for(int i=len-;i>=;i--)printf("%08lld",s[i]);}
    }
};
HPint a,b,res;
int cnt;
HPint gcd(HPint x,HPint y){
    while(!(x==y)){
        if(x<y)swap(x,y);
        if(x.check()&&y.check()) x.DIV2(),y.DIV2(),cnt++;
        else if(x.check()&&!y.check()) x.DIV2();
        else if(!x.check()&&y.check()) y.DIV2();
        else if(!x.check()&&!y.check()) x=x-y;
    }
    return x;
}
int main(){
    a.read();b.read();
    res=gcd(a,b);
    while(cnt--) res.MUL2();
    res.print();
    return ;
}

---

　　来自PaperCloud的博客，未经允许，请勿转载，TKS！