[SDOI2009]SuperGCD

题目链接

这题。高精度。恶心。难受。

那么高精度的gcd怎么做呢？

若a=b gcd(a,b)=a

①a偶b偶 gcd(a,b)=2*gcd(a/2,b/2)

②a偶b奇 gcd(a,b)=gcd(a/2,b)

③a奇b奇 gcd(a,b)=gcd(a-b,b)

嗯。这玩意就这样了。

#include<cmath>
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#define ll long long
using namespace std;
int a[],b[],c[],f[],s0,T,kpl[],kyl;
char s[];
void bigscanf(int *a){
    scanf("%s",s);
    int len=strlen(s),i,j;
    for(i=;i<len;++i){
        j=(len-i+)/;
        a[j]=a[j]*+s[i]-'';
    }
    a[]=(len+)/;
}
void bigprintf(int *a){
 cout<<a[a[]];
 for(int i=a[]-;i>;--i){
     for(int j=;j>;j/=)cout<<a[i]/j%;
    }
    cout<<endl;
}
int bigcmp(int *a,int *b){
    if(a[]>b[])return ;
    if(a[]<b[])return -;
    for(int i=a[];i>;i--){
        if(a[i]>b[i])return ;
        if(a[i]<b[i])return -;
    }
    return ;
}
void bigsub(int *a,int *b){
    int k=a[],g=;
    for(int i=;i<=k;i++){
        a[i]=a[i]-b[i];
        if(a[i]<){
            a[i+]--;
            a[i]=a[i]+;
        }
        else g=;
    }
    while(k>&&a[k]==)k--;
    a[]=k;
}
void bigmul2(int *a,int *b,int *c){
    int k=a[]+b[];
    for(int i=;i<=a[];i++){
        for(int j=;j<=b[];j++){
            c[i+j-]=c[i+j-]+a[i]*b[j];
            c[i+j]=c[i+j]+c[i+j-]/;
            c[i+j-]=c[i+j-]%;
        }
    }
    while(k>&&c[k]==)k--;
    c[]=k;
}
void bigmul1(int *a,int b){
    int k=a[],g=;
    for(int i=;i<=a[];i++){
        a[i]=a[i]*b+g;
        g=a[i]/;
        a[i]=a[i]%;
    }
    while(g>){
        k++;
        a[k]=g%;
        g=g/;
    }
    while(k>&&a[k]==)k--;
    a[]=k;
}
void bigdiv1(int *a,int b){
    int k=a[],d=;
    for(int i=a[];i>=;i--){
        d=d*+a[i];
        a[i]=d/b;
        d=d%b;
    }
    while(k>&&a[k]==)k--;
    a[]=k;
}
void copy(int *a,int *b){
    for(int i=;i<=b[];++i){
        a[i]=b[i];
    }
}
void Gcd(int *a,int *b,int t){
    int u=bigcmp(a,b);//cout<<u<<endl;
    if(u==){T=t;return;}
    if(u<){Gcd(b,a,t);return;}
    int ta=,tb=;
    if(a[]%==){
        bigdiv1(a,);
        ta=;
    }
    if(b[]%==){
        bigdiv1(b,);
        tb=;
    }
    if(ta&&tb)Gcd(a,b,t+);
    else if(! ta&&! tb){bigsub(a,b);Gcd(a,b,t);}
    else Gcd(a,b,t);
}
int main(){
    bigscanf(a);
    bigscanf(b);
    Gcd(a,b,);
    if(T==)bigprintf(a);
    else {
        f[]=f[]=;
        for(int i=;i<=T;i++)bigmul1(f,);
        bigmul2(f,a,kpl);
        copy(f,kpl);
        bigprintf(f);
    }
return ;
}