线性同余方程组模板

void exgcd(int a,int b,int &d, int &x,int &y){///扩展欧几里德算法
    if(!b){x=1;y=0;d=a;}
    else{exgcd(b,a%b,d,y,x);  y-=x*(a/b);}
}
int n,m;///m代表方程组的数量
int a[11],b[11];///x%a[i]=b[i]
///求解同余方程组
bool Chi_r(int &minnum,int &mod){///minnum表示最小正整数。mod表示模
    int m1=a[0],b1=b[0],x,y,d;
    for(int i=1;i<m;++i){
        int m2=a[i],b2=b[i];
        exgcd(m1,m2,d,x,y);
        if((b2-b1)%d!=0) return false;///不存在解
        int t=m2/d;
        x=(x*((b2-b1)/d)%t+t)%t;
        b1=m1*x+b1;
        m1=m1*(m2/d);
    }
    minnum=b1;
    mod=m1;
    return true;///存在解
}

posted on 2018-03-21 02:02  yjbjingcha  阅读(120)  评论(0编辑  收藏  举报

导航