线性同余方程组模板
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) 编辑 收藏 举报