这里写代码片
///用于计算ax==b(mod m)的一个特解,如果没有特解返回m本身
LL gcd_ex(LL a,LL b,LL &x,LL &y)
{
if(b==0) {x=1;y=0;return a;}
LL d=gcd_ex(b,a%b,y,x);
y=y-a/b*x;
return d;
}
LL solve(LL a,LL b,LL m)
{
LL x,y;LL d=gcd_ex(a,m,x,y);
if(b%d==0){x%=m;while(x<0) x+=m;x%=m;return x*(b/d)%(m/d);}
else return m;
}
///用于计算ax==b(mod m)在[0,m)的所有解,所有解存在ans容器中,返回ans;
///输出时的定义 vector<LL>ans=solve(a,b,m);输出循环for(int i=0;i<ans.size();i++) 输出ans[i];
LL gcd_ex(LL a,LL b,LL &x,LL &y)
{
if(b==0) {x=1;y=0;return a;}
LL d=gcd_ex(b,a%b,y,x);
y=y-a/b*x;
return d;
}
vector<LL> solve(LL a,LL b,LL m)
{
LL x,y;LL gcd=gcd_ex(a,m,x,y);
vector<LL>ans;ans.clear();
if(b%gcd==0)
{
x%=m;while(x<0) x+=m;x%=m;
ans.push_back(x*(b/gcd)%(m/gcd));
for(LL i=1;i<gcd;i++)
ans.push_back((ans[0]+i*m/gcd)%m);
}
return ans;
}
///中国剩余定理
/// 用于计算x==a[i](mod m[i]) 的一个特解,要求所有的模两两互质。
///简单的求解过程如下:
///1.M0=m[1]*m[2]*...m[n];
///c[i]是方程m[i]*x==1(mod m[i])的一个特解,其中M[i]=M0/m[i];
///x==a[1]*c[1]*M[1]+a[2]*c[2]*M[2]+....a[n]*c[n]*M[n](mod M0);
///时间复杂度nlogM;
///返回特解的值;
void gcd_ex(LL a,LL b,LL &x,LL &y)
{
if(b==0) {x=1;y=0;}
else { gcd_ex(b,a%b,y,x); y-=a/b*x;}
}
LL CRT(LL *a,LL *m,LL n)
{
LL M=1,res=0;
for(int i=0;i<n;i++) M*=m[i];
for(int i=0;i<n;i++)
{
LL x,y,temp=M/m[i];
gcd_ex(temp,m[i],x,y);
res=(res+temp*x*a[i])%M;
}
return (res+M)%M;
}
求解ax===b(mod m)的特解和[0,m)的所有解
最新推荐文章于 2022-10-30 19:51:43 发布