题目链接:http://www.51nod.com/Challenge/Problem.html#!#problemId=1079

一个正整数K,给出K Mod 一些质数的结果,求符合条件的最小的K。例如,K % 2 = 1, K % 3 = 2, K % 5 = 3。符合条件的最小的K = 23。

收起

 

输入

第1行:1个数N表示后面输入的质数及模的数量。(2 <= N <= 10)
第2 - N + 1行,每行2个数P和M,中间用空格分隔,P是质数,M是K % P的结果。(2 <= P <= 100, 0 <= K < P)

输出

输出符合条件的最小的K。数据中所有K均小于10^9。

输入样例

3
2 1
3 2
5 3

输出样例

23

解题思路:中国剩余定理模板题
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
ll n,a[],m[];
//ax+by=gcd(a,b);
//x=y1,y=x1-a/b*y1;
void exgcd(ll a,ll b,ll &x,ll &y,ll &c)
{
if(!b){
x=; y=; c=a;
return;
}
exgcd(b,a%b,y,x,c);
y-=a/b*x;
}
ll China()
{
ll x,y,c,lcm=,ans=;
for(int i=;i<=n;i++)
lcm*=m[i];
for(int i=;i<=n;i++)
{
exgcd(lcm/m[i],m[i],x,y,c);
x=(x%m[i]+m[i])%m[i];
ans=(ans+x*lcm/m[i]*a[i])%lcm;
}
return (ans+lcm)%lcm;
}
int main()
{
while(cin>>n){
for(int i=;i<=n;i++)cin>>m[i]>>a[i];
cout<<China()<<endl;
}
return ;
}

题目链接:http://poj.org/problem?id=2891

给定 2n2n 个正整数 a_1,a_2,\cdots ,a_na1​,a2​,⋯,an​ 和 m_1,m_2,\cdots ,m_nm1​,m2​,⋯,mn​,求一个最小的正整数 xx,满足 \forall i\in[1,n],x\equiv a_i\ (\bmod m_i\ )∀i∈[1,n],x≡ai​ (modmi​ ),或者给出无解。

输入格式

多组数据。

每组数据第一行一个整数 nn;
接下来 nn 行,每行两个整数 m_i,a_imi​,ai​。

输出格式

对于每组数据,若无解,输出 -1−1;否则输出一个非负整数,若有多解,输出最小的满足条件的答案。

样例

样例输入

2
8 7
11 9

样例输出

31

数据范围与提示

对于全部数据,所有的输入都是非负的,并且可以用 6464 位有符号整数表示。保证 1\le n\le 10^5,m_i\gt a_i1≤n≤105,mi​>ai​。

解题思路:这种是一般情形,需要用扩展中国剩余定理。

代码:

#include<iostream>
using namespace std;
typedef long long ll;
ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
ll n,a[],m[];
//ax+by=gcd(a,b);
//x=y1,y=x1-a/b*y1;
void exgcd(ll a,ll b,ll &x,ll &y,ll &c)
{
if(!b){
x=; y=; c=a;
return;
}
exgcd(b,a%b,y,x,c);
y-=a/b*x;
}
ll inv(ll a,ll b)
{
ll x,y,c;
exgcd(a,b,x,y,c);
x=(x%(b/c)+(b/c))%(b/c);
return x;
}
ll exCRT()
{
for(int i=;i<=n;i++)
{
ll m1=m[i-],m2=m[i],a1=a[i-],a2=a[i],c=gcd(m1,m2);
if((a2-a1)%c!=)return -;
m[i]=m1*m2/c;
a[i]=(inv(m1/c,m2/c)*(a2-a1)/c)%(m2/c)*m1+a1;
a[i]=(a[i]%m[i]+m[i])%m[i];
}
return a[n];
}
int main()
{
while(cin>>n){
for(int i=;i<=n;i++)cin>>m[i]>>a[i];
cout<<exCRT()<<endl;
}
return ;
}

