解题报告（十三）中国剩余定理（ACM / OI）

整理的算法模板合集： ACM模板

点我看算法全家桶系列！！！

实际上是一个全新的精炼模板整合计划

---

繁凡出品的全新系列：解题报告系列 —— 超高质量算法题单，配套我写的超高质量的题解和代码，题目难度不一定按照题号排序，我会在每道题后面加上题目难度指数（$1\sim 5$），以模板题难度 $1$ 为基准。

---

这样大家在学习算法的时候就可以执行这样的流程：

%
阅读【学习笔记】 / 【算法全家桶】学习算法 $\Rightarrow$ 阅读相应算法的【解题报告】获得高质量题单 $\Rightarrow$ 根据一句话题解的提示尝试自己解决问题 $\Rightarrow$ 点开详细题解链接学习巩固（好耶）
%
要是26个英文字母用完了我就接上24个希腊字母，我就不信50道题不够我刷的hhh
%
解题报告系列合集：【解题报告系列】超高质量题单 + 题解（ICPC / CCPC / NOIP / NOI / CF / AT / NC / P / BZOJ）

本题单前置知识：《算法竞赛中的初等数论》（二）正文 0x20同余（ACM / OI / MO）（十五万字符数论书）

A、（P3868 [TJOI2009]）猜数字

Weblink

https://www.luogu.com.cn/problem/P3868

Problem

Solution

$$\begin{array}{rl}& b_i|(n-a+i)\\ & n-a_i\%b_i=0\\ & n-a_i\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}{b}_i)\\ & n\equiv a_i(\mathrm{m}\mathrm{o}\mathrm{d}{b}_i)\end{array}$$

显然有 $n$ 个同余方程，直接CRT解方程组即可。

注意数据保证 $\prod _{i=0}^k{b}_i\le 10^{18}$，即 $M\le 10^{18}$，那么CRT的过程中随便一乘就会爆 long long ，所以要用快速乘。用 $logn$ 的快速乘竟然 T 了…所以需要加一些经典优化或者用 $O(1)$ 的快速乘，可以处理小于 $1.7\times 10^{308}$ 的数据

Code

#include <bits/stdc++.h>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <map>
#include <queue>
using namespace std;
#define int long long
typedef long long ll;
typedef int itn;
typedef pair<int, int> PII;
typedef pair<double, int> PDI;
const int N = 50 + 7, mod = 1e18 + 7, INF = 0x3f3f3f3f;
const double PI = acos(-1.0), eps = 1e-8;

int n, t, kcase;
int M, m[N], a[N], k;
/*
int mul(int a, int b, int mod)
{
    int res = 0;
    while(b) {
        if(b & 1) res = (res + a) % mod;
        a = (a + a) % mod;
        b >>= 1;
    }
    return res;
}*/

int mul(int x, int y, int p)
{
   
    int z = (long double)x / p * y;
    ll res = (unsigned long long)x * y - (unsigned long long) z * p;
    return (res + p) % p;
}

ll exgcd(int a, int b, int &x, int &y)
{
   
    if(b == 0) {
   
        x = 1, y = 0;
        return a;
    }
    ll d = exgcd(b, a % b, x, y);
    ll z = x;
    x = y;
    y = z - y * (a / b);
    return d;
}

void solve()
{
   
    ll M = 1;
    scanf("%lld"