高精度gcd

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #define inf 1000000000
 using namespace std;
 char ch1[],ch2[];
 int la,lb,cnt;
 struct data{int a[],l;}a,b;
 bool com()
 {
     if(a.l<b.l)return ;
     if(a.l>b.l)return ;
     for(int i=a.l;i>;--i)
         if(a.a[i]>b.a[i])return ;
         else if(a.a[i]<b.a[i])return ;
     return ;
 }
 void print(data a)
 {
     while(a.a[a.l]==)a.l--;
     for(int i=a.l;i>;--i)
         if(i==a.l)printf("%d",a.a[i]);
         else printf("%09d",a.a[i]);
 }
 inline data sub(data a,data b)
 {
     int k;
     data c;
     for(int i=;i<=;++i)
     {
        if(i<=b.l)c.a[i]=a.a[i]-b.a[i];
        else if(i<=a.l)c.a[i]=a.a[i];
        else c.a[i]=;
        if(c.a[i]<)
        {
           c.a[i]+=inf;
           a.a[i+]--;
        }
     }
     c.l=a.l;
     while(c.a[c.l]==&&c.l)c.l--;
     return c;
 }
 void diva()
 {
     for(int i=;i<=a.l;i++)
     {
         if(a.a[i]&)a.a[i-]+=inf/;
         a.a[i]>>=;
     }
     if(!a.a[a.l])a.l--;
 }
 void divb()
 {
     for(int i=;i<=b.l;i++)
     {
         if(b.a[i]&)b.a[i-]+=inf/;
         b.a[i]>>=;
     }
     if(!b.a[b.l])b.l--;
 }
 void mul()
 {
     for(int i=a.l;i>;i--)
     {
        a.a[i]<<=;
        a.a[i+]+=a.a[i]/inf;
        a.a[i]%=inf;
     }
     while(a.a[a.l]>)a.l++;
     for(int i=b.l;i>;i--)
     {
        b.a[i]<<=;
        b.a[i+]+=b.a[i]/inf;
        b.a[i]%=inf;
     }
     while(b.a[b.l]>)b.l++;
 }
 int main()
 {
     scanf("%s%s",ch1+,ch2+);
     la=strlen(ch1+);lb=strlen(ch2+);
     if(la%)a.l=la/+;
     else a.l=la/;
     if(lb%)b.l=lb/+;
     else b.l=lb/;
     for(int i=;i<=a.l;++i)
     {
         int k1=max(,la-i*+),k2=la-(i-)*;
         for(int j=k1;j<=k2;++j)
             a.a[i]=a.a[i]*+ch1[j]-'';
     }
     for(int i=;i<=b.l;++i)
     {
         int k1=max(,lb-i*+),k2=lb-(i-)*;
         for(int j=k1;j<=k2;++j)
             b.a[i]=b.a[i]*+ch2[j]-'';
     }
     while()
     {
         if((a.a[]%==)&&(b.a[]%==)){diva();divb();cnt++;}
         else if((a.a[]%==))diva();
         else if((b.a[]%==))divb();
         if(com()){a=sub(a,b);if(!a.l){while(cnt--)mul();print(b);break;}}
         else {b=sub(b,a);if(!b.l){while(cnt--)mul();print(a);break;}}
     }
     return ;
 }