51Nod－1247－可能的路径

最新推荐文章于 2022-10-06 14:15:45 发布

f_zyj

最新推荐文章于 2022-10-06 14:15:45 发布

阅读量1k

点赞数

分类专栏：数论数学相关 51Nod-题解集锦文章标签：数论

本文链接：https://blog.csdn.net/f_zyj/article/details/52157897

版权

51Nod-题解集锦同时被 3 个专栏收录

460 篇文章 35 订阅

订阅专栏

数论

235 篇文章 1 订阅

订阅专栏

数学相关

182 篇文章 3 订阅

订阅专栏

ACM模版

描述

题解

神级模拟推导题……经过模拟推导可以推出：

(a, b) -> (b, a) ··············································(1)
if (a, b) -> (x, y) to (x, y) -> (a, b) ··········(2)

那么我们可以逐次推导出(假设a>b)：

(a, b) -> (a - b, b) -> (a - 2b, b) -> … -> (a - nb, b) 其中，n = a / b
(a, b) - > (a % b, b) -> (b, a % b)
由此可以联想到欧几里得算法求解最大公约数！
如果(a, b)的最大公约数是c，那么(a, b)一定可以达到(c, c)
当(x, y)的最大公约数也是c时，说明(x, y)也可以到达(c, c)
如此，说明两者存在路径连通，输出”Yes”；反之输出”No”

代码

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cmath>

typedef long long ll;
using namespace std;

ll kgcd(ll a, ll b)
{
    if (!a)
    {
        return b;
    }
    else if (!b)
    {
        return a;
    }
    else if (!(a & 1) && !(b & 1))
    {
        return kgcd(a >> 1, b >> 1) << 1;
    }
    else if (!(b & 1))
    {
        return kgcd(a, b >> 1);
    }
    else if (!(a & 1))
    {
        return kgcd(a >> 1, b);
    }
    else
    {
        return kgcd(abs(a - b), min(a, b));
    }
}

int main(int argc, const char * argv[])
{
    int T;
    cin >> T;

    ll a, b, x, y;
    while (T--)
    {
//        scanf("%lld %lld %lld %lld", &a, &b, &x, &y);
        cin >> a >> b >> x >> y;
        if (kgcd(a, b) == kgcd(x, y))
        {
            printf("Yes\n");
        }
        else
        {
            printf("No\n");
        }
    }

    return 0;
}