51nod1079 poj2891 中国剩余定理与其扩展的更多相关文章

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

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

  2. CRT&EXCRT 中国剩余定理及其扩展

    前言: 中国剩余定理又名孙子定理.因孙子二字歧义,常以段子形式广泛流传. 中国剩余定理并不是很好理解,我也理解了很多次. CRT 中国剩余定理 中国剩余定理,就是一个解同余方程组的算法. 求满足n个条 ...

  3. 扩展中国剩余定理(扩展CRT)详解

    今天在$xsy$上翻题翻到了一道扩展CRT的题,就顺便重温了下(扩展CRT模板也在里面) 中国剩余定理是用于求一个最小的$x$,满足$x\equiv c_i \pmod{m_i}$. 正常的$CRT$ ...

  4. HDU 5768Lucky7(多校第四场)容斥+中国剩余定理(扩展欧几里德求逆元的)+快速乘法

    地址:http://acm.hdu.edu.cn/showproblem.php?pid=5768 Lucky7 Time Limit: 2000/1000 MS (Java/Others)    M ...

  5. LOJ2721 [NOI2018] 屠龙勇士 【扩展中国剩余定理】

    好久没写了,写一篇凑个数. 题目分析: 这题不难想,讲一下中国剩余定理怎么扩展. 考虑$$\left\{\begin{matrix}x \equiv a\pmod{b}\\ x \equiv c\pm ...

  6. 学习笔记 - 中国剩余定理&扩展中国剩余定理

    中国剩余定理&扩展中国剩余定理 NOIP考完回机房填坑 ◌ 中国剩余定理 处理一类相较扩展中国剩余定理更特殊的问题: 在这里要求 对于任意i,j(i≠j),gcd(mi,mj)=1 (就是互素 ...

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

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

  8. 中国剩余定理+扩展中国剩余定理 讲解+例题(HDU1370 Biorhythms + POJ2891 Strange Way to Express Integers)

    0.引子 每一个讲中国剩余定理的人,都会从孙子的一道例题讲起 有物不知其数,三三数之剩二,五五数之剩三,七七数之剩二.问物几何? 1.中国剩余定理 引子里的例题实际上是求一个最小的x满足 关键是,其中 ...

  9. [poj2891]Strange Way to Express Integers(扩展中国剩余定理)

    题意:求解一般模线性同余方程组 解题关键:扩展中国剩余定理求解.两两求解. $\left\{ {\begin{array}{*{20}{l}}{x = {r_1}\,\bmod \,{m_1}}\\{ ...

随机推荐

  1. python3 selenium webdriver 元素定位xpath定位骚操作

    源文http://www.cnblogs.com/qingchunjun/p/4208159.html By.xpath() 这个方法是非常强大的元素查找方式,使用这种方法几乎可以定位到页面上的任意元 ...

  2. MyEclipse10 复制之前的项目部署到tomcat时项目名称对不上,还是复制前的项目名称,哪里修改设置

    工程 -- 右键属性 -- Myeclispse -- web修改一下发布名字就可以了.

  3. bootstrap模态框关闭后清除模态框的数据

    https://segmentfault.com/q/1010000008789123 bootstrap模态框第二次打开时如何清除之前的数据? 我用了bootstrap模态框的remote功能,在弹 ...

  4. 版本控制--git+idea

  5. IntelliJ IDEA -- 破解

    ① 到这个地方下载 IntelliJ IDEA 注册码:http://idea.lanyus.com/  就是这个jar包:JetbrainsCrack-2.6.10-release-enc.jar ...

  6. Mysql Router 的集群

    1. c:\mysql-router, c:\mysql-5.7.23, 这两个目录的bin都要加入path 2. c:\mysql-shell,在bin下,有一个 mysqlsh.exe, 双击,打 ...

  7. Python——进程通信之间数据共享

    from multiprocessing import Manager,Process,Lock def main(dic,lock): lock.acquire() dic['count'] -= ...

  8. css last

  9. mybatis-spring-1.2.2.jar下载地址

    http://www.java2s.com/Code/Jar/m/Downloadmybatisspring120jar.htm

  10. codeforces347B

    Fixed Points CodeForces - 347B A permutation of length n is an integer sequence such that each integ ...