中国剩余定理
互质

长春网站建设公司成都创新互联,长春网站设计制作,有大型网站制作公司丰富经验。已为长春成百上千家提供企业网站建设服务。企业网站搭建\外贸网站建设要多少钱,请找那个售后服务好的长春做网站的公司定做!
ll m[N],p[N];//p是质数,m是余数
ll Exgcd(ll a,ll b,ll &x,ll &y)
{
if(b==0)
{
x=1;
y=0;
return a;
}
ll x1,y1;
ll d=Exgcd(b,a%b,x1,y1);
x=y1;
y=x1-a/b*y1;
return d;
}
ll Crt(ll *m,ll *p,int l)
{
ll res=0,n=1,x,y;
for(int i=0;i非互质
#include//*m模数 *a余数
#include
#include
#include
#include
using namespace std;
typedef long long ll;
void exgcd(ll a,ll b,ll &x,ll &y)
{
if(!b)
{
x=1;y=0;return ;
}
exgcd(b,a%b,x,y);
ll temp=x;
x=y;y=temp-(a/b)*y;
}
ll inv(ll a,ll b)
{
ll d=__gcd(a,b);
if(d!=1)return -1;
ll x,y;
exgcd(a,b,x,y);
return (x%b+b)%b;
}
bool merge(ll a1,ll m1,ll a2,ll m2,ll &a3,ll &m3)
{
ll d=__gcd(m1,m2);
ll c=a2-a1;
if(c%d)return false;
c=(c%m2+m2)%m2;
m1/=d;m2/=d;c/=d;
c*=inv(m1,m2);c%=m2;
c*=m1*d;
c+=a1;m3=m1*m2*d;
a3=(c%m3+m3)%m3;
return true;
}
ll CRT(ll *a,ll *m,ll n)
{
ll a1=a[1],m1=m[1];
for(ll i=2;i<=n;i++)
{
ll a2=a[i],m2=m[i];
ll m3,a3;
if(!merge(a1,m1,a2,m2,a3,m3))
return -1;
a1=a3;m1=m3;
}
return (a1%m1+m1)%m1;
}
文章名称:中国剩余定理
网站链接:http://www.jxjierui.cn/article/ispiso.html


咨询
建站咨询
