bzoj2187: fraction（类欧几里得）

传送门
题意简述：多组询问，每次给出$a,b,c,d$，求满足$\frac{a}{b}<\frac{p}{q}<\frac{c}{d}$的所有二元组$(p,q)$中$q$为第一关键字，$p$为第二关键字排出来字典序最小的那一对。

---

思路：
设计函数$f(a,b,p,q,c,d)$。
按照题目中保证$q$最小的要求考虑该函数几个边界：

1. $\lfloor \frac{a}{b}\rfloor +1\le \lceil \frac{c}{d}\rceil -1$，这个时候$p=\lfloor \frac{a}{b}\rfloor +1,q=1$时字典序最小。
2. $a=0$，这个时候$0<\frac{p}{q}<\frac{c}{d}\Rightarrow q>\frac{dp}{c}$，显然$p=1,q=\lfloor \frac{c}{d}\rfloor +1$的时候字典序最小。

然后考虑如何辗转缩小问题规模：

1. $a>b$ $or$ $c>d:原式\iff \frac{a\%b}{b}<\frac{p}{q}-\lfloor \frac{a}{b}\rfloor <\frac{c}{d}-\lfloor \frac{a}{b}\rfloor$
   $\Rightarrow f(a,b,p,q,c,d)=f(a\%b,b,p,q,c-\lfloor \frac{a}{b}\rfloor d,d),p+=\lfloor \frac{a}{b}\rfloor q$
2. $a\le b$ $and$ $c\le d:原式\iff \frac{d}{c}<\frac{q}{p}<\frac{b}{a}$
   $\Rightarrow f(a,b,p,q,c,d)=f(d,c,q,p,b,a)$这样会回到第一步。

代码：

#include<bits/stdc++.h>
#define ri register int
using namespace std;
typedef long long ll;
ll a,b,c,d,p,q;
inline ll gcd(ll a,ll b){while(b){ll t=a;a=b,b=t-t/a*a;}return a;}
inline void calc(ll&a,ll&b){ll d=gcd(a,b);a/=d,b/=d;}
inline void solve(ll a,ll b,ll&p,ll&q,ll c,ll d){
	calc(a,b),calc(c,d);
	if(!a){p=1,q=d/c+1;return;}
	ll x=a/b+1,y=c/d+(c%d>0)-1;
	if(x<=y){p=x,q=1;return;}
	if(a<=b&&c<=d)return solve(d,c,q,p,b,a);
	solve(a-a/b*b,b,p,q,c-d*(a/b),d),p+=q*(a/b);
}
int main(){
	while(~scanf("%lld%lld%lld%lld",&a,&b,&c,&d))solve(a,b,p,q,c,d),printf("%lld/%lld\n",p,q);
	return 0;
}