#include using namespace std; #if __has_include() #include using namespace atcoder; templateistream &operator>>(istream &is,static_modint &a){long long b;is>>b;a=b;return is;} istream &operator>>(istream &is,modint &a){long long b;cin>>b;a=b;return is;} #endif #ifdef LOCAL #include "debug.h" #else #define debug(...) static_cast(0) #define debugg(...) static_cast(0) templateostream &operator<<(ostream &os,const pair&p){os<; templateusing minque=priority_queue,greater>; templatebool chmax(T &a,const T &b){return (abool chmin(T &a,const T &b){return (a>b?(a=b,true):false);} templateistream &operator>>(istream &is,pair&p){is>>p.first>>p.second;return is;} templateistream &operator>>(istream &is,vector &a){for(auto &i:a)is>>i;return is;} templatevoid operator++(pair&a,int n){a.first++,a.second++;} templatevoid operator--(pair&a,int n){a.first--,a.second--;} templatevoid operator++(vector&a,int n){for(auto &i:a)i++;} templatevoid operator--(vector&a,int n){for(auto &i:a)i--;} #define reps(i,a,n) for(int i=(a);i<(n);i++) #define rep(i,n) reps(i,0,n) #define all(x) x.begin(),x.end() #define pcnt(x) __builtin_popcount(x) ll myceil(ll a,ll b){return (a+b-1)/b;} template auto vec(const int (&d)[n],const T &init=T()){ if constexpr (id(d,init)); else return init; } void SOLVE(); int main(){ ios::sync_with_stdio(false); cin.tie(nullptr); #ifdef LOCAL clock_t start=clock(); #endif int testcase=1; cin>>testcase; for(int i=0;i T floorsum(T n,T m,T a,T b){ assert(n>=0&&m>=1); T ret=0; if(a<0){ T aa=a%m+m; ret-=n*(n-1)/2*((aa-a)/m); a=aa; } if(b<0){ T bb=b%m+m; ret-=n*((bb-b)/m); b=bb; } while(true){ if(a>=m){ ret+=n*(n-1)/2*(a/m); a%=m; } if(b>=m){ ret+=n*(b/m); b%=m; } T ymx=a*n+b; if(ymx>p>>q>>k; ll g=gcd(p,q); p/=g,q/=g; ll ok=1e18,ng=0; while(ok-ng>1){ ll mid=(ok+ng)/2; ll n=min(mid/p+1,q); if(floorsum<__int128_t>(n,q,p,mid-p*(n-1))+n-1>=k)ok=mid; else ng=mid; } cout<