HDU-3579 Hello Kiki（中国剩余定理）

题意

给定一个实数 $n$ ，一些硬币 $MiMi$ 个数还剩下 $Ai$ 个 $(i\in[1,n])$，求硬币的最小数目。

思路

中国剩余定理模板题，从头梳理一遍。

exgcd函数

求的是方程 $ax+by=gcd(a,b)$ 的一组特解 $(x,y)$。
首先根据 $gcd$ 函数的原理 $gcd(a,b)=gcd(b,a\%b)$
那么列出两组方程并联立
$ax_1+by_1=gcd(a,b)$ 且 $ax_2+by_2=gcd(b,a\%b)$
故 $ax_1+by_1=bx_2+a\%by_2$
$ax_1+by_1=bx_2+(a-[a/b]*b)y_2$
$a(x_1-y_2)+b(y_1-x2+[a/b]*y_2)=0$
求出上式的一组特解，不难看出 $x_1=y_2, y_1=x_2-[a/b]*y_2$ 时方程成立。
把前式代入后式，得到 $x_1=y_2, y_1=x_2-[a/b]*x_1$
这是最容易代码实现的方程，由于 $x_1,y_1$ 的结果来源于 $x_2,y_2$ 所以递归时把 $x,y$ 引用往下传。

Exgcd函数

求的是方程 $Ax+By=C$ 中 $x$ 的最小非负整数解。
在这里进行解的存在性判断，令 $g=gcd(A,B)$ 。首先，若 $C$ 不能整除 $g$ ，那这个方程很明显是无解的。
将 $A,B,C$ 同除以 $g$。
由上面的 $exgcd$ 函数，不难求出 $Ax+By=gcd(A,B)=1$ 的特解（ $A,B$ 已经互质）$x',y'$。
那么通解就是 $x=x'+Bk$ ， $y=y'-Ak$ 。
方程 $Ax+By=C$ 的通解，就是 $x=Cx'+Bk$ ， $y=Cy'-Ak$ 了。
那么 $x$ 的最小非负整数解，即为 $x=(Cx'\%B+B)\%B$ 。

kase_insert函数

这个函数用来合并两个形如 $X=Pk+A$ 的方程（即 $X$ 模 $P$ 得 $A$）。
令这两个方程为
$X=P_1k_1+A_1$
$X=P_2k_2+A_2$
合并后的方程为
$X=P_3k_3+A_3$
联立化简，最后结果为
$P_1k_1-P_2k_2=A2-A1$
使用 $Exgcd$ 函数，求得 $k_1$ 的最小非负整数解 $k_1'$ 。
$A_3$ 即为 $P_1k_1+A_1$ ， $P_3$ 比较明显，是 $lcm(P_1,P_2)$。

自此，整个模板的函数已大致介绍完毕。

代码

#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define FOR(i,x,y) for(int i=(x);i<=(y);i++)
#define DOR(i,x,y) for(int i=(x);i>=(y);i--)
typedef long long LL;
using namespace std;
struct CRT
{
    LL A,P;
    LL gcd(LL a,LL b){return b?gcd(b,a%b):a;}
    void exgcd(LL a,LL b,LL &x,LL &y)
    {
        if(!b){x=1,y=0;return;}
        exgcd(b,a%b,y,x);y-=a/b*x;
        return;
    }
    LL Exgcd(LL A,LL B,LL C)
    {
        LL k1,k2,g=gcd(A,B);
        if(C%g)return -1;
        A/=g,B/=g,C/=g;
        exgcd(A,B,k1,k2);
        return (k1*C%B+B)%B;
    }
    bool kase_insert(CRT _)
    {
        LL res=Exgcd(P,_.P,_.A-A);
        if(!~res)return 0;
        A+=P*res;
        P=P/gcd(P,_.P)*_.P;
        return 1;
    }
    LL ask(){return A?A:P;}
};

int main()
{
    int T;
    scanf("%d",&T);
    FOR(Ti,1,T)
    {
        int n;
        LL M[10],A[10];
        scanf("%d",&n);
        FOR(i,1,n)scanf("%lld",&M[i]);
        FOR(i,1,n)scanf("%lld",&A[i]);
        CRT sum=(CRT){A[1],M[1]};
        bool flag=0;
        FOR(i,2,n)if(!sum.kase_insert((CRT){A[i],M[i]}))
        {
            flag=1;
            break;
        }
        printf("Case %d: %lld\n",Ti,flag?-1:sum.ask());
    }
    return 0;
}