题目大意:给你一些关于$x$的方程组:
$$
\begin{cases}
x\equiv a_1\pmod{mod_1}\\
x\equiv a_2\pmod{mod_2}\\
\vdots\\
x\equiv a_n\pmod{mod_n}
\end{cases}
$$
求解$x$的最小非负整数解($\gcd(mod_1,mod_2,\cdots,mod_n)\not=1$)

题解:$EXCRT$,假设有两个方程
$$
\begin{cases}
x\equiv x_1\pmod{A}\\
x\equiv x_2\pmod{B}
\end{cases}\\
设g=\gcd(A,B),k=\left\lfloor\dfrac xg\right\rfloor,c=x\bmod g\\
x=kg+c\\
所以x_1\equiv x_2\equiv c\pmod g\\
\begin{cases}
kg+c\equiv x_1\pmod A\\
kg+c\equiv x_2\pmod B\\
\end{cases}\\
令k_1=\left\lfloor\dfrac {x_1}g\right\rfloor,k_2=\left\lfloor\dfrac {x_2}g\right\rfloor\\
\begin{cases}
kg+c\equiv k_1g+c\pmod A\\
kg+c\equiv k_2g+c\pmod B\\
\end{cases}\\
\begin{cases}
kg\equiv k_1g\pmod A\\
kg\equiv k_2g\pmod B\\
\end{cases}\\
\begin{cases}
k\equiv k_1\pmod{\left\lfloor\dfrac Ag\right\rfloor}\\
k\equiv k_2\pmod{\left\lfloor\dfrac Bg\right\rfloor}\\
\end{cases}\\
这样就可以用CRT解决了
$$
$CRT$部分见博客

卡点:

C++ Code:

#include <cstdio>
long long gcd(const long long a, const long long b) {
if (!b) return a;
return gcd(b, a % b);
}
inline long long lcm(const long long a, const long long b) { return a / gcd(a, b) * b; }
inline long long mul(const long long x, const long long y, const long long mod) {
static long long res;
res = x * y - static_cast<long long> (static_cast<long double> (x) * y / mod + 0.5) * mod;
return res + (res >> 63 & mod);
}
void exgcd(const long long a, const long long b, long long &x, long long &y) {
if (!b) x = 1, y = 0;
else exgcd(b, a % b, y, x), y -= a / b * x;
}
inline long long inv(const long long a, const long long mod) {
static long long x, y;
exgcd(a, mod, x, y);
return x + (x >> 63 & mod);
}
inline long long getreduce(const long long x, const long long mod) { return x + (x >> 63 & mod); } long long CRT(const long long A, const long long mod1, const long long B, const long long mod2) {
const long long mod = mod1 * mod2, inv1 = inv(mod1, mod2);
return getreduce(mul(mul(getreduce((B - A) % mod2, mod2), inv1, mod2), mod1, mod) + A, mod);
}
long long EXCRT(const long long A, const long long mod1, const long long B, const long long mod2) {
if (A % mod2 == B) return A;
const long long GCD = gcd(mod1, mod2), t = CRT(A / GCD, mod1 / GCD, B / GCD, mod2 / GCD);
return t * GCD + (A % GCD);
} int n;
int main() {
scanf("%d", &n);
long long LCM = 1, ans = 0;
for (int i = 0; i < n; ++i) {
static long long a, b;
scanf("%lld%lld", &a, &b); b %= a;
ans = EXCRT(ans, LCM, b, a);
LCM = lcm(LCM, a);
}
printf("%lld\n", ans);
}

