一本通1627:【例 3】最大公约数

#include <bits/stdc++.h>
using namespace std;
int a[2000], b[2000], maxnum = 1e8;
void jtll(int tmp[]) {
    string s;
    cin >> s;
    int len = s.length(), cnt = 0, ans = 0, t = 1, y = 1;

    for (int i = len - 1; i >= 0; i--) {
        cnt++;
        ans += t * (s[i] - 48);
        t *= 10;

        if (cnt == 8) {
            tmp[y++] = ans;
            ans = 0;
            cnt = 0;
            t = 1;
        }
    }

    if (cnt > 0)
        tmp[y] = ans;

    tmp[0] = (len + 7) / 8;
}
bool dlb() {
    if (a[0] > b[0])
        return true;

    if (a[0] < b[0])
        return false;

    for (int i = a[0]; i > 0; i--)
        if (a[i] > b[i])
            return true;
        else if (a[i] < b[i])
            return false;

    return true;
}
void SB(int a[], int b[]) {
    for (int i = 1; i <= a[0]; i++) {
        a[i] -= b[i];

        if (a[i] < 0)
            a[i] += maxnum, a[i + 1]--;
    }

    while (a[a[0]] == 0 and a[0] > 0)
        a[0]--;
}
void jntm(int a[]) {
    for (int i = a[0]; i >= 1; i--) {
        a[i - 1] += (a[i] % 2) * maxnum;
        a[i] /= 2;
    }

    while (a[a[0]] == 0 and a[0] > 0)
        a[0]--;
}
void muti(int a[], int m) {
    int g = 0;

    for (int i = 1; i <= a[0]; i++) {
        int x = a[i] * m + g;
        a[i] = x % maxnum;
        g = x / maxnum;

        if (i == a[0] and g > 0)
            a[0]++;
    }
}
void print(int a[]) {
    for (int i = a[0]; i > 0; i--)
        if (i == a[0])
            cout << a[i];
        else
            cout << setw(8) << setfill('0') << a[i];
}
int main() {
    jtll(a);
    jtll(b);
    int timer = 1;

    while (true) {
        if (a[1] % 2 == 0 and b[1] % 2 == 0) {
            jntm(a);
            jntm(b);
            timer *= 2;
            continue;
        } else if (a[1] % 2 == 0) {
            jntm(a);
            continue;
        } else if (b[1] % 2 == 0) {
            jntm(b);
            continue;
        }

        if (dlb()) {
            SB(a, b);

            if (a[0] == 0) {
                muti(b, timer);
                print(b);
                break;
            }
        } else {
            SB(b, a);

            if (b[0] == 0) {
                muti(a, timer);
                print(a);
                break;
            }
        }
    }

    return 0;
}

  • 9
    点赞
  • 8
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值