2015 ACM/ICPC 上海区域赛 L—LCM Walk【数论LCM】

https://icpcarchive.ecs.baylor.edu/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=5433

题意:

开始时有一个格子(x,y),每次移动会移动到(x,y+lcm(x,y))或(x+lcm(x,y),y),现在给你最终的点,让你求可能的起点个数。

x,y<=1e9

假设前一个格子是(x,y) ,令 x=a*k,y=b*k,(k=gcd(x,y),a,b,互质)

则lcm(x,y)=abk。不妨设我们下一步往右走,

则下一个点坐标为 (x,y+lcm(x,y)) = (a*k,b*k+a*b*k) = k*(a,b*(a+1))  ,

如果(aa,bb)存在bb>aa,坐标可以表示为(a*k,(b+a*b)*k),那么肯定有上一步(a,b)

k=gcd(aa,bb)。如果yy/k%(a+1)==0,一定有上一点 a=xx/k,bb=yy/k/(a+1),

注意: x==y时特判为1。

代码:

#include<bits/stdc++.h>
#define ll long long
using namespace std;
ll ans,x,y;
ll gcd(ll x,ll y){return y==0?x:gcd(y,x%y);}
ll dfs(ll x,ll y)
{
    if(x==y) return 1;
    if(x>y) swap(x,y);
    ll g=gcd(x,y);
    x/=g;y/=g;
    if(y%(x+1)==0) return 1+dfs(g*x,g*y/(x+1));
    else return 1;
}
int main()
{
    int T,cas=1;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%lld%lld",&x,&y);
        printf("Case #%d: %lld\n",cas++,dfs(x,y));
    }
}

 

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值