[洛谷P4777]【模板】扩展中国剩余定理(EXCRT)的更多相关文章

  1. [洛谷P4777] [模板] 扩展中国剩余定理

    扩展中国剩余定理,EXCRT. 题目传送门 重温一下中国剩余定理. 中国剩余定理常被用来解线性同余方程组: x≡a[1] (mod m[1]) x≡a[2] (mod m[2]) ...... x≡a ...

  2. 中国剩余定理(CRT) & 扩展中国剩余定理(ExCRT)总结

    中国剩余定理(CRT) & 扩展中国剩余定理(ExCRT)总结 标签:数学方法--数论 阅读体验:https://zybuluo.com/Junlier/note/1300035 前置浅讲 前 ...

  3. 扩展中国剩余定理 (exCRT) 的证明与练习

    原文链接https://www.cnblogs.com/zhouzhendong/p/exCRT.html 扩展中国剩余定理 (exCRT) 的证明与练习 问题模型 给定同余方程组 $$\begin{ ...

  4. 扩展中国剩余定理 (ExCRT)

    扩展中国剩余定理 (ExCRT) 学习笔记 预姿势: 扩展中国剩余定理和中国剩余定理半毛钱关系都没有 问题: 求解线性同余方程组: \[ f(n)=\begin{cases} x\equiv a_1\ ...

  5. 扩展中国剩余定理 exCRT 学习笔记

    前言 由于 \(\{\mathrm{CRT}\}\subseteq\{\mathrm{exCRT}\}\),而且 CRT 又太抽象了,所以直接学 exCRT 了. 摘自 huyufeifei 博客 这 ...

  6. [洛谷P4720] [模板] 扩展卢卡斯

    题目传送门 求组合数的时候,如果模数p是质数,可以用卢卡斯定理解决. 但是卢卡斯定理仅仅适用于p是质数的情况. 当p不是质数的时候,我们就需要用扩展卢卡斯求解. 实际上,扩展卢卡斯=快速幂+快速乘+e ...

  7. 扩展中国剩余定理(EXCRT)快速入门

    问题 传送门 看到这个问题感觉很难??? 不用怕,往下看就好啦 假如你不会CRT也没关系 EXCRT大致思路 先考虑将方程组两两联立解开,如先解第一个与第二个,再用第一个与第二个的通解来解第三个... ...

  8. [Luogu P4777] 【模板】扩展中国剩余定理(EXCRT) (扩展中国剩余定理)

    题面 传送门:洛咕 Solution 真*扩展中国剩余定理模板题.我怎么老是在做模板题啊 但是这题与之前不同的是不得不写龟速乘了. 还有两个重点 我们在求LCM的时候,记得先/gcd再去乘另外那个数, ...

  9. P4777 【模板】扩展中国剩余定理(EXCRT)/ poj2891 Strange Way to Express Integers

    P4777 [模板]扩展中国剩余定理(EXCRT) excrt模板 我们知道,crt无法处理模数不两两互质的情况 然鹅excrt可以 设当前解到第 i 个方程 设$M=\prod_{j=1}^{i-1 ...

  10. P4777 【模板】扩展中国剩余定理(EXCRT)&& EXCRT

    EXCRT 不保证模数互质 \[\begin{cases} x \equiv b_1\ ({\rm mod}\ a_1) \\ x\equiv b_2\ ({\rm mod}\ a_2) \\ ... ...

随机推荐

  1. Qt 独立运行时伴随CMD命令窗口

    用Qt写了一个小软件,在把程序release后,打包分装后,发现程序运行的时候会伴随cmd命令窗口,可把我愁怀了 不过功夫不负有心人,在老师和我网友的帮助下,终于搞完了 CONFIG:指定工程配置和编 ...

  2. 使用httpClient获取请求cookie

    package mytest; import java.util.ArrayList; import java.util.List; import org.apache.http.NameValueP ...

  3. C++clock()延时循环

    函数clock(),返回程序开始执行后所用的系统时间,但是有两个复制问题. 1.clock()返回时间的单位不一定是秒 2.该函数的返回类型在某些系统上可能是Long,也可能是unsigned lon ...

  4. 1208: [HNOI2004]宠物收养所

    1208: [HNOI2004]宠物收养所 Time Limit: 10 Sec Memory Limit: 162 MB Submit: 12030 Solved: 4916 Description ...

  5. Kali信息收集-DNS

    1.whois查询 直接在终端输入whois 域名 2.查找dns服务器 (1)host (2)dig (3)nslookup 3.域传输 4.域名枚举 (1)dnsdict6 kali没有集成这款工 ...

  6. NTP错误总结

    Ntp错误总结 解决ntp的错误 no server suitable for synchronization found 当用ntpdate -d 来查询时会发现导致 no server suita ...

  7. POJ 2631 Roads in the North(求树的直径,两次遍历 or 树DP)

    题目链接:http://poj.org/problem?id=2631 Description Building and maintaining roads among communities in ...

  8. POJ 3348 Cows(凸包+多边形面积)

    Description Your friend to the south is interested in building fences and turning plowshares into sw ...

  9. C++多态实现与继承

    面向对象的三个基本特征 面向对象的三个基本特征是:封装.继承.多态.其中, 封装可以隐藏实现细节,使得代码模块化: 继承可以扩展已存在的代码模块(类),它们的目的都是为了——代码重用: 而多态则是为了 ...

  10. activemq控制面板里的NumberOfPendingMessages、MessagesEnqueued、MessagesDequeued含义

    Number Of Consumers 消费者 这个是消费者端的消费者数量.Number Of Pending Messages 等待消费的消息 这个是当前未出队列的数量.可以理解为总接收数-总出队列 ...