【poj2891】 Strange Way to Express Integers

poj.org/problem?id=2891 (题目链接)

题意：求解线性同余方程组，不保证模数一定两两互质。

Solotion

　　用exgcd将俩个同余方程合并成一个
　　
　　如合并n%M=R,n%m=r
　　
　　即M*x+R=m*y+r
　　
　　M*x-m*y=r-R
　　
　　设a=M/t,b=m/t,c=(r-R)/t，t=gcd(a,b)
　　
　　若(r-R)%t!=0则无解
　　
　　用exgcd得到a*x+b*y=c的解x0，
　　
　　通解x=x0+k*b，k为整数
　　
　　带入M*x+R=n
　　
　　M*b*k+M*x0+R=n
　　
　　所以R=R+M*x0，M=M*b
　　

　　蒯自hzwer。
　　注意当最后发现方程无解直接退出时，会导致有数据没有读完，然后就会Re，所以现用数组将所有数据存下来。

代码：

// poj2891
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<cmath>
#define LL long long
#define inf 2147483640
#define Pi 3.1415926535898
#define free(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
using namespace std;const int maxn=100010;
LL mm[maxn],rr[maxn],n;LL gcd(LL a,LL b) {return a%b==0 ? b : gcd(b,a%b);
}
void exgcd(LL a,LL b,LL &x,LL &y) {if (!b) {x=1,y=0;return;}exgcd(b,a%b,y,x);y-=x*(a/b);   //注意此处a/b一定要到打括号，因为取的是a/b的整数部分
}
int main() {while (scanf("%lld",&n)!=EOF) {LL M,R,flag=1;scanf("%lld%lld",&M,&R);for (int i=1;i<n;i++) scanf("%lld%lld",&mm[i],&rr[i]);for (int i=1;i<n;i++) {LL m,r,x,y;m=mm[i],r=rr[i];LL d=gcd(M,m);if ((r-R)%d!=0) {printf("-1\n");flag=0;break;}exgcd(M/d,m/d,x,y);   //x*(M/d)+y*(m/d)=d,解出x，yx=((r-R)/d*x%(m/d)+(m/d))%(m/d);   //求出最小xR+=M*x;M*=m/d;   //lcm(M,m)}if (flag) printf("%lld\n",R);}return 0;
}

【poj2891】 Strange Way to Express Integers相关推荐

【poj 2891】Strange Way to Express Integers（数论--拓展欧几里德求解同余方程组模版题）...
题意:Elina看一本刘汝佳的书(O_O*),里面介绍了一种奇怪的方法表示一个非负整数 m .也就是有 k 对 ( ai , ri ) 可以这样表示--m%ai=ri.问 m 的最小值. 解法:拓展欧 ...
数论练习1 （曹冲养猪 + [POJ 2891]Strange Way to Express Integers + 乘法逆元【带证明】）
虽然作业还没有做完,但是我还是放不下它,对此,我只想说: 今天你对作业爱理不理,明天它就让你补到飞起 DP先放放,我们要雨露均沾练习上手:乘法逆元题目题解代码实现曹冲养猪?(互质的中国剩余定 ...
POJ2891 Strange Way to Express Integers【扩展中国剩余定理】
题目大意就是模板...没啥好说的思路因为模数不互质,所以直接中国剩余定理肯定是不对的然后就考虑怎么合并两个同余方程 \(ans = a_1 + x_1 * m_1 = a_2 + x_2 * ...
[POJ2891] Strange Way to Express Integers
题目描述 FJ正在读佳佳写的一本书,书中描述一种表示非负整数的方法:选择k个不同的正整数a1,a2,-,ak,对于某个整数m分别对ai求余对应整数ri,如果适当选择a1,a2,-,ak,那么整数m可由 ...
poj-2891(Strange Way to Express Integers)--中国剩余定理扩展欧几里得
题意:找到一个m,使得m%ai=ri,并且这个m最小 m = a1*x + r1; m = a2*y + r2; 可得:a1*x - a2*y = r2 - r1 可得:a1*x ≡ r2-r1(mo ...
POJ2891 Strange Way to Express Integers (扩展欧几里德)
本文为博主原创文章,欢迎转载,请注明出处 www.cnblogs.com/yangyaojia 题目大意求解一组同余方程 x ≡ r1 (mod a1) x ≡ r2 (mod a2) x ≡ r ...
POJ2891 Strange Way to Express Integers
线性同余方程题目传送门解线性同余方程的简单应用,注意要把数据读完再输出-1. 代码: #include<cstdio> #include<cstring> #include ...
POJ2891——Strange Way to Express Integers（一元线性同余方程组）
由若干个一元线性同余方程组构成的方程组,叫做一元线同余方程组. 求解我们可以将其统一划成a*x ≡ b(mod m)的形式,这样有利于算法的实现. 对于同于方程组的求解,其实质过程就是对于其中的同余 ...
poj 2891 Strange Way to Express Integers 2012-09-05
http://poj.org/problem?id=2891 解线性模方程组. 比较坑爹,数据比较大,很容易溢出. 1 Program poj2891; 2 3 var m:int64; 4 5 a, ...

【poj2891】 Strange Way to Express Integers

【poj2891】 Strange Way to Express Integers相关推荐

最新文章

热门文章