The Erdös-Straus Conjecture 暴力

题目链接：http://acm.csu.edu.cn:20080/csuoj/problemset/problem?pid=2276

题目大意：找到按字典序找到最小的三元组 $(a,b,c)$使其满足$\frac{1}{a}+\frac{1}{b}+\frac{1}{c}=\frac{4}{n}$

比赛的时候想了很多，没做出来，没想到赛后有人说是暴力。
开始以为是贪心，a肯定是$\lfloor \frac{n}{4}\rfloor +1$,然后再怎么弄，但是后续并没有什么进展。

下面给粗的做法似乎有点过于玄学。首先，枚举a的范围$[\lfloor \frac{n}{4}\rfloor +1,\lfloor \frac{n}{4}\rfloor +10]$,得到$\frac{4}{n}-\frac{1}{a}=\frac{x}{y}$,接着枚举b的范围$[\lfloor \frac{y}{x}\rfloor +1,\lfloor \frac{2\ast y}{x}\rfloor ]$，得到c验证即可。

a的范围大概率是猜出来的，由于我们假设b<=c,因此b的范围就是上面?那个。那么问题来了，x不会超过40，y大概是1e9的级别，就算x,y除了gcd那也得是1e7级别的，b在这么大的范围内枚举，而且还有1000个case，居然都没有T，不得不让我相信这是玄学(或者说这里面有什么理论依据？)…

完整代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll gcd(ll a,ll b){
	return b==0?a:gcd(b,a%b);
}
int main(){
	ll t,n;
	scanf("%lld",&t);
	for(int i=1;i<=t;i++){
		scanf("%d%lld",&i,&n);
		ll a,b,c,ok=0;
		for(a=n/4+1;a<=n/4+10;a++){
			ll x=4*a-n,y=n*a;
			ll g=gcd(x,y);
			x/=g,y/=g;
			for(b=y/x+1;b<=2*y/x;b++){
				if((y*b)%(x*b-y)==0){
					c=(y*b)/(x*b-y);
					ok=1;
					break;
				}
			}
			if(ok)
				break;
		}
		printf("%d %lld %lld %lld\n",i,a,b,c);
	}
	return 0;
}

另外吐槽一波__gcd在本地居然不能用…