POJ2115-C Looooops【扩欧,同余】

最新推荐文章于 2019-08-14 11:43:11 发布

QuantAsk

最新推荐文章于 2019-08-14 11:43:11 发布

阅读量231

点赞数

分类专栏：数论and数学文章标签：扩展欧几里得算法同余方程 POJ 数论

本文链接：https://blog.csdn.net/mr_wuyongcong/article/details/80741579

版权

数论and数学专栏收录该内容

261 篇文章 5 订阅

订阅专栏

正题

链接：
http://poj.org/problem?id=2115

大意

就是给出个循环

$f o r (i = A; i! = B; i = (i + C) m o d 2 k)$ $for (i=A;i!=B;i=(i+C)mod2^k)$

求需要循环次数

解题思路

我们定义 $l=2^k$
首先可以推出：

$C x + A \equiv B (m o d l)$ $Cx+A≡B(mod l)$

然后解mod

$C x + A = B + l k$ $Cx+A=B+lk$

然后定义 $y=-l$ ，移项

$C x + l y = B - A$ $Cx+ly=B-A$

然后我们定义 $d=gcd(C,l)$ ，之后同时除去d

$C x / d + l y / d = (A - B) / d$ $Cx/d+ly/d=(A-B)/d$

因为d是C和l的最大公约数那么因为 $C$ % $d=0$ , $l$ % $d=0$ ，所以只要 $(A-B)$ % $d=0$ 这个方程就有解
之后因为这样求出的x不是最大解所以我们要：
定义 $g=(l/d)$
然后

$(x*((B-A)/d)$ % $g+g)$ % g <script type="math/tex" id="MathJax-Element-28">g</script>

求出最小解

代码

#include<cstdio>
using namespace std;
long long x,y,d,a,b,c,k;
long long gcdup(long long a,long long b)
{
    if (b==0)
    {x=1;y=0;return a;}
    d=gcdup(b,a%b);
    long long k=x;
    x=y;
    y=k-a/b*y;
    return d;
}
int main()
{
    while (true)
    {
        scanf("%lld%lld%lld%lld",&a,&b,&c,&k);
        if (a==0 && b==0 && c==0 && k==0) break;
        k=1ll<<k;
        d=gcdup(c,k);
        if ((b-a)%d) printf("FOREVER\n");
        else printf("%lld\n",(x*((b-a)/d)%(k/d)+(k/d))%(k/d));
    }
}