HDU 6624 fraction【数论】

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=6624
题目大意：给定$p,x$,存在两整数$a,b$，求解满足$a<b,a\equiv bx(modp)$的最小值$b$。

思路：
令$a=bx-cp$，则存在不等式$0\le bx-cp\le b$，
变换一下得到：$\frac{p}{x}<\frac{b}{c}<\frac{p}{x-1}$,我们要解出这个不等式中的$min(b)$值。
这时候可以想到运用类似于辗转相除法的形式来解决这个问题:
（1）令$z_1=\lfloor \frac{p}{x}\rfloor ,z_2=\lfloor \frac{p}{x-1}\rfloor$，
（2）如果$z_1!=z_2$，则令$c=1,b=z_1+1$(因为不含等号，所以要加1)跳出循环。否则就执行（3）。
（3）不等式同步减掉$z_1$,不等式变为：$\frac{p-z_1\ast x}{x}<\frac{b-z_1\ast c}{c}<\frac{p-z_1\ast (x-1)}{x-1}$，分子分母翻转一下得到：$\frac{x-1}{p-z_1\ast (x-1)}<\frac{c}{b-z_1\ast c}<\frac{x}{p-z_1\ast x}$，这个式子和最初得到的原式$(\frac{p}{x}<\frac{b}{c}<\frac{p}{x-1})$是不是有点相像，两式对应一下递归就$ok$了。

证明：算法正确性菜鸡无法证明，很像辗转相除法。
AC代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll p,xx;
ll b,c;
void solve(ll p,ll x,ll &b,ll &c,ll pp,ll y)//y=x-1
{
    ll z1=p/x;
    ll z2=pp/y;
    if(z1!=z2){
       c=1;
       b=z1+1;//因为不等式不带等号
       return ;
    }
    p-=z1*x;
    pp-=z1*y;
    b-=c*z1;
    solve(y,pp,c,b,x,p);
    b+=c*z1;
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--){
       scanf("%lld%lld",&p,&xx);
       b=0;
       c=0;
       solve(p,xx,b,c,p,xx-1);//按照不等式从左到右依次传参
       ll a=b*xx-c*p;
       printf("%lld/%lld\n",a,b);
    }
    return 0;
}
/*
3
11 7
998244353 554580197
998244353 998244352

2/5
8/9
499122176/499122177
*/