本文为博主原创文章,欢迎转载,请注明出处 www.cnblogs.com/yangyaojia

题目大意

求解一组同余方程

x ≡ r1 (mod a1)

x ≡ r2 (mod a2)

x ≡ r3 (mod a3)

......

x ≡ rk (mod ak)

的解x(a1,a2,a3,.....ak 并不一定互质)。如果不存在则输出-1.

输入格式

有多组数据,每组数组第一行为k,后面有k行,每行两个数,代表ai,ri。

输出格式

每一行对应每一个询问的解x。

样例输入

2

8 7

11 9

样例输出

31

分析

看到这一题,可以发现好像问的就是中国剩余定理的问题,可是题目中a并不互质,无法用中国剩余定理来解决。不过我们可以从扩展欧几里德算法入手。

先来分析规模小一点的。

对于一组同余方程

\(\begin {cases} x \mod a_1=r_1\\ x\mod a_2=r_2\end{cases}\to\begin{cases}k_1*a_1+r_1=x\\ k_2* a_2+r_2=x\end {cases}\)

我们可以上下两式相减得出

$ k_1* a_1-k_2* a_2=r_2-r_1$

我们可以发现,这个式子可以用扩展欧几里德来求解。于是我们求出了\(k_1\)

便可以将其带入原来的式子,求出\(x=k_1*a_1+r_1\)。这是两个方程的求解,面对多个方程,我们可以这样做:

假设我们刚才求出的是\(x_1\),而为了求出满足三个方程的解\(x_2\)我们可以得到这样的式子

\(\begin {cases} x_2\mod lcm(a_1,a_2)=x_1\\ x_2\mod a_3=r_3\end{cases}\)

证明:

\(x_2=k_3*lcm(a_1,a_2)+x_1\) ,因为\(k_3*lcm(a_1,a_2)\)这一块肯定能整除\(a_1\)或\(a_2\),所以肯定会剩下一个\(c_1\),而\(c_1\)恰好满足前面的方程

用同样的方法,我们可以求出满足所有式子的解。

#include <cstdio>
#include <iostream>
#include <cmath>
#include <queue>
#include <algorithm>
#include <cstring>
#include <climits>
#define MAXN 10000+10
#define abs(a) a>0?a:-a
using namespace std;
long long e_gcd(long long a,long long b,long long& x,long long& y)
{
if(!b)
{
x=1;y=0;
return a;
}
long long ans=e_gcd(b,a%b,x,y);
long long tmp=x;
x=y;
y=tmp-a/b*y;
return ans;
}
long long a1,a2,b1,b2,x,y,gcd,t,c;
int main()
{
while(scanf("%lld",&t)!=EOF)
{
scanf("%lld%lld",&a1,&b1);
int flag=0;
if(t==1)
printf("%lld\n",b1);
else
{
for(int i=2;i<=t;i++)
{
scanf("%lld%lld",&a2,&b2);
if(flag) continue;
gcd=e_gcd(a1,a2,x,y);
if((b2-b1)%gcd!=0) {flag=1;}
x*=(b2-b1)/gcd;
x%=a2/gcd;
if(x<0) x+=abs(a2/gcd);
c=a1*x+b1;
a1=a1/gcd*a2; b1=c;
}
}
if(flag==1) {printf("-1\n");continue;}
printf("%lld\n",c); } return 0;
}

POJ2891 Strange Way to Express Integers (扩展欧几里德)的更多相关文章

  1. POJ2891 Strange Way to Express Integers 扩展欧几里德 中国剩余定理

    欢迎访问~原文出处——博客园-zhouzhendong 去博客园看该题解 题目传送门 - POJ2891 题意概括 给出k个同余方程组:x mod ai = ri.求x的最小正值.如果不存在这样的x, ...

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

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

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

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

  4. POJ2891——Strange Way to Express Integers(模线性方程组)

    Strange Way to Express Integers DescriptionElina is reading a book written by Rujia Liu, which intro ...

  5. POJ2891 Strange Way to Express Integers

    题意 Language:Default Strange Way to Express Integers Time Limit: 1000MS Memory Limit: 131072K Total S ...

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

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

  7. POJ2891 Strange Way to Express Integers【扩展中国剩余定理】

    题目大意 就是模板...没啥好说的 思路 因为模数不互质,所以直接中国剩余定理肯定是不对的 然后就考虑怎么合并两个同余方程 \(ans = a_1 + x_1 * m_1 = a_2 + x_2 * ...

  8. POJ.2891.Strange Way to Express Integers(扩展CRT)

    题目链接 扩展中国剩余定理:1(直观的).2(详细证明). [Upd:]https://www.luogu.org/problemnew/solution/P4774 #include <cst ...

  9. POJ - 2891 Strange Way to Express Integers (扩展中国剩余定理)

    题目链接 扩展CRT模板题,原理及证明见传送门(引用) #include<cstdio> #include<algorithm> using namespace std; ty ...

随机推荐

  1. jquery日历插件FullCalendar使用技巧

    原文链接:http://blog.csdn.net/u013493957/article/details/44920341   FullCalendar是一款基于jquery的日历控件,它有着很强大的 ...

  2. oracle 解除锁表sql

    select b.owner,b.object_name,a.session_id,a.locked_mode from v$locked_object a,dba_objects b where b ...

  3. spring的启动过程就是创建ioc容器的过程

    1. spring简介 spring的最基本的功能就是创建对象及管理这些对象之间的依赖关系,实现低耦合.高内聚.还提供像通用日志记录.性能统计.安全控制.异常处理等面向切面的能力,还能帮我们管理最头疼 ...

  4. POJ 1944

    明天补上... 这道题的思路确实很精致.考虑到连的边肯定不会是一个环,所以至少有一个断点.于是,可以枚举这个断点.断点一确定,那么连边的走向也就确定了.用D[i]表示由i开始可以到达的最远点即可.对于 ...

  5. Cocos2D-x设计模式发掘之中的一个:单例模式

    http://www.tuicool.com/articles/NBRn2murl=pVtZACoQFKXC3u3uGwMLnTy4YDWihcVg0ata5gy506pmPpQEc0PO9hm6wG ...

  6. 更改python字符编码以便使用UTF-8的编码url路径

    url编码分两种, 一种是unicode, 另一种是gb2312, 今天遇到的一个网站是要将字符编码按照gb2312来编码,用来得到一个先填写blanks后再返回页面的数据,废话少说,需要做的就是先查 ...

  7. Caffe solver.prototxt学习

    在solver解决下面的四个问题: a.训练的记录(bookkeeping),创建用于training以及test的网络结构: b.使用前向以及反向过程对training网络参数学习的过程: c.对t ...

  8. spring《四》自动装配

    byName模式<bean autowire="byName"> Spring会查找一个叫做date的bean定义. byType模式<bean autowire ...

  9. Re:从0开始的微服务架构--(二)快速快速体验微服务架构?--转

    原文地址:https://mp.weixin.qq.com/s/QO1QDQWnjHZp8EvGDrxZvw 这是专题的第二篇文章,看看如何搭建一个简单模式的微服务架构. 记得好久之前看到一个大牛说过 ...

  10. 入门 IT 行业,该具备哪些技能?

    对于刚开始进入IT的新人来说,“必备技能”往往意味着一个长长的.标有重要度的学习列表,但是过长的列表通常会导致新人不知如何开始学习,压力倍增.本文尝试列举出最重要的几个技能,也期望通过此列表能给新人一 ...