码迷,mamicode.com
首页 > 其他好文 > 详细

poj2891--Strange Way to Express Integers(不互素的中国剩余定理)

时间:2015-08-05 10:39:54      阅读:128      评论:0      收藏:0      [点我收藏+]

标签:

题目链接:点击打开链接

题目大意:有一个数x,x%ai = ri ,给出n对ai和ri,问x的最小非负整数是什么,如果不存在输出-1

不互素的中国剩余定理:

x%a1= r1 ; x%a2 = r2 ; 设k1,k2得到x = a1*k1 + r1 , x = a2*k2+r2 

那么a1*k1+r1 = a2*k2+r2 --> a1*k1 = (r2-r1) + a2*k2---->对整个式子进行a2取余,得到(a1*k1)%a2 = (r2-r1)%a2,这里面只有一个未知量k1,用扩展gcd求出k1,计算出x = a1*k1+r1

这个x只是满足x%a1= r1 ; x%a2 = r2 ;不一定满足x%a3 = r3,所以求出的x只是真正的解的一个特解,真正的解围ans

那么ans ≡ x%lca(a1,a2)  --> 其中只有ans是一个未知数,所以就将两个式子转换为了一个式子,这样不断的合并,求最终的结果。

判断是不是有解,在扩展gcd中求k1时判断。


#include <cstdio>
#include <cstring>
#include <algorithm>
#define LL __int64
using namespace std;
void gcd(LL a,LL b,LL &d,LL &x,LL &y) {
    if(b == 0) {
        d = a ;
        x = 1 ;
        y = 0 ;
    }
    else {
        gcd(b,a%b,d,y,x);
        y = (a/b)*x ;
    }
    return ;
}
int main()
{
    LL k , a1 , a2 , r1 , r2 , d , x , y ;
    while(scanf("%I64d", &k)!=EOF)
    {
        LL flag = 1 ;
        scanf("%I64d %I64d", &a1, &r1);
        k-- ;
        while(k--)
        {
            scanf("%I64d %I64d", &a2, &r2);
            gcd(a1,a2,d,x,y);
            if( (r2-r1)%d )
                flag = 0 ;
            if( flag )
            {
                x = (r2-r1)/d*x ;
                y = a2/d ;
                x = ( x%y +y)%y ;
                r1 = x*a1 + r1 ;
                a1 = (a1*a2)/d ;
            }
        }
        gcd(1,a1,d,x,y);
        if( r1%d )
            flag = 0 ;
        if(flag == 0)
            printf("-1\n");
        else
        {
            x = r1/d*x ;
            y = a1 / d ;
            x = ( x%y+y )%y ;
            printf("%I64d\n", x);
        }
    }
    return 0;
}


版权声明:本文为博主原创文章,未经博主允许不得转载。

poj2891--Strange Way to Express Integers(不互素的中国剩余定理)

标签:

原文地址:http://blog.csdn.net/winddreams/article/details/47291795

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!