扩展欧几里得算法，中国剩余定理（Two Arithmetic Progressions，cf 710D）

理论知识

http://blog.csdn.net/xiaofengsheng/article/details/4813170

参考代码

http://blog.csdn.net/miracle_ma/article/details/52290256

http://blog.csdn.net/wust_zzwh/article/details/52287769

WA在了一组大质数，原因是解法不够优秀导致爆long long。

一度想用高精度解决，后来取了一下余就OK了。

感觉自己自从会了高精度，啥都想直接上高精度。

事实上高精度是绝对不需要滥用的，一个非常普遍而又简单的方法就是取余。

至于在哪里可以取余，就看自己去发现了。

本题中x是等间距的，因此只要找到任意一个x，以及间距，就能找出区间内所有的x。

如果只是求区间内x的个数的话，只需要对x稍加变换，然后带入公式就可以求得。

间距即是a1,a2的最小公倍数。

x可用中国剩余定理+扩展欧几里得算法求得。

由于我求x的方法不够好，所以爆了long long，但取余后得以解决。

x=a1*k+b1   ==>   x=b1%a1

x=a2*l+b2   ==>   x=b2%a2

一号方程组          二号方程组

由一号方程组得：

a1*k+b1=a2*l+b2   ==>   b1-b2=a2*l-a1*k

由中国剩余定理得：

方程组有一个小于[a1,a2]（a1，a2的最小公倍数）的非负整数解的充分必要条件是（b1-b2）%（a1，a2）==0

令bei=（b1-b2）/（a1，a2）

我们可以通过扩展欧几里得算法求得k与l使（a1,a2）=a2*l-a1*k

再将k和l扩大bei倍，就可以得到新的k与l。使得   b1-b2=a2*l-a1*k

此时的k与l就是一号方程组的某一个解。

带入一号方程组中任意一个方程。

x=a2*l+b2

就可以求出x。

红色的地方就是可能爆的地方。

bei是有可能爆int的，原来的l也可以十分接近int_max。

那么扩大后的l就可能爆long long了。

就算没爆，稍微乘个小东西也爆了。比如a2。

事实上我们只用求出任意一个x即可。

那么新的l显然没有必要那么大。

x是以   lcm(a1,a2)   为周期的。

容易得出在   x=a2*l+b2   中的l，也就是新的l，则是以   lcm(a1,a2)/a2   为周期的。

因此只要在   l=l*bei   后面加一句   l%=a1   这样求得的l就不会超过int_max了。

怕中途爆了可以这样 l=((l%a1)*(bei%a1))%a1。

代码

#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
ll a1,b1,a2,b2,L,R;

ll gcd(ll a,ll b)
{
    return a%b==0?b:gcd(b,a%b);
}

ll lcm(ll a,ll b)
{
    return a/gcd(a,b)*b;
}

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()
{
    scanf("%I64d %I64d %I64d %I64d %I64d %I64d",&a1,&b1,&a2,&b2,&L,&R);
    if(a1<a2)
    {
        swap(a1,a2);
        swap(b1,b2);
    }
    ll gcda1a2=gcd(a1,a2);
    ll lcma1a2=lcm(a1,a2);
    if((b1-b2)%gcda1a2==0)
    {
        ll bei=(b1-b2)/gcda1a2;
        ll k,l;
        exgcd(a1,a2,k,l);
        k=-k;
        k=((k%a2)*(bei%a2))%a2;
        l=((l%a1)*(bei%a1))%a1;
        ll x=a2*l+b2;
        //ll x=a1*k+b1;
        L=max(L,max(b1,b2));
        if(x>=L) x=x-lcma1a2*((x-L)/lcma1a2+1);
        ll ans=(R-x)/lcma1a2-(L-x-1)/lcma1a2;
        if(ans<=0) puts("0");
        else printf("%I64d\n",ans);
    }
    else puts("0");
    return 0;
}