結果
問題 | No.2262 Fractions |
ユーザー | 👑 potato167 |
提出日時 | 2023-04-07 22:59:21 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
TLE
|
実行時間 | - |
コード長 | 3,554 bytes |
コンパイル時間 | 2,423 ms |
コンパイル使用メモリ | 215,952 KB |
実行使用メモリ | 40,612 KB |
最終ジャッジ日時 | 2024-10-02 20:18:51 |
合計ジャッジ時間 | 8,911 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | TLE | - |
testcase_01 | RE | - |
testcase_02 | RE | - |
testcase_03 | RE | - |
testcase_04 | RE | - |
testcase_05 | RE | - |
testcase_06 | RE | - |
testcase_07 | RE | - |
testcase_08 | RE | - |
testcase_09 | RE | - |
testcase_10 | RE | - |
testcase_11 | RE | - |
testcase_12 | TLE | - |
testcase_13 | RE | - |
testcase_14 | RE | - |
testcase_15 | RE | - |
testcase_16 | AC | 123 ms
6,820 KB |
testcase_17 | AC | 124 ms
6,816 KB |
testcase_18 | AC | 121 ms
6,820 KB |
testcase_19 | TLE | - |
testcase_20 | -- | - |
testcase_21 | -- | - |
testcase_22 | -- | - |
testcase_23 | -- | - |
testcase_24 | -- | - |
testcase_25 | -- | - |
testcase_26 | -- | - |
testcase_27 | -- | - |
testcase_28 | -- | - |
testcase_29 | -- | - |
testcase_30 | -- | - |
testcase_31 | -- | - |
testcase_32 | -- | - |
testcase_33 | -- | - |
testcase_34 | -- | - |
testcase_35 | -- | - |
testcase_36 | -- | - |
testcase_37 | -- | - |
testcase_38 | -- | - |
testcase_39 | -- | - |
testcase_40 | -- | - |
testcase_41 | -- | - |
testcase_42 | -- | - |
testcase_43 | -- | - |
testcase_44 | -- | - |
testcase_45 | -- | - |
ソースコード
#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; 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]==0) break; for(auto x:table[i]) res+=M[x]*(p[i]/x); } return res; }; if(calc(N,1)<K){ cout<<"-1\n"; return; } auto out=[&](ll a,ll b)->void{ //cout<<C<<endl; cout<<a<<"/"<<b<<"\n"; }; if(calc(N,1)==K){ out(N,1); return; } if(K==1){ out(1,N); return; } pair<ll,ll> L={1,N},R={N,1}; while(true){ pair<ll,ll> med={-1,-1}; rep(i,1,N+1){ ll X=(L.first*i+L.second)/L.second; ll n_X=max(1ll,X); ll Y=(i*(R.first)-1)/R.second; ll n_Y=min(N,Y); if(n_X<=n_Y){ med={sqrt(n_X*n_Y),i}; break; } } if(C==2) med={1,1}; //cout<<L.first<<" "<<L.second<<" "<<R.first<<" "<<R.second<<endl; assert(med.first!=-1); auto tmp=calc(med.first,med.second); if(tmp==K){ out(med.first,med.second); return; } if(tmp<K) L=med; else R=med; } }