結果

問題 No.2963 Mecha DESU
ユーザー 遭難者
提出日時 2024-10-26 00:58:58
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 186 ms / 2,000 ms
コード長 574 bytes
コンパイル時間 520 ms
コンパイル使用メモリ 82,368 KB
実行使用メモリ 105,056 KB
最終ジャッジ日時 2024-10-26 11:01:44
合計ジャッジ時間 8,893 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 57
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys

input = sys.stdin.readline
n, m, k = map(int, input().split())
a = list(map(int, input().split()))
d0 = [0 for i in range(n + 1)]
for i in a:
    if i <= n:
        d0[i] += 1
d1 = [0 for i in range(n + 1)]
for i in range(1, n + 1):
    for j in range(i, n + 1, i):
        d1[j] += d0[i]
ans = 0
mod = 998244353
m_inv = pow(m, mod - 2, mod)
m_cnt = [0 for i in range(m + 1)]
for i in range(1, n + 1):
    m_cnt[d1[i]] += 1
for i in range(1, m + 1):
    if m_cnt[i] == 0:
        continue
    ans += (1 - pow((m - i) * m_inv, k, mod)) * m_cnt[i]
print(ans % mod)
0