luogu P2152 [SDOI2009]SuperGCD

题目传送门：https://www.luogu.org/problemnew/show/P2152

题意：

给两个正整数n,m（0 < n,m ≤ 10 ^ 10000），求n,m的最大公约数。

思路：

辗转相减法+优化。

即：对于求a,b的GCD(a,b)有：

若a为奇数，b为偶数，GCD(a,b)=GCD(a,b/2)

若a为偶数，b为奇数，GCD(a,b)=GCD(a/2,b)

若a为偶数，b为偶数，GCD(a,b)=2*GCD(a/2,b/2)

若a为奇数，b为奇数，GCD(a,b)=GCD(a-b,b) (a>b) 或 GCD(b-a,a) (b>a)

代码：

#include
#include
#define R register
#define I inline
    char s[10001];
    int a[10001],b[10001],d[10001];
    int la,lb,tot=0;
I void print(int a[10001],int p)
{
    for(R int i=p;i>=1;i--)
        printf("%c",a[i]+48);
}
I void chu(int a[10001],int &p)
{
    memset(d,0,sizeof(d));
    int r=0;
    for(R int i=p;i>=1;i--)
    {
        r=r*10+a[i];
        d[i]=(r>>1);
        r&=1;
    }
    while(!d[p]&&p>1)
        p--;
    memcpy(a,d,sizeof(d));
}
I void cheng()
{
    a[la+1]=0;
    for(R int i=la;i>=1;i--)
    {
        a[i]<<=1;
        a[i+1]+=a[i]/10;
        a[i]%=10;
    }
    while(a[la+1])
        la++;
}
I void jian(int a[10001],int b[10001],int &p)
{
    memset(d,0,sizeof(d));
    for(R int i=1;i<=p;i++)
    {
        if(a[i]1)
        p--;
    memcpy(a,d,sizeof(d));
}
I int pd()
{
    if(la>lb) return 1;
    if(la=1;i--)
    {
        if(a[i]>b[i]) return 1;
        if(a[i]