結果
| 問題 |
No.2045 Two Reflections
|
| コンテスト | |
| ユーザー |
👑 |
| 提出日時 | 2022-08-20 14:30:23 |
| 言語 | C (gcc 13.3.0) |
| 結果 |
AC
|
| 実行時間 | 4 ms / 2,000 ms |
| コード長 | 1,320 bytes |
| コンパイル時間 | 230 ms |
| コンパイル使用メモリ | 31,744 KB |
| 実行使用メモリ | 6,820 KB |
| 最終ジャッジ日時 | 2024-10-09 07:18:53 |
| 合計ジャッジ時間 | 1,154 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 27 |
ソースコード
#include <stdio.h>
const int Mod = 998244353;
int main()
{
int N, p, q;
scanf("%d %d %d", &N, &p, &q);
if (p == 1 || q == 1) {
if (p == 1 && q == 1) printf("1\n");
else printf("2\n");
fflush(stdout);
return 0;
}
int i, j, A[100001];
for (i = 1; i <= N; i++) A[i] = i;
for (i = 1, j = p; i < j; i++, j--) {
A[i] ^= A[j];
A[j] ^= A[i];
A[i] ^= A[j];
}
for (i = N - q + 1, j = N; i < j; i++, j--) {
A[i] ^= A[j];
A[j] ^= A[i];
A[i] ^= A[j];
}
int jj, n, k = 0, prime[100001], prime_inv[100001], mul[100001] = {}, flag[100001] = {};
for (i = 2; i <= N; i++) {
if (flag[i] != 0) continue;
prime[k] = i;
prime_inv[i] = k++;
for (j = i * 2; j <= N; j += i) flag[j] = 1;
}
for (i = 1; i <= N; i++) flag[i] = 0;
for (i = 1; i <= N; i++) {
if (flag[i] != 0) continue;
for (j = i, jj = A[i], flag[i] = 1; flag[jj] == 0; j = jj, jj = A[j]) flag[jj] = flag[j] + 1;
if (flag[j] == 1) continue;
for (n = flag[j], j = 0; j < k && prime[j] * prime[j] <= n; j++) {
for (jj = 0; n % prime[j] == 0; n /= prime[j], jj++);
if (mul[j] < jj) mul[j] = jj;
}
if (n > 1 && mul[prime_inv[n]] == 0) mul[prime_inv[n]] = 1;
}
long long ans = 2;
for (i = 0; i < k; i++) for (j = 0; j < mul[i]; j++) ans = ans * prime[i] % Mod;
printf("%lld\n", ans);
fflush(stdout);
return 0;
}