hdu 3579 Hello Kiki 中国剩余定理(不互质的情况) 对互质的情况,处理起来比较方便,可以直接套模板 本题给出不互质的模线性方程组,求出满足方程的最小正整数解 方案:对于不互质的模线性方程组,可以进行方程组合并,求出合并后的方程的解,这样就可以很快地推出方程的最终解。 两个方程合并的一种方法: x = c1 (mod b1) x = c2(mod b2) 此时b1,b2不必互质的。 显然可以得到x = k1 b1 + c1 x = k2 b2 + c2, 两个方程合并一下就可以得到:k1 b1 = c2 - c1 (mod b2), 这样可以设g=gcd(b1,b2),于是就有b1/gk1-b2/gk2=(c2-c1)/g, 显然判断(c2-c1)/g是否为整数就能判断是否存在解, 这样在经过类似的变换就能得到k1 = K (mod (b2/g)), 最后得到x = Kb1 + c1 (mod (b1 * b2/g))。 对于题目所给正整数的要求,只有一种反例,就是结果输出为0的情况, 这个可以特殊考虑,只需要考虑所有数的最小公倍数即可。
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
using namespace std;
__int64 x, y, t;
__int64 egcd(__int64 a, __int64 b)
{
if (b==0)
{
x=1;
y=0;
return a;
}
else
{
__int64 e=egcd(b,a % b);
t=x;
x=y;
y=t-a/b*y;
return e;
}
}
__int64 gcd(__int64 x, __int64 y)
{
if (!x || !y)
return x > y ? x : y;
for (__int64 t; t = x % y; x = y, y = t);
return y;
}
__int64 mm[10],rr[10];
int main()
{
int T,Case,N;
__int64 m1,m2,r1,r2,d,c,t;
bool flag;
scanf ("%d",&T);
for(Case = 1; Case <= T; Case++)
{
scanf ("%d",&N);
flag=0;
for (__int64 i=0;i<N;i++)
{
scanf ("%I64d",&mm[i]);
}
for (__int64 i=0;i<N;i++)
{
scanf ("%I64d",&rr[i]);
}
m1=mm[0];
r1=rr[0];
for (__int64 i = 0; i < N - 1; i++)
{
m2=mm[i+1];
r2=rr[i+1];
if (flag)
continue;
d = egcd(m1, m2);
c = r2 - r1;
if (c % d)
{
flag = 1;
continue;
}
t=m2/d;
x=(c/d*x%t+t)%t;
r1=m1*x+r1;
m1=m1*m2/d;
}
if (flag)
printf ("Case %d: -1\n",Case);
else
{
if (r1==0&&N>1)
{
r1=mm[0];
__int64 ans=1;
for (int i=1;i<N;i++)
r1=gcd(mm[i],r1);
for (int i=0;i<N;i++)
ans*=mm[i];
r1=ans/r1;
}
if (r1==0&&N==1)
r1=mm[0];
printf ("Case %d: %I64d\n",Case,r1);
}
}
return 0;
}