POJ 2891

线性同余方程组,模板了。但要注意读完数据才跳出循环啊

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

void exgcd(__int64 a,__int64 b,__int64 &d,__int64 &xx,__int64 &yy){
	if(b==0){
		xx=1;yy=0;d=a;
		return ;
	}
	else{
		exgcd(b,a%b,d,xx,yy);
		__int64 tmp=xx;
		xx=yy;
		yy=tmp-(a/b)*yy;
	}
}

int main(){
	__int64 n;
	__int64 a1,r1,a2,r2; __int64 d,xx,yy;
	bool ifhave;
	while(scanf("%I64d",&n)!=EOF){
		scanf("%I64d%I64d",&a1,&r1);
		ifhave=true;
		for(__int64 i=1;i<n;i++){
			scanf("%I64d%I64d",&a2,&r2);
			if(!ifhave) continue;
			__int64 a=a1,b=a2,c=r2-r1;
			exgcd(a,b,d,xx,yy);
			if(c%d!=0){
				ifhave=false;
			}
			__int64 t=b/d;
			xx=(xx*(c/d)%t+t)%t;
			r1=a1*xx+r1;
			a1=a1*(a2/d);
		}
		if(!ifhave){
			printf("-1\n");
		}
		else printf("%I64d\n",r1);
	}
	return 0;
}

  

时间: 2024-10-31 02:51:39

POJ 2891的相关文章

poj 2891 Strange Way to Express Integers

http://poj.org/problem?id=2891 这道题的题意是:给你多个模性方程组:m mod ai=ri 求最小的m: 中国剩余定理 1 #include <cstdio> 2 #include <cstring> 3 #include <algorithm> 4 #define ll long long 5 using namespace std; 6 7 ll gcd(ll a,ll b,ll &x,ll &y) 8 { 9 if(!

poj 2891 Strange Way to Express Integers (扩展gcd)

题目链接 题意:给k对数,每对ai, ri.求一个最小的m值,令m%ai = ri; 分析:由于ai并不是两两互质的, 所以不能用中国剩余定理. 只能两个两个的求. a1*x+r1=m=a2*y+r2联立得:a1*x-a2*y=r2-r1;设r=r2-r2; 互质的模线性方程组m=r[i](mod a[i]).两个方程可以合并为一个,新的a1为lcm(a1,a2), 新的r为关于当前两个方程的解m,然后再和下一个方程合并--.(r2-r1)不能被gcd(a1,a2)整除时无解. 怎么推出的看了好

poj 2891 Strange Way to Express Integers (非互质的中国剩余定理)

Strange Way to Express Integers Time Limit: 1000MS   Memory Limit: 131072K Total Submissions: 9472   Accepted: 2873 Description Elina is reading a book written by Rujia Liu, which introduces a strange way to express non-negative integers. The way is

poj 2891 Strange Way to Express Integers (解模线性方程组)

链接:poj 2891 题意:有一个数x,给定k组ai和ri,使得x%ai=ri 求x最小为多少 分析:求解模线性方程组 x = a1(mod m1) x = a2(mod m2) x = a3(mod m3) 先求解方程组前两项. x=m1*k1+a1=m2*k2+a2 -> m1*k1+m2*(-k2)=a2-a1 这个方程可以通过欧几里得求解出最小正整数的k1 则x=m1*k1+a1 显然x为两个方程的最小正整数解. 则这两个方程的通解为 X=x+k*LCM(m1,m2) -> X=x(

【POJ 2891】Strange Way to Express Integers(扩展欧几里得)

[POJ 2891]Strange Way to Express Integers(扩展欧几里得) Time Limit: 1000MS   Memory Limit: 131072K Total Submissions: 12934   Accepted: 4130 Description Elina is reading a book written by Rujia Liu, which introduces a strange way to express non-negative in

POJ 2891 Strange Way to Express Integers 中国剩余定理MOD不互质数字方法

http://poj.org/problem?id=2891 711323 97935537 475421538 1090116118 2032082 120922929 951016541 15898 418373 161478614 149488440 1748022751 21618619576 810918992 241779667 1772616743 1953316358 125248280 2273149397 3849022001 2509433771 3885219405 35

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,r:array[1..30000000]of int64; 6 7 ans,x,y,lcm:int64; 8 9 10 Procedure init; 11 12 var i,j:longint; 13 14 begin 15 16 m:=0; 17 18 readln(m); 19 20 for

poj 2891 Strange Way to Express Integers(中国剩余定理)

http://poj.org/problem?id=2891 题意:求解一个数x使得 x%8 = 7,x%11 = 9; 若x存在,输出最小整数解.否则输出-1: ps: 思路:这不是简单的中国剩余定理问题,由于输入的ai不一定两两互质,而中国剩余定理的条件是除数两两互质. 这是一般的模线性方程组,对于 X mod m1=r1 X mod m2=r2 ... ... ... X mod mn=rn 首先,我们看两个式子的情况 X mod m1=r1-----------------------(

poj——2891 Strange Way to Express Integers

Strange Way to Express Integers Time Limit: 1000MS   Memory Limit: 131072K Total Submissions: 16839   Accepted: 5625 Description Elina is reading a book written by Rujia Liu, which introduces a strange way to express non-negative integers. The way is

POJ 2891 中国剩余定理(不互素)

Strange Way to Express Integers Time Limit: 1000MS   Memory Limit: 131072K Total Submissions: 17877   Accepted: 6021 Description Elina is reading a book written by Rujia Liu, which introduces a strange way to express non-negative integers. The way is