結果
問題 |
No.2045 Two Reflections
|
ユーザー |
![]() |
提出日時 | 2023-08-03 21:13:28 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 7 ms / 2,000 ms |
コード長 | 1,786 bytes |
コンパイル時間 | 6,014 ms |
コンパイル使用メモリ | 319,612 KB |
実行使用メモリ | 6,824 KB |
最終ジャッジ日時 | 2024-10-13 18:21:09 |
合計ジャッジ時間 | 6,684 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 27 |
ソースコード
#include<bits/stdc++.h> using namespace std; #include <atcoder/all> using namespace atcoder; using mint=modint998244353; //1000000007; using ll=long long; using pp=pair<int,int>; #define sr string #define vc vector #define fi first #define se second #define rep(i,n) for(int i=0;i<(int)n;i++) #define pb push_back #define all(v) v.begin(),v.end() #define pque priority_queue #define bpc(a) __builtin_popcount(a) struct divs{ int n; vector<int> f, primes; divs(int n=1):n(n), f(n+1) { f[0] = f[1] = -1; for (ll i = 2; i <= n; ++i) { if (f[i]) continue; primes.push_back(i); f[i] = i; for (ll j = i*i; j <= n; j += i) { if (!f[j]) f[j] = i; } } } bool isPrime(int x) { return f[x] == x;} vector<int> factorList(int x) { vector<int> res; while (x != 1) { res.push_back(f[x]); x /= f[x]; } return res; } vector<pair<int,int>> factor(int x) { vector<int> fl = factorList(x); if (fl.size() == 0) return {}; vector<pair<int,int>> res(1, pair<int,int>(fl[0], 0)); for (int p : fl) { if (res.back().first == p) { res.back().second++; } else { res.emplace_back(p, 1); } } return res; } }; int main(){ int n,l,r;cin>>n>>l>>r; if(l==1&&r==1){ cout<<1; return 0; } if(l==1||r==1){ cout<<2; return 0; } vc<int>a(n); rep(i,n)a[i]=i; rep(i,l)a[i]=l-1-i; auto pre=a; rep(i,r)a[n-r+i]=pre[n-1-i]; dsu d(n); rep(i,n)d.merge(i,a[i]); map<int,int>m; divs di(n); rep(i,n)if(d.leader(i)==i){ int t=d.size(i); for(auto [x,y]:di.factor(t)){ m[x]=max(m[x],y); } } mint ans=2; for(auto [x,y]:m)ans*=mint(x).pow(y); //rep(i,n)cout<<a[i]<<' '; cout<<ans.val(); }