CodeForces - 710D Two Arithmetic Progressions 扩欧同余方程

题意：给一个式子x=a1*k+b1=a2*l+b2.(k,l>=0) ,问有多少满足条件的x属于[L,R].a1,a2,b1,b2,L,R已知

思路：化为同余方程

x%a1=b1 （1

x%a2=b2 （2

化为a1*n1+a2*n2=b2-b1.用扩展欧几里得解出最小正整数解n0.

通解N=n0+k1*(a2/d).将通解代入之后得X=a1*N+b1=a1*n0+b1+k1*lcm(a1,a2) （k1>=0）为方程（1，（2的通解。

这样就将式子合并了。弱化版的X问题

题目链接：http://codeforces.com/problemset/problem/710/D

#include<cstdio>
#include<algorithm>
#define ll long long
using namespace std;
void exgcd(ll a,ll b,ll &d,ll &x,ll &y)
{
	if(!b)
	{
		d=a;
		x=1;y=0;
	}
	else
	{
		exgcd(b,a%b,d,y,x);
		y-=x*(a/b);
	}
}
int main()
{
	ll a1,b1,a2,b2,L,R;//a1*x+a2*y=b2-b1
//	freopen("in.txt","r",stdin);
	while(~scanf("%I64d%I64d%I64d%I64d%I64d%I64d",&a1,&b1,&a2,&b2,&L,&R))
	{
		ll d,x,y;
		exgcd(a1,a2,d,x,y);
		ll c=b2-b1;
		if(c%d)
		{
			printf("0\n");
			continue;
		}
		else
		{
			x*=(c/d);
			ll t=a2/d;
			x=(x%t+t)%t;
			L=max(L,max(b1,b2));
			ll b=x*a1+b1;
			ll lcm=a1*(a2/d);
			ll ans=0;
			if(L>R)
			{
				printf("0\n");
				continue;
			}
			if(b<=R)ans+=1+(R-b)/lcm;  //R-b是排除了b，也就是[b+1,R]宽度.+1是b本身也要被算进来的
			if(b<L)ans-=1+(L-1-b)/lcm;//减去b到L-1的步数 
			printf("%I64d\n",ans);
		}
	}
}