POJ 2891 扩展欧几里德

这个题乍一看跟剩余定理似的,但是它不满足两两互素的条件,所以不能用剩余定理,也是给了一组同余方程,找出一个X满足这些方程,如果找不到的话就输出-1

因为它不满足互素的条件,所以两个两个的合并,最后合成一个。

题目给定的是

M % m1 = r1

M % m2 = r2

......

M % mn = rn

只需将两个式子合并成一个式子,那么这个合并的这个式子就可以继续和下面的式子继续合并,知道合到最后一个式子。

首先来看下两个式子怎么合并。

M % m1 = r1    可以写成  M = k1 * m1 + r1;

同理M =  k2 * m2 + r2

所以k1 * m1 + r1  = k2 * m2 + r2

k1 * m1 - k2 * m2 = r2 - r1

将k1 看成x, m1看成a, k2看成y, m2看成b, r2 - r1看成c

那么这个式子就变成了

ax-by=c;不定方程

这时候就可以用不定方程来解了。

因为题目给的两两可能不互素,所以gcd(a, b) 有可能不等于1,所以,这个方程不一定有解,所以就是判断这个题有没有解的条件。

如果有解的话,用扩展欧几里德可以的出这个方程的解,然后将x带入到 k1 * m1 + r1中得到M,这时候M就是这两个方程组的一个特解,通解就是M‘ = M + K * LCM(m1, m2);

这个式子可以变化成M‘ % LCM(m1, m2) = M;所以又可以继续跟下面的式子合并了。所以这个式子对应的余数r就是M, 模数(就是式子中的mi)就是LCM(m1,m2);

到这里之后就合并好了一个式子。下面就是继续按照这个过程合并到第n个式子了。最后合并完之后就剩下一个式子,求出来的M就是方程组的特解,通解为M‘ = M + k * LCM; 其中LCM = LCM(m1, m2...mn), 最小正整数解为(M %  LCM + LCM) % LCM;

还有一个坑就是如果边输入边处理的话,不要找到无解立马就跳出循环,因为后面的数据还没输入完。

代码如下:

    /*************************************************************************
        > File Name: 5.cpp
        > Author: Howe Young
        > Mail: [email protected]
        > Created Time: 2015年08月01日 星期六 15时39分47秒
     ************************************************************************/

    #include <iostream>
    #include <cstdio>
    using namespace std;
    const int maxn = 1e5 + 5;
    typedef long long ll;
    ll m[maxn], r[maxn];
    ll exgcd(ll a, ll b, ll &x, ll &y)
    {
        if (b == 0)
        {
            x = 1;
            y = 0;
            return a;
        }
        ll r = exgcd(b, a % b, x, y);
        ll t = x;
        x = y;
        y = t - a / b * y;
        return r;
    }
    int main()
    {
        int k;
        while (~scanf("%d", &k))
        {
            for (int i = 1; i <= k; i++)
                scanf("%lld %lld", &m[i], &r[i]);
            ll m1, m2, r1, r2, c, d;
            m1 = m[1]; r1 = r[1];
            bool flag = true;
            for (int i = 2; i <= k; i++)
            {
                ll x, y;
                d = exgcd(m1, m[i], x, y);
                if ((r[i] - r1) % d != 0)
                {
                    flag = false;
                    break;
                }
                c = r[i] - r1;
                x = c / d * x % (m[i] / d);
                r1 = m1 * x + r1;
                m1 = m1 / d * m[i];
                r1 %= m1;
            }
            if (!flag)
            {
                printf("-1\n");
                continue;
            }
            r1 = (r1 + m1) % m1;
            printf("%lld\n", r1);
        }

        return 0;
    }

时间: 2024-10-21 09:48:04

POJ 2891 扩展欧几里德的相关文章

poj 1061 扩展欧几里德同余方程求解

摘要写在一瞪眼. #include<iostream> using namespace std; long long exgcd(long long a,long long b,long long &k,long long &t) { if (b==0) { k=1; t=0; return a; } else { long long tp_gcd; tp_gcd=exgcd(b,a%b,k,t); long long temp; temp=k; k=t; t=temp-(a/

数学#扩展欧几里德 POJ 1061&amp;2115&amp;2891

寒假做的题了,先贴那时写的代码. POJ 1061 #include<iostream> #include<cstdio> typedef long long LL; using namespace std; void extend_gcd(LL a,LL b,LL &d,LL &x,LL &y) { if(b==0) { d=a; x=1,y=0; } else { extend_gcd(b,a%b,d,y,x); y-=x*(a/b); } } int

POJ 2891-Strange Way to Express Integers(扩展欧几里德)

题目地址:POJ 2891 题意:给你k组同余关系,每组包含一个ai和ri,让你找出一个最小的数m,满足m%a1=r1,m%a2=r2.......m%ak=rk. 思路:纵观上述公式,很熟悉,其实就是求两两公式之间的最小值,例如K=3,那么先求第一组和第二组的最小,然后合并第一组和第二组,然后用合并之后的再和第三组找最小,最后的结果就是最终的结果.也就是这个题分两部分来完成. 1.找出两组最小.对于m%a1=r1和m%a2=r2可以得出两个公式m=a1*x+r1,m=a2*y+r2(x,y相当

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 2142-The Balance(扩展欧几里德)

题目地址:POJ 2142 题意:有两种类型的砝码质量分别为a和b,要求称出质量为d的物品,要求a的数量x和b的数量y的和x+y最小,若有多个x+y的值,取ax+by最小的. 思路:我们应该求ax+by=d.这里我们应用扩展欧几里德求出ax+by=gcd(a,b),那么ax/gcd(a,b)+by/gcd(a,b)=1,然后求出来特解,令x=x*n,把x转化为最小正值,即x=(x%b+b)%b,求出此时的y=(d-ax)/b,若求出的y是负值,把y变成正的,因为砝码的位置涉及左右之分.同理求出y

POJ 1061青蛙的约会(扩展欧几里德)

对欧几里德不太熟悉,参考了网上的一些讲解又学习了一下 利用扩展欧几里德算法求线性方程的一般过程:a*x + b*y = m 令a1 = a/gcd(a,b) b1 = b/gcd(a,b) m1 = m/gcd(a,b) a*x + b*y = m两边同除以m1a*x/m1 + b*y/m1 = m/m1 = gcd(a,b)设x1 = x/m1 ,y1 = y/m1 则原式变为a*x1 + b*y1 = gcd(a,b)若求出这个方程中的x1,y1,那么x = x1*m1, y = y1*m1

poj 1061 青蛙约会(扩展欧几里德)

题目链接: http://poj.org/problem?id=1061 题目大意: 中文题目,题意一目了然,就是数据范围大的出奇. 解题思路: 假设两只青蛙都跳了T次,可以列出来不定方程:p*l + (n-m)*T == x - y.列出等式以后,利用扩展欧几里德计算不定方程的解.在求出整数最小解的地方卡了好久,好久. 想具体了解扩展欧几里德的用法和证明的话,可以看一下神牛的博文,我自认弱绞尽脑汁也写不来这么好,附上链接:http://www.cnblogs.com/frog112111/ar

【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 2142 The Balance【扩展欧几里德】

题意:有两种类型的砝码,每种的砝码质量a和b给你,现在要求称出质量为c的物品,要求a的数量x和b的数量y最小,以及x+y的值最小. 用扩展欧几里德求ax+by=c,求出ax+by=1的一组通解,求出当x取最小合法正整数解时y的取值,当y小于0时,说明应该放在a的另一边,变为正值.同理当y取最小时,可得到另一组解,比较两组解,取最小即可. #include<stdio.h> int ex_gcd(int a,int b,int &x,int &y){ if(!b){ x=1,y=