結果
問題 | No.1931 Fraction 2 |
ユーザー | 👑 potato167 |
提出日時 | 2022-05-06 23:34:48 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
RE
|
実行時間 | - |
コード長 | 2,539 bytes |
コンパイル時間 | 2,361 ms |
コンパイル使用メモリ | 210,648 KB |
実行使用メモリ | 57,236 KB |
最終ジャッジ日時 | 2023-09-20 22:39:31 |
合計ジャッジ時間 | 19,659 ms |
ジャッジサーバーID (参考情報) |
judge14 / judge13 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 118 ms
52,952 KB |
testcase_01 | AC | 115 ms
53,344 KB |
testcase_02 | AC | 117 ms
53,024 KB |
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 | RE | - |
testcase_13 | RE | - |
testcase_14 | WA | - |
testcase_15 | RE | - |
testcase_16 | RE | - |
testcase_17 | RE | - |
testcase_18 | AC | 429 ms
56,984 KB |
testcase_19 | AC | 428 ms
56,924 KB |
testcase_20 | AC | 425 ms
57,052 KB |
testcase_21 | AC | 421 ms
57,008 KB |
testcase_22 | AC | 430 ms
57,236 KB |
testcase_23 | AC | 425 ms
56,996 KB |
testcase_24 | AC | 424 ms
57,000 KB |
testcase_25 | AC | 417 ms
57,188 KB |
testcase_26 | AC | 421 ms
57,100 KB |
testcase_27 | AC | 419 ms
56,892 KB |
testcase_28 | AC | 437 ms
56,992 KB |
testcase_29 | AC | 432 ms
57,092 KB |
testcase_30 | AC | 436 ms
57,028 KB |
testcase_31 | AC | 439 ms
57,128 KB |
testcase_32 | AC | 435 ms
57,040 KB |
testcase_33 | AC | 424 ms
56,928 KB |
testcase_34 | AC | 422 ms
57,036 KB |
testcase_35 | AC | 432 ms
57,028 KB |
testcase_36 | AC | 435 ms
57,044 KB |
testcase_37 | AC | 443 ms
57,028 KB |
testcase_38 | AC | 113 ms
53,092 KB |
ソースコード
#include <bits/stdc++.h> #pragma GCC optimize("Ofast") #define _GLIBCXX_DEBUG 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 ll mod=998244353; #define rep(i,a) for (ll i=0;i<a;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";} ll jyo(ll x,ll y,ll z){ ll H=y; //ここから ll a=1,b=(x%z+z)%z,c=1; while(H>0){ a*=2; if(H%a!=0){ H-=a/2; c*=b; c%=z; } b*=b; b%=z; } //ここまで return c; } void solve(); // oddloop int main() { ios::sync_with_stdio(false); cin.tie(nullptr); int t=1; //cin>>t; rep(i,t) solve(); } void solve(){ ll N; cin>>N; ll L=2e5+10; vector<vector<int>> table(L); rep(i,L){ if(i<2) continue; if(!table[i].empty()) continue; for(int j=i;j<L;j+=i) table[j].push_back(i); } vector<ll> A(N),B(N); vector<int> p(N); vector<vector<ll>> dp(L,vector<ll>(20)); ll ans=0; rep(i,N) cin>>A[i]>>B[i]; rep(i,N){ //cout<<N<<" "<<i<<endl; ans+=(A[i]*jyo(B[i],mod-2,mod))%mod; ans%=mod; for(auto x:table[B[i]]){ //cout<<x<<endl; int tmp=B[i],C=0,X=1; while(tmp%x==0) tmp/=x,C++,X*=x; //cout<<X<<endl; chmax(p[x],C); if(x==2){ if(X==2) dp[x][C]++; else if(X==4) dp[x][C]+=(A[i]*jyo(tmp,1,4))%4; else dp[x][C]+=(A[i]*jyo(tmp,X/4-1,X))%4; } else dp[x][C]+=(A[i]*jyo(tmp,X-X/x-1,X))%X; //cout<<"#<"<<endl; } //cout<<"###"<<endl; } ll D=1; rep(i,L){ if(table[i].size()!=1) continue; for(int k=19;k>=0;k--){ if(dp[i][k]%i!=0){ D=(D*jyo(i,k,mod))%mod; break; } if(k!=0){ dp[i][k-1]+=dp[i][k]/i; } } } cout<<(ans*D)%mod<<" "<<D%mod<<"\n"; }