https://vjudge.net/problem/HDU-3579
套模板
先可以先找两个同余方程 设通解为N;
N=r1(mod(m1)),N=r2(mod(m2));
显然可以化为k1*m1+r1=k2*m2+r2;—>k1*m1+(-k2*m2)=r2-r1;
设a=m1,b=m2,x=k1,y=(-k2),c=r2-r1方程可写为ax+by=c;
由欧几里得解得x即可,那么将x化为原方程的最小正整数解,(x*(c/d)%(b/d)+(b/d))%(b/d);
这里看不懂的去看解模线性方程。那么这个x就是原方程的最小整数解。
所以N=a*(x+n*(b/d))+r1====N=(a*b/d)*n+(a*x+r1),
这里只有n为未知数所以又是一个N=(a*x+r1)(mod(a*b/d))的式子,
然后只要不断的将两个式变成一个式子,最后就能解出这个方程组的解。
#include<iostream>
#include<cstdio>
#include<cstring>
u