中国剩余定理

最新推荐文章于 2024-04-27 09:09:22 发布

Owaski

最新推荐文章于 2024-04-27 09:09:22 发布

阅读量394

点赞数

分类专栏： Learning

本文链接：https://blog.csdn.net/qq_20118433/article/details/46406553

版权

Learning 专栏收录该内容

3 篇文章 0 订阅

订阅专栏

$~~~~~~$ 这玩意一直以为很难的 $...$ 结果看了几十分钟后发现就是个很简单的东西 $...$
$~~~~~~$ 先看一个题目：
$~~~~~~$ $\begin{matrix}x=a_1(mod~m_1)\\x=a_2(mod~m_2)\\...\\x=a_n(mod~m_n)\end{matrix}$
$~~~~~~$ 给出 $a_i,m_i$ ，求出所有的 $x$ 。
$~~~~~~$ 中国剩余定理证明了如果 $\forall m_i,m_j(i\neq j),\gcd(m_i,m_j)=1$ ，那么对于任意整数 $a_i$ ， $x$ 有解。
$~~~~~~$ 且解的形式如下：
$~~~~~~$ 设 $M=\prod m_i(1\leq i\leq n)$ ， $M_i=M/m_i$ ， $T_i(T_iM_i\equiv 1(mod~m_i))$
$~~~~~~$ 则有 $x=kM+ \sum a_iT_iM_i(1\leq i\leq n),k\in \mathbb{Z}$ 。

$~~~~~~$ 证明：

$~~~~~~$ 令 $x= \sum a_iT_iM_i(1\leq i\leq n)$ 。
$~~~~~~$ 则有 $a_iT_iM_i\equiv a_i(mod~m_i),a_jT_jM_j\equiv 0(mod~m_i)(i\neq j)$ 。
$~~~~~~$ 因此 $x$ 是满足上述方程组的。

$~~~~~~$ 令 $x_1, x_2$ 为方程组两个不同的解，那么：
$~~~~~~$ $\forall i,x_2-x_1\equiv 0(mod~m_i)\Rightarrow (所有的m_i互质)~x_2-x_1\equiv 0(mod~M)$ 。
$~~~~~~$ 所以任意两解之间都相差 $M$ 的整数倍，因此解的形式就是：
$~~~~~~$ $x=kM+ \sum a_iT_iM_i(1\leq i\leq n),k\in \mathbb{Z}$ 。

$~~~~~~$ 顺便扔道裸题上来，poj1006：

AC code：

#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <cstring>
#include <cctype>
#include <algorithm>
#include <string>
#include <sstream>
#include <iostream>
#include <map>
#include <set>
#include <list>
#include <stack>
#include <queue>
#include <vector>
#define pb push_back
#define mp make_pair
typedef long long LL;
typedef double DB;
typedef long double LD;
using namespace std;

const int MODP = 23;
const int MODE = 28;
const int MODI = 33;

int pow_mod(int a, int b, int p)
{
    int ret = 1, base = a;
    while(b)
    {
        if(b&1) ret = ret*base%p;
        base = base*base%p;
        b >>= 1;    
    }
    return ret;
}

int inv(int x, int p)
{
    if(p == 23) return pow_mod(x, p-2, p);
    else if(p == 28) return pow_mod(x, 11, p);
    else return pow_mod(x, 19, p);
}

int main()
{
    #ifndef ONLINE_JUDGE
    freopen("input.txt", "r", stdin);
    freopen("output.txt", "w", stdout);
    #endif

    int Test = 1;
    int p, e, i, d, x;
    int M, MP, ME, MI;
    int tp, te, ti;
    M = MODP*MODE*MODI;
    MP = M/MODP, ME = M/MODE, MI = M/MODI;
    tp = MP*inv(MP, MODP), te = ME*inv(ME, MODE), ti = MI*inv(MI, MODI);
    while(true)
    {
        scanf("%d%d%d%d", &p, &e, &i, &d);
        if(p == -1) break;
        x = p*tp+e*te+i*ti;
        if(x > d) x -= (x-d-1)/M*M;
        else x += (d-x+M)/M*M;
        x -= d;
        printf("Case %d: the next triple peak occurs in %d days.\n", Test++, x);
    }

    #ifndef ONLINE_JUDGE
    fclose(stdin);
    fclose(stdout);
    #endif
    return 0;
}