https://codeforces.com/problemset/problem/773/A
一开始二分枚举d,使得(x+d)/(y+d)>=p/q&&x/(y+d)<=p/q,错在这些数是离散的,不能由两边异号判定一定存在这个交点。
然后改成枚举d,使得y=d*q,这样就一定是倍数了。然后就是要想清楚了,找不到这样卡在中间的d,其实都是因为d不够大的原因,d够大保证是可以的除非正确率是100%。
然后就是二分的上界,按道理q的最大值是1e9,y的最大值也是1e9,他们的公倍数肯定在1e18范围内(p和q任意组合能得到的比值最多就是1e18/2,把1e18都枚举完肯定能遍历所有能构造出的情况),那么最大值肯定是1e18/q,多个1都不行。
二分,老朋友了,while(1),当l==m的时候是边界,特判就好。
#includeusing namespace std;#define ll long longstring a,b;string c;int main(){ int t; scanf("%d",&t); while(t--){ ll x,y,p,q; scanf("%lld%lld%lld%lld",&x,&y,&p,&q); ll l=1,r=1e18/q,m; ll ans; while(1){ m=(l+r)>>1; //cout< <<" " < <<" "< < =0&&(x+del)>=l*p&&x<=l*p){ ans=l; } else if((r*q-y)>=0&&(x+(r*q-y))>=r*p&&x<=r*p){ ans=r; } else{ ans=-1; } break; } if(m*q =m*p&&x<=m*p){ r=m; //cout< <<" ok"<
其实还有直接用公式解的方法,设最终的状态为pt/qt,只要保证增量为正数即可,当然是pt>=x,qt>=y,还有错的题不会再变对,qt-pt>=y-x,三个式子直接解出来。
这个方法要注意特判p==q的情况,这种时候不能作除法,还有p==0和q==0的时候。