結果
問題 | No.1931 Fraction 2 |
ユーザー | ytqm3 |
提出日時 | 2022-04-27 22:26:52 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 401 ms / 2,000 ms |
コード長 | 1,792 bytes |
コンパイル時間 | 4,860 ms |
コンパイル使用メモリ | 267,972 KB |
実行使用メモリ | 29,288 KB |
最終ジャッジ日時 | 2023-09-20 22:29:50 |
合計ジャッジ時間 | 17,869 ms |
ジャッジサーバーID (参考情報) |
judge14 / judge13 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 11 ms
14,696 KB |
testcase_01 | AC | 12 ms
14,696 KB |
testcase_02 | AC | 11 ms
14,840 KB |
testcase_03 | AC | 13 ms
14,960 KB |
testcase_04 | AC | 13 ms
14,772 KB |
testcase_05 | AC | 13 ms
14,792 KB |
testcase_06 | AC | 13 ms
14,888 KB |
testcase_07 | AC | 14 ms
14,956 KB |
testcase_08 | AC | 12 ms
14,904 KB |
testcase_09 | AC | 13 ms
14,724 KB |
testcase_10 | AC | 12 ms
14,896 KB |
testcase_11 | AC | 14 ms
14,768 KB |
testcase_12 | AC | 13 ms
14,844 KB |
testcase_13 | AC | 81 ms
17,504 KB |
testcase_14 | AC | 255 ms
23,908 KB |
testcase_15 | AC | 122 ms
19,288 KB |
testcase_16 | AC | 32 ms
15,560 KB |
testcase_17 | AC | 38 ms
15,840 KB |
testcase_18 | AC | 352 ms
28,216 KB |
testcase_19 | AC | 347 ms
27,616 KB |
testcase_20 | AC | 356 ms
28,060 KB |
testcase_21 | AC | 364 ms
27,936 KB |
testcase_22 | AC | 364 ms
27,588 KB |
testcase_23 | AC | 367 ms
28,192 KB |
testcase_24 | AC | 349 ms
27,876 KB |
testcase_25 | AC | 357 ms
27,904 KB |
testcase_26 | AC | 359 ms
27,764 KB |
testcase_27 | AC | 352 ms
27,544 KB |
testcase_28 | AC | 389 ms
29,020 KB |
testcase_29 | AC | 391 ms
29,116 KB |
testcase_30 | AC | 390 ms
28,968 KB |
testcase_31 | AC | 394 ms
29,000 KB |
testcase_32 | AC | 394 ms
29,288 KB |
testcase_33 | AC | 395 ms
29,040 KB |
testcase_34 | AC | 396 ms
29,116 KB |
testcase_35 | AC | 397 ms
29,032 KB |
testcase_36 | AC | 393 ms
28,880 KB |
testcase_37 | AC | 401 ms
29,172 KB |
testcase_38 | AC | 13 ms
14,812 KB |
ソースコード
#include<bits/stdc++.h> #include<atcoder/all> using i64=long long; using u64=unsigned long long; using namespace std; template<u64 mod> using modint=atcoder::static_modint<mod>; struct primefact{ int n; vector<i64> lpf; primefact(int n_):n(n_),lpf(n_+1,-1){ for(i64 p=2;p<=n;++p){ if(lpf[p]!=-1){ continue; } lpf[p]=p; for(i64 q=p;q<=n;q+=p){ if(lpf[q]==-1){ lpf[q]=p; } } } } vector<pair<i64,i64>> operator()(int m){ vector<i64> v; while(m!=1){ v.emplace_back(lpf[m]); m/=lpf[m]; } if(v.size()==0){ return {}; } vector<pair<i64,i64>> res; res.emplace_back(v[0],1); for(size_t i=1;i<v.size();++i){ if(v[i-1]!=v[i]){ res.emplace_back(v[i],1); } else{ res.back().second++; } } return res; } }; int main(){ constexpr u64 mod=998244353; typedef modint<mod> mint; primefact pf(300'001); int N; cin>>N; vector<i64> A(N),B(N); vector<vector<pair<i64,i64>>> idx(300'001); vector<i64> m(300'001); mint sm=0; for(int i=0;i<N;++i){ cin>>A[i]>>B[i]; sm+=mint(A[i])/B[i]; auto tmp=pf(B[i]); for(auto [p,e]:tmp){ idx[p].emplace_back(i,e); m[p]=max(m[p],e); } } mint d=1; for(i64 p=2;p<=300'001;++p){ if(m[p]==0){ continue; } vector<i64> pw(21,1); for(int i=1;i<=20;++i){ pw[i]=pw[i-1]*p; } i64 sm=0; for(auto [i,e]:idx[p]){ sm+=A[i]*pw[m[p]-e]*atcoder::inv_mod(B[i]/pw[e],pw[m[p]]); sm%=pw[m[p]]; } i64 prd=pw[m[p]]; for(int i=0;i<m[p];++i){ if(sm%p==0){ sm/=p; prd/=p; } else{ break; } } d*=prd; } cout<<(sm*d).val()<<" "<<d.val()<<endl; }