結果

問題 No.2428 Returning Shuffle
ユーザー tnakao0123
提出日時 2023-08-22 14:50:34
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 250 ms / 2,000 ms
コード長 3,493 bytes
コンパイル時間 565 ms
コンパイル使用メモリ 59,092 KB
実行使用メモリ 17,060 KB
最終ジャッジ日時 2024-12-17 14:01:20
合計ジャッジ時間 4,452 ms
ジャッジサーバーID
(参考情報)
judge4 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

/* -*- coding: utf-8 -*-
*
* 2428.cc: No.2428 Returning Shuffle - yukicoder
*/
#include<cstdio>
#include<vector>
#include<algorithm>
#include<utility>
using namespace std;
/* constant */
const int MAX_N = 1000000;
const int MAX_P = 1000000;
const int MOD = 998244353;
/* typedef */
typedef vector<int> vi;
typedef pair<int,int> pii;
typedef vector<pii> vpii;
template<const int MOD>
struct MI {
int v;
MI(): v() {}
MI(int _v): v(_v % MOD) {}
MI(long long _v): v(_v % MOD) {}
MI operator+(const MI m) const { return MI(v + m.v); }
MI operator-(const MI m) const { return MI(v + MOD - m.v); }
MI operator*(const MI m) const { return MI((long long)v * m.v); }
MI &operator+=(const MI m) { return (*this = *this + m); }
MI &operator-=(const MI m) { return (*this = *this - m); }
MI &operator*=(const MI m) { return (*this = *this * m); }
bool operator==(const MI m) const { return v == m.v; }
bool operator!=(const MI m) const { return v != m.v; }
MI pow(int n) const { // a^n % MOD
MI pm = 1, a = *this;
while (n > 0) {
if (n & 1) pm *= a;
a *= a;
n >>= 1;
}
return pm;
}
MI inv() const { return pow(MOD - 2); }
MI operator/(const MI m) const { return *this * m.inv(); }
MI &operator/=(const MI m) { return (*this = *this / m); }
};
typedef MI<MOD> mi;
/* global variables */
bool primes[MAX_P + 1];
int ps[MAX_N], ss[MAX_N], qs[MAX_N];
bool used[MAX_N];
/* subroutines */
int gen_primes(int maxp, vi &pnums) {
fill(primes, primes + maxp + 1, true);
primes[0] = primes[1] = false;
int p;
for (p = 2; p * p <= maxp; p++)
if (primes[p]) {
pnums.push_back(p);
for (int q = p * p; q <= maxp; q += p) primes[q] = false;
}
for (; p <= maxp; p++)
if (primes[p]) pnums.push_back(p);
return (int)pnums.size();
}
bool prime_decomp(int n, vi &pnums, vpii& pds) {
pds.clear();
int pn = pnums.size();
for (int i = 0; i < pn; i++) {
int pi = pnums[i];
if (pi * pi > n) {
if (n > 1) pds.push_back(pii(n, 1));
return true;
}
if (n % pi == 0) {
int fi = 0;
while (n % pi == 0) n /= pi, fi++;
pds.push_back(pii(pi, fi));
}
}
return false;
}
/* main */
int main() {
vi pnums;
gen_primes(MAX_P, pnums);
int n, m;
scanf("%d%d", &n, &m);
for (int i = 0; i < n; i++) ps[i] = i;
for (int i = 0; i < m; i++) {
int ti;
scanf("%d", &ti);
for (int j = 0; j < ti; j++) scanf("%d", ss + j), ss[j]--;
for (int j = 0; j < ti; j++) qs[j] = ps[ss[j]];
for (int j = 0; j < ti; j++) ps[ss[(j + 1) % ti]] = qs[j];
}
//for (int i = 0; i < n; i++) printf("%d ", ps[i]); putchar('\n');
vpii lpds;
for (int st = 0; st < n; st++)
if (! used[st]) {
int l = 0;
for (int u = st; ! used[u]; u = ps[u])
used[u] = true, l++;
if (l > 1) {
vpii pds, lpds1;
prime_decomp(l, pnums, pds);
int i = 0, j = 0, n0 = lpds.size(), n1 = pds.size();
while (i < n0 && j < n1) {
if (lpds[i].first < pds[j].first)
lpds1.push_back(lpds[i++]);
else if (lpds[i].first > pds[j].first)
lpds1.push_back(pds[j++]);
else {
int e = max(lpds[i].second, pds[j].second);
lpds1.push_back(pii(lpds[i].first, e));
i++, j++;
}
}
while (i < n0)
lpds1.push_back(lpds[i++]);
while (j < n1)
lpds1.push_back(pds[j++]);
swap(lpds, lpds1);
}
}
mi p = 1;
for (auto &pd: lpds)
p *= mi(pd.first).pow(pd.second);
printf("%d\n", p.v);
return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0