poj2115 同余方程

最新推荐文章于 2022-02-03 23:23:34 发布

zxc106

最新推荐文章于 2022-02-03 23:23:34 发布

阅读量341

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/u013983192/article/details/37771679

版权

数论专栏收录该内容

11 篇文章 0 订阅

订阅专栏

这道题要是有扩展的欧几里得算法基础的，不难想出来思路。

没有的就不要瞎想了，先看看扩展的欧几里得算法，在看看同余方程再来AC吧~

ax=b (mod n)

该方程有解的充要条件为 gcd(a,n) | b ，即 b% gcd(a,n)==0

令d=gcd(a,n)

有该方程的最小整数解为 x = e (mod n/d)

其中e = [x0 mod(n/d) + n/d] mod (n/d) ，x0为方程的最小解

那么原题就是要计算b% gcd(a,n)是否为0，若为0则计算最小整数解，否则输出FOREVER

好，我们目标出现了，就是求x0，方程的最小解。这就是扩展的欧几里得算法要做的事情：

对于：a*x+b*y=gcd(a,b)

1.当b==0时，gcd(a,b)=a，此时x=1,y=0;

2.先求出 a*x+b*y=gcd(a,b) 的一组解。

因为 a*x1+b*y1=gcd(a,b)

b*x2+a%by2=gcd(b,a%b)

且 gcd(a,b)=gcd(b,a%b);

所以有a*x1+b*y1=b*x2+(a-(a/b)*b)*y2

从而得x1=y2，y1=x2-(a/b)*y2

然后执行程序段：
void expgcd(int a,int b,int &x,int &y)
{
if(b==0)

     {
          x=1;
          y=0;
          return ;
     }
     expgcd(b,a%b,x,y);
     int t=x;
     x=y;
     y=t-(a/b)*y;
}
得出一组解x0，y0;

剩下的事情就是把这些求出来的参数往一起凑凑，就有结果了。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
long long X,Y;
long long gcd(long long a,long long b)  //求最大公约数
{
    if(a==0) return b;
    return gcd(b%a,a);
}
void expgcd(long long a,long long b,long long &x,long long &y)  //求d=ax+by 的x和y。
{
    if(b==0)
    {
        x=1;
        y=0;
        return ;
    }
    expgcd(b,a%b,x,y);
    long long t=x;
    x=y;
    y=t-a/b*y;
}
int main()
{
    long long a,b,c,d,k,sum,B;
    long long m;
    while(cin>>a>>b>>c>>k)
    {
        sum=a+b+c+k;
        if(sum==0) break;
        m=1LL<<k;
        B=(b-a+m)%m;
        d=gcd(c,m);
        if(B%d!=0) {cout<<"FOREVER"<<endl;continue;}
        expgcd(c,m,X,Y);
        X=X*(B/d)%m;
        X=(X%(m/d)+m/d)%(m/d); //方程cx=B(mod m)的最小整数解
        cout<<X<<endl;
    }
    return 0;
}