CF1182F Maximum Sine【类欧，扩欧】

题目链接：洛谷

题目描述：求整数$x\in [a,b]$使得$|2px \ mod \ 2q-q|$最小，如果有多个$x$输出最小的。

数据范围：$1\leq a,b,p,q\leq 10^9$

---

第一道类欧的不是模板的题？？

首先一看就尝试一下二分，如何判断$2px \ mod \ 2q \in [l,r]$呢？我们发现

$$[2px \ mod \ 2q \in [l,r]]=\lfloor\frac{2px-l}{2q}\rfloor-\lfloor\frac{2px-r-1}{2q}\rfloor$$

所以存在$x\in [a,b]$使得$2px \ mod \ 2q \in [l,r]$，只需要判断$\sum_{x=a}^b(\lfloor\frac{2px-l}{2q}\rfloor-\lfloor\frac{2px-r-1}{2q}\rfloor)$是否大于0就可以了，这个使用类欧计算。

找到最小的偏差值$l$之后可以使用扩欧解同余方程。时间复杂度$O(\log^2 n)$。

 #include<bits/stdc++.h>
 #define Rint register int
 using namespace std;
 typedef long long LL;
 LL t, a, b, p, q, len;
 inline LL calc(LL a, LL b, LL c, LL n){
     if(!a || !n) return (n + ) * (b / c);
     if(n < ) return ;
     LL m = (a * n + b) / c;
     if(a >= c || b >= c) return calc(a % c, b % c, c, n) + (n * (n + ) >> ) * (a / c) + (n + ) * (b / c);
     return m * n - calc(c, c - b - , a, m - );
 }
 inline LL exgcd(LL a, LL b, LL &x, LL &y){
     if(!b){x = ; y = ; return a;}
     LL gcd = exgcd(b, a % b, y, x);
     y -= a / b * x;
     return gcd;
 }
 int main(){
     scanf("%lld", &t);
     while(t --){
         scanf("%lld%lld%lld%lld", &a, &b, &p, &q);
         LL l = , r = q, mid, P = p << , Q = q << ;
         while(l < r){
             mid = l + r >> ;
             LL L = q - mid, R = q + mid + ;
             if(calc(P, Q - L, Q, b) + calc(P, Q - R, Q, a - ) - calc(P, Q - L, Q, a - ) - calc(P, Q - R, Q, b) > ) r = mid;
             else l = mid + ;
         }
         LL x, y, gcd = exgcd(P, Q, x, y), ans = 1e9; P /= gcd; Q /= gcd;
         if((q - l) % gcd == ){
             LL xx = (q - l) / gcd * x; xx += (a - xx) / Q * Q;
             while(xx >= a) xx -= Q; while(xx < a) xx += Q;
             ans = min(ans, xx);
         }
         if((q + l) % gcd == ){
             LL xx = (q + l) / gcd * x; xx += (a - xx) / Q * Q;
             while(xx >= a) xx -= Q; while(xx < a) xx += Q;
             ans = min(ans, xx);
         }
         printf("%lld\n", ans);
     }
 }

CF1182F

后来看了一波官方题解，结果发现是一个类似BSGS的分块，时间复杂度$O(\sqrt{n})$。（写得不好还要带一个log...）

（还是上面那个方法更好一些）

（最主要是看起来更高大上一点。。。）