codeforces gym102059 Fractions 简单数学

https://codeforces.com/gym/102059/problem/H

题目大意：给出$A、B、C、D$(范围在$[1,10^{12}]$)，$x$满足$A<=x<=B$，$y$满足$C<=y<=D$，求出满足在最简形式下$x+y<1000$的$x/y$的个数。

思路：暴力计算出满足$x+y<1000$的所有最简分数，然后遍历这个序列开始统计答案：计算$A、B$分别是分子的多少倍数，计算$C、D$分别是分母的多少倍数，然后对两段新区间取交集即可。

#include<bits/stdc++.h>
#define INF 0x3f3f3f3f
typedef long long ll;
using namespace std;

const int maxn=1e6+5;

struct node
{
	int x,y;
	node(){}
	node(int xx,int yy)
	{
		x=xx,y=yy;
	}
}a[maxn];

int len;
ll A,B,C,D;

int gcd(int x,int y)
{
	return y==0?x:gcd(y,x%y);
}

void init()
{
	for(int i=1;i<999;i++)
	{
		for(int j=1;j<999;j++)
		{
			if(i+j>=1000)
				break;
			if(gcd(i,j)==1)
				a[len++]=node(j,i);
		}
	}
}

int main()
{
	init();
	scanf("%I64d %I64d %I64d %I64d",&A,&B,&C,&D);
	ll l1,r1,l2,r2,ans=0;
	for(int i=0;i<len;i++)
	{
		l1=ceil(A*1.0/a[i].x);
		r1=B/a[i].x;
		l2=ceil(C*1.0/a[i].y);
		r2=D/a[i].y;
		if(r1<l2||r2<l1)
			continue;
		ans+=min(r1,r2)-max(l1,l2)+1;
	}
	printf("%I64d\n",ans);
	return 0;
}