GCD
题目意思:
从 (x,y) 点出发,可以向上或向右走,每次走的距离为 LCM(x,y),即移动后的坐标为 (x+LCM(x,y),y) 或 (x,y+LCM(x,y)),
移动若干次后停止在 (ex,ey) 点,问有多少个点可以移动到点 (ex,ey)
本题要点:
1、加入当前的坐标是 (x1, y1), 那么下一个坐标的位置是 (x1 + lcm(x1, y1), y1) 或者 (x1, lcm(x1, y1) + y1),
记 k = gcd(x1, y1), 那么lcm(x1, y1) = kx1y1。 现在以第2个坐标为例子,
假设(x1, lcm(x1, y1) + y1) = (x, y), x和 y, 较大的值是有 lcm 累加到的。
gcd(x1, lcm(x1, y1) + y1) = k * gcd(x1, y1 + + x1 * y1), 注意到 gcd(x1, y1 + x1 * y1) = 1,
从而 x / k = x1, y / k = y1 + x1 * y1。
2、解法: 由当前的坐标 (x, y), (假如x < y)递推 回去上一个坐标(x1, y1)的过程
k = gcd(x, y), x1 = x / k, y1 = (y / k) / (x1 + 1);
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
int T, x, y, ans;
int gcd(int a, int b)
{
return b == 0 ? a : gcd(b, a % b);
}
int lcm(int a, int b)
{
return a / gcd(a, b) * b;
}
void solve()
{
ans = 1;
if(x > y)
swap(x, y);
int k = gcd(x, y), x1 = x / k, Y = y / k, y1;
while(Y % (x1 + 1) == 0)
{
++ans;
y1 = Y / (x1 + 1);
x = x1, y = y1;
if(x > y)
swap(x, y);
k = gcd(x, y), x1 = x / k, Y = y / k;
}
printf("%d\n", ans);
}
int main()
{
scanf("%d", &T);
for(int t = 1; t <= T; ++t)
{
scanf("%d%d", &x, &y);
printf("Case #%d: ", t);
solve();
}
return 0;
}
/*
3
6 10
6 8
2 8
*/
/*
Case #1: 1
Case #2: 2
Case #3: 3
*/