[codevs3990]中国剩余定理（rt）

最新推荐文章于 2018-03-13 11:13:29 发布

wwyx2001

最新推荐文章于 2018-03-13 11:13:29 发布

阅读量296

点赞数 1

分类专栏：中国剩余定理

本文链接：https://blog.csdn.net/Blue_CuSO4/article/details/78190814

版权

中国剩余定理专栏收录该内容

4 篇文章 0 订阅

订阅专栏

题目：

我是超链接

题解：

裸题？下面讲解见咯

代码：

#include <cmath>
#include <cstdio>
#define LL long long
using namespace std;
int m[12],c[12];
void exgcd(LL a,LL b,LL &x,LL &y){if (!b) x=1,y=0;else exgcd(b,a%b,y,x),y-=(a/b)*x;}
int main()
{
    LL M=1,ans=0,a,b;
    int n,i;
    scanf("%d%lld%lld",&n,&a,&b);
    for (i=1;i<=n;i++)
    {
        scanf("%d%d",&m[i],&c[i]);
        M=M*(LL)m[i];
    }
    for (i=1;i<=n;i++)
    {
        LL x,y;
        exgcd(M/m[i],m[i],x,y);
        x=(x%m[i]+m[i])%m[i];
        if (!x) x+=m[i];
        LL ni=M/m[i]*x;
        ans+=(LL)c[i]*ni;
    }
    ans%=M;if (!ans) ans+=M;    
    if (ans>b || (ans<a && (ceil((double)(a-ans)/(double)M)>floor((double)(b-ans)/(double)M)) ) ) printf("0\n0");
    else
    {
        if (ans<a) 
        {
            int x=(int)ceil((double)(a-ans)/(double)M);
            ans+=x*M;
        }
        printf("%lld\n",(b-ans)/M+1);
        printf("%lld",ans);
    }
}