HDU 5584 LCM Walk

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5584

题意：有一只青蛙在(x,y)点，每次它只能跳向（x+z,y）或者(x,y+z)  z = lcm(x,y)。现在知道它跳了若干步后在(ex,ey)点，求可能的起点个数。

思路：设x = p * k , y = q * k ( k = gcd(x,y)  , p  q  互质 )  那么z = p*q*k. 下一步会跳到( k * p * (1+q)  , q * k  ) 或者 ( p * k  ,  k * q * ( 1 + p )  )，我们注意到下一步的点x,y的gcd还是k，因为p,q都是质数，所以我们可以根据当前步推出上一步的x,y。先把gcd都除掉，然后用较大的值/(较小的值+1)，另一个值则不变，这样一直循环下去即可。

#include <cstdio>
#include <cmath>
#include <cstring>
#include <string>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <stack>
#include <map>
#include <set>
#include <vector>
#include <sstream>
#include <queue>
#include <utility>
using namespace std;

#define rep(i,j,k) for (int i=j;i<=k;i++)
#define Rrep(i,j,k) for (int i=j;i>=k;i--)

#define Clean(x,y) memset(x,y,sizeof(x))
#define LL long long
#define ULL unsigned long long
#define inf 0x7fffffff
#define mod 100000007

int gcd(int a,int b)
{
    return b==0?a:gcd(b,a%b);
}
int a,b;

int solve()
{
    int ans = 0;
    while(1)
    {
        int k = gcd( a , b );
        if ( a > b ) swap(a,b);
        a/=k;
        b/=k;
        if ( b % (a+1) == 0 )
        {
            ans++;
            b/=(a+1);
        }
        else break;
    }

    return ans+1;
}

int main()
{
    int T;
    cin>>T;
    rep(kase,1,T)
    {
        scanf("%d%d",&a,&b);
        printf("Case #%d: %d\n",kase,solve());
    }
    return 0;
}