結果
問題 | No.2066 Simple Math ! |
ユーザー | gs15120 |
提出日時 | 2022-09-14 12:40:43 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 396 ms / 2,000 ms |
コード長 | 1,783 bytes |
コンパイル時間 | 2,755 ms |
コンパイル使用メモリ | 186,988 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-12-14 23:54:08 |
合計ジャッジ時間 | 12,085 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 31 |
ソースコード
#include<bits/stdc++.h> using namespace std; #define ll long long #define lll long long//__int128 #define pii pair<int,int> //#define f first #define s second # define pi 3.14159265358979323846 #include <ext/pb_ds/assoc_container.hpp> using namespace __gnu_pbds; struct chash_key { ll operator()(pii x) const { return (ll)x.first*1000000000 + x.second; } }; //gp_hash_table<pii,int,chash_key> hs[2],p[2]; //gp_hash_table<ll,int> dp; /* ios_base::sync_with_stdio(false); cin.tie(0),cout.tie(0); cout << fixed; cout.precision(3); cout << x << '\n'; */ struct P { int x,y,z; bool operator <(const P &a)const { return y>a.y; }; }; int a,b,c; const ll mod=998244353; //sum (ai+b)/c for i=0 to n long long FloorSumAP(long long a, long long b, long long c, long long n){ //cout<<a<<b<<c<<n<<"\n"; if(!a) return (b / c) * (n + 1); if(a >= c or b >= c) return ( ( n * (n + 1) ) / 2) * (a / c) + (n + 1) * (b / c) + FloorSumAP(a % c, b % c, c, n); long long m = (a * n + b) / c; return m * n - FloorSumAP(c, c - b - 1, a, m - 1); } int main() { ios_base::sync_with_stdio(false); cin.tie(0),cout.tie(0); int ttt=100; cin>>ttt; a=3,b=4; for(int tt=1;tt<=ttt;tt++) { c=tt; cin>>a>>b>>c; ll x=1,y=(ll)a*c; int g=__gcd(a,b); ll L=(ll)a*b/g; for(;x<y;) { ll h=(x+y)/2; ll d=0; if(h<L) { d=FloorSumAP(a,h%a+b,b,h/a)-1; } else { d=FloorSumAP(a,h-L+a+b,b,L/a-1)-1; } if(d<c) x=h+1; else y=h; } cout<<x<<"\n"; } //cout<<FloorSumAP(3,4%3+4,4,4/3)<<" "<<FloorSumAP(3,5%3+4,4,5/3); }