POJ 2115 C Loooooops

最新推荐文章于 2021-07-19 01:00:41 发布

skajre

最新推荐文章于 2021-07-19 01:00:41 发布

阅读量325

点赞数

分类专栏：数论：GCD

本文链接：https://blog.csdn.net/qq_33901573/article/details/52062079

版权

数论：GCD 专栏收录该内容

3 篇文章 0 订阅

订阅专栏

//
//  main.cpp
//  Richard
//
//  Created by 邵金杰 on 16/7/29.
//  Copyright © 2016年 邵金杰. All rights reserved.
//



#include<cstdio>
#include<cmath>
typedef long long LL;
LL a,b,c,k;
void gcd(LL a,LL b,LL &d,LL &x,LL &y)
{
    if(!b) {d=a;x=1;y=0;}
    else {gcd(b,a%b,d,y,x);y-=(a/b)*x;}
}
void cal(LL a,LL b,LL c)
{
    LL d,x,y;
    gcd(a,b,d,x,y);
    if(c%d!=0){ printf("FOREVER\n");return ;}
    LL s=b/d;
    LL k=c/d;
    LL ans=k*x;
    printf("%lld\n",(ans%s+s)%s);
}
int main()
{
    while(scanf("%lld%lld%lld%lld",&a,&b,&c,&k))
    {
        if(a==0&&b==0&&c==0&&k==0) break;
        if(b==0) {printf("0\n");continue;}
        cal(c,(1LL<<(k)),b-a);//这里不能用(1<<k)否则会溢出，应使用(1LL<<k),那么1代表的是long long int
    }
    return 0;
}