[模板]求解线性同余方程

这一生的挚爱 提交于 2019-11-29 00:18:11

求解a*x = b (mod p)

#include<bits/stdc++.h>
using namespace std;
int exgcd(int a,int b,int& x,int& y)
{
    int d=a;
    if(b){
        d=exgcd(b,a%b,y,x);
        y-=(a/b)*x;
    }
    else{
        x=1;y=0;
    }
    return d;
}
int gcd(int a,int b)
{
    return !b?a:gcd(b,a%b);
}
int main()
{
    int a,b,p;
    scanf("%d%d%d",&a,&b,&p);//计算 a*x = b (mod p)
    int x,y;
    int d=exgcd(a,p,x,y);
    if(b%d!=0)    printf("-1\n");
    else
    {
        x=(x)*b/d;
        x+=p;
        x%=p;
        for(int i=1;i<10;i++)
        {
            cout<<(x+(i-1)*p/d)<<endl;
        }
    }
}

 

易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!