結果
問題 | No.2262 Fractions |
ユーザー | 👑 potato167 |
提出日時 | 2023-04-07 23:29:41 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 1,533 ms / 2,000 ms |
コード長 | 3,447 bytes |
コンパイル時間 | 2,185 ms |
コンパイル使用メモリ | 216,908 KB |
実行使用メモリ | 44,192 KB |
最終ジャッジ日時 | 2024-10-02 20:39:38 |
合計ジャッジ時間 | 39,797 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 820 ms
30,112 KB |
testcase_01 | AC | 103 ms
5,376 KB |
testcase_02 | AC | 107 ms
5,376 KB |
testcase_03 | AC | 113 ms
5,376 KB |
testcase_04 | AC | 132 ms
5,376 KB |
testcase_05 | AC | 125 ms
5,376 KB |
testcase_06 | AC | 135 ms
5,376 KB |
testcase_07 | AC | 135 ms
5,376 KB |
testcase_08 | AC | 128 ms
5,376 KB |
testcase_09 | AC | 138 ms
5,376 KB |
testcase_10 | AC | 136 ms
5,376 KB |
testcase_11 | AC | 534 ms
5,376 KB |
testcase_12 | AC | 525 ms
5,376 KB |
testcase_13 | AC | 538 ms
5,376 KB |
testcase_14 | AC | 544 ms
5,376 KB |
testcase_15 | AC | 531 ms
5,376 KB |
testcase_16 | AC | 95 ms
5,376 KB |
testcase_17 | AC | 94 ms
5,376 KB |
testcase_18 | AC | 95 ms
5,376 KB |
testcase_19 | AC | 1,112 ms
28,564 KB |
testcase_20 | AC | 1,030 ms
26,672 KB |
testcase_21 | AC | 1,148 ms
29,632 KB |
testcase_22 | AC | 1,022 ms
26,236 KB |
testcase_23 | AC | 848 ms
21,604 KB |
testcase_24 | AC | 1,419 ms
43,308 KB |
testcase_25 | AC | 1,447 ms
44,056 KB |
testcase_26 | AC | 1,408 ms
43,072 KB |
testcase_27 | AC | 1,423 ms
43,372 KB |
testcase_28 | AC | 1,408 ms
42,880 KB |
testcase_29 | AC | 1,476 ms
43,932 KB |
testcase_30 | AC | 1,459 ms
44,188 KB |
testcase_31 | AC | 1,533 ms
44,060 KB |
testcase_32 | AC | 1,456 ms
44,064 KB |
testcase_33 | AC | 1,463 ms
44,064 KB |
testcase_34 | AC | 1,452 ms
44,064 KB |
testcase_35 | AC | 307 ms
43,936 KB |
testcase_36 | AC | 300 ms
43,976 KB |
testcase_37 | AC | 286 ms
44,024 KB |
testcase_38 | AC | 299 ms
43,904 KB |
testcase_39 | AC | 1,094 ms
23,804 KB |
testcase_40 | AC | 1,097 ms
23,812 KB |
testcase_41 | AC | 1,101 ms
23,680 KB |
testcase_42 | AC | 1,112 ms
23,804 KB |
testcase_43 | AC | 1,093 ms
23,692 KB |
testcase_44 | AC | 1,462 ms
44,064 KB |
testcase_45 | AC | 1,483 ms
44,192 KB |
ソースコード
#include <bits/stdc++.h> #pragma GCC optimize("Ofast") #pragma GCC optimize("unroll-loops") using namespace std; using std::cout; using std::cin; using std::endl; using ll=long long; using ld=long double; ll ILL=2167167167167167167; const int INF=2100000000; const int mod=998244353; #define rep(i,a,b) for (ll i=a;i<b;i++) #define all(p) p.begin(),p.end() template<class T> using _pq = priority_queue<T, vector<T>, greater<T>>; template<class T> ll LB(vector<T> &v,T a){return lower_bound(v.begin(),v.end(),a)-v.begin();} template<class T> ll UB(vector<T> &v,T a){return upper_bound(v.begin(),v.end(),a)-v.begin();} template<class T> bool chmin(T &a,const T &b){if(a>b){a=b;return 1;}else return 0;} template<class T> bool chmax(T &a,const T &b){if(a<b){a=b;return 1;}else return 0;} template<class T> void So(vector<T> &v) {sort(v.begin(),v.end());} template<class T> void Sore(vector<T> &v) {sort(v.begin(),v.end(),[](T x,T y){return x>y;});} void yneos(bool a){if(a) cout<<"Yes\n"; else cout<<"No\n";} template<class T> void vec_out(vector<T> &p){for(int i=0;i<(int)(p.size());i++){if(i) cout<<" ";cout<<p[i];}cout<<"\n";} template<class T> T vec_min(vector<T> &a){assert(!a.empty());T ans=a[0];for(auto &x:a) chmin(ans,x);return ans;} template<class T> T vec_max(vector<T> &a){assert(!a.empty());T ans=a[0];for(auto &x:a) chmax(ans,x);return ans;} template<class T> T vec_sum(vector<T> &a){assert(!a.empty());T ans=a[0]-a[0];for(auto &x:a) ans+=x;return ans;} //Mobius(n)=0 exists p n|p*p // =-1 if n=p_1*p_2*...*p_{2a+1} // =1 if n=p_1*p_2*...*p_{2a} std::vector<int> Mobius(int N){ std::vector<int> p(N+1,-2); p[1]=1; for(int i=2;i<=N;i++){ if(p[i]!=-2) continue; p[i]=-1; for(int j=i*2;j<=N;j+=i){ if(p[j]==-2) p[j]=-1; else p[j]*=-1; } if(N/i<i) continue; for(int j=i*i;j<=N;j+=i*i) p[j]=0; } for(int i=1;i<=N;i++){ if(p[i]==-2) p[i]=-1; } return p; } void solve(); // oddloop int main() { ios::sync_with_stdio(false); cin.tie(nullptr); int t=1; cin>>t; rep(i,0,t) solve(); } void solve(){ ll N,K; cin>>N>>K; vector<vector<int>> table(N+1); auto M=Mobius(N); rep(i,1,N+1){ for(int j=i;j<=N;j+=i) table[j].push_back(i); } int C=0; vector<ll> q(N+1,-1); vector<ll> r(N+1); auto calc=[&](ll A,ll B)->ll{ C++; vector<ll> p(1+N); rep(i,1,N+1) p[i]=min(N,(A*i)/B); ll res=0; rep(i,1,N+1){ if(p[i]!=q[i]){ r[i]=0; for(auto x:table[i]) r[i]+=M[x]*(p[i]/x); } res+=r[i]; } swap(p,q); return res; }; if(calc(N,1)<K){ cout<<"-1\n"; return; } bool rev=0; auto out=[&](ll a,ll b)->void{ //cout<<C<<endl; if(rev) swap(a,b); cout<<a<<"/"<<b<<"\n"; }; if(calc(N,1)==K){ out(N,1); return; } if(K==1){ out(1,N); return; } ll h=calc(1,1); if(h==K){ out(1,1); return; } if(h<K){ K=2*h-K; rev=1; } ll L=N,R=N*N; while(R-L>1){ ll med=(R+L)/2; auto tmp=calc(med,N*N); if(tmp<K) L=med; else R=med; } rep(i,1,N+1){ ll ans=(R*i)/(N*N); if(L*i<N*N*ans&&ans<N){ out(ans,i); return; } } }