#include #if __has_include() #endif using namespace std; #define LL(...) ll __VA_ARGS__;lin(__VA_ARGS__) #define FO(n) for(ll IJK=n;IJK-->0;) #define defpp templatevoid pp(const auto&...a){[[maybe_unused]]const char*c="";((o<(a...);} #define entry defpp void main();void main2();}int main(){my::io();my::main();}namespace my{ #define multiple_testcases LL(T);FO(T)main2();}void main2(){ namespace my{ void io(){cin.tie(nullptr)->sync_with_stdio(0);cout<>...>>a);} constexpr auto abs(auto x){return x<0?-x:x;} constexpr ll size10(auto x){x|=1;ll r=0;while(x>0)x/=10,++r;return r;} templateconstexpr ll maxsize10(){return size10(numeric_limits::max());} auto floor(auto x,auto y){if(y<0)x=-x,y=-y;return x>=0?x/y:(x+1)/y-1;} bool amin(auto&a,const auto&b){return b>{a...});} auto min(auto...a){return min(initializer_list>{a...});} templatecommon_type_tgcd(T a,U b){return b?gcd(b,a%b):abs(a);} auto binary_search(ll o,ll x,const auto&f){while(abs(o-x)>1)(f((o+x)/2)?o:x)=(o+x)/2;return o;} templatestruct infinity{ templatestatic constexpr T ones(size_t n){return n?ones(n-1)*10+1:0;} templateconstexpr operator T()const{return ones(maxsize10())*(1-is_negative*2);} }; constexpr infinity oo; } namespace my{ templateT floor_sum(lll n,lll m,lll a,lll b){ if(n<=0)return 0; if(a<0)return-floor_sum(n,m,-a,-b+m-1); if(b<0)return floor_sum(n,m,a,b*(1-m))+(T)b*n; if(a>=m||b>=m)return floor_sum(n,m,a%m,b%m)+(T)(n-1)*n/2*(a/m)+(T)n*(b/m); if(a==0)return 0; return floor_sum((a*n+b)/m,a,m,(a*n+b)%m); } } namespace my{ templateT impl_ax_by_le_c_count(lll a,lll b,lll c,lll A,lll B,lll C,lll D){ if(A>B||C>D)return 0; if(a<0)return impl_ax_by_le_c_count(-a,b,c,-B,-A,C,D); if(b<0)return impl_ax_by_le_c_count(a,-b,c,A,B,-D,-C); assert(a>=0); assert(b>=0); if(a>0&&b==0)return impl_ax_by_le_c_count(b,a,c,C,D,A,B); if(a==0){ if(b==0)return c>=0?(B-A+1)*(D-C+1):0; else return (B-A+1)*max(min(floor(c,b),D)-C+1,0); } assert(a>0); assert(b>0); lll X=B-A,Y=D-C; c-=a*A+b*C; if(c<0)return 0; lll X2=min(X,floor(c,a)); if(c-b*Y<0){ return floor_sum(X2+1,b,-a,c)+(X2+1); }else{ lll Y2=floor(c-b*Y,a); amin(Y2,X2); return(Y2+1)*(Y+1)+(floor_sum(X2+1,b,-a,c)-floor_sum(Y2+1,b,-a,c))+(X2-Y2); } } templateT ax_by_le_c_count(lll a,lll b,lll c,lll A,lll B,lll C,lll D){return impl_ax_by_le_c_count(a,b,c,A,B-1,C,D-1);} } namespace my{entry void main(){ multiple_testcases LL(a,b,K);K++; ll g=gcd(a,b); a/=g,b/=g; ll X=binary_search(oo,0,[&](ll X){ // X以下でas+bt(0<=s=0)表示できる非負整数がK個以上あるか return ax_by_le_c_count(a,b,X,0,b,0,1e12)>=K; }); pp(X*g); }}