結果

問題 No.3398 Accuracy of Integer Division Approximate Function 2
コンテスト
ユーザー 👑 Mizar
提出日時 2025-12-04 21:28:55
言語 C++17
(gcc 13.3.0 + boost 1.89.0)
結果
AC  
実行時間 238 ms / 2,000 ms
コード長 508 bytes
記録
コンパイル時間 5,731 ms
コンパイル使用メモリ 426,936 KB
実行使用メモリ 7,720 KB
最終ジャッジ日時 2025-12-04 23:34:06
合計ジャッジ時間 8,174 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<bits/stdc++.h>
#include<boost/multiprecision/cpp_int.hpp>
using namespace std;using namespace boost::multiprecision;int main(){int t;cin>>t;while(t--){cpp_int D,A,B,K,n,a,b,c,d,m,x,g,r,l=0,h,q,s,v;cin>>D>>A>>B>>K;D/=g=gcd(D,A);A/=g;if(A*B%D<1&&D*K+2>A){cout<<-1<<'\n';}else{h=A*B*K+2;while(l+1<h){x=n=(l+h)/2;a=B;b=-A*B/D;c=D;m=A;s=v=d=0;while(a+=c/m*b,c%=m,s+=d/m*b,d%=m,v=max(v,max(s,s+a*~-n+b*(q=(c*~-n+d)/m))),n=q){s+=(a+b)*(a<0);d=m+~d;swap(a,b);swap(c,m);}(v<=B*K?l:h)=x;}cout<<g*D*l<<'\n';}}}
0