1. 程式人生 > >中國剩餘定理 (非互質情況)

中國剩餘定理 (非互質情況)

在前邊的時候介紹了普通剩餘定理的模板,但在很多時候是這些數並不互質,那次是有應該咋整呢,下邊我們就介紹一下不互質情況下的中國剩餘定理

首先不互質的情況就是兩個和成一個的過程


模板;

#include <iostream>
#include <cstdio>
using namespace std;
int a[10],b[10];
int gcd(int a,int b)
{
    return b?gcd(b,a%b):a;
}
void exgcd(int a,int b,int &x,int &y)
{
   if(b==0)
   {
       x=1;
       y=0;
       return ;
   }
   int x1,y1;
   exgcd(b,a%b,x1,y1);
   x=y1;
   y=x1-(a/b)*y1;
}
int main()
{
   int n,m,x,y;
    cin>>n;
    for(int h=1;h<=n;h++)
    {
        scanf("%d",&m);
        for(int i=1;i<=m;i++)
        scanf("%d",&a[i]);
        for(int j=1;j<=m;j++)
        scanf("%d",&b[j]);
        int m1=b[1],n1=a[1];
        int flag=1;
        for(int i=2;i<=m;i++)
        {
          int m2=b[i],n2=a[i];
          int d=gcd(n1,n2);
            if((m2-m1)%d!=0)
            {
                flag=0;
                break;
            }
            exgcd(n1,n2,x,y);
            int t=n2/d;
            x=x*((m2-m1)/d);
            x=(x%t+t)%t;
            m1=m1+x*n1;
            n1=(n1*n2)/d;
        }
        if(m1==0)
        {
            m1=1;
            for(int i=1;i<=m;i++)
                m1=m1*a[i]/gcd(m1,a[i]);
        }
        printf("Case %d: ",h);
        if(flag==0)
        printf("-1\n");
        else
        printf("%d\n",m1);
    }
    return 0;
}