poj2115扩展欧几里得

题目： 点击打开链接

原文： 点击打开链接

大致题意：

对于C的for(i=A ; i!=B ;i +=C)循环语句，问在k位存储系统中循环几次才会结束。

若在有限次内结束，则输出循环次数。

否则输出死循环。

 

解题思路：

题意不难理解，只是利用了 k位存储系统 的数据特性进行循环。

例如int型是16位的，那么int能保存2^16个数据，即最大数为65535（本题默认为无符号），

当循环使得i超过65535时，则i会返回0重新开始计数

如i=65534，当i+=3时，i=1

其实就是 i=(65534+3)%(2^16)=1

 

有了这些思想，设对于某组数据要循环x次结束，那么本题就很容易得到方程：

x=[(B-A+2^k)%2^k] /C

即 Cx=(B-A)(mod 2^k)  此方程为 模线性方程，本题就是求X的值。

#include<iostream>
#include<cstdio>
using namespace std;
#define LL long long
void ex_gcd(LL a, LL b, LL &d, LL &x, LL &y)
{
    if(!b)
    {
        d = a;
        x = 1;
        y = 0;
    }
    else
    {
        ex_gcd(b, a%b, d, y, x);
        y -= x*(a/b);
    }
}
int main()
{
    LL A,B,C,k;
    while(~scanf("%lld%lld%lld%lld",&A,&B,&C,&k))
    {
        if(A==0&&B==0&&C==0&&k==0)
            break;
        LL a,b,d,x,y;
        a = C;
        b = 1LL<<k;
        ex_gcd(a,b,d,x,y);
        if((B-A)%d != 0)
        {
            printf("FOREVER\n");
            continue;
        }
        x = x * (B-A) / d;
        x = (x%(b/d)+b/d)%(b/d);//求最小正整数解
        printf("%lld\n",x);
    }
    return 0;
}