結果

問題 No.2318 Phys Bone Maker
ユーザー cleanttedcleantted
提出日時 2023-02-26 02:56:00
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 2,143 ms / 3,000 ms
コード長 1,934 bytes
コンパイル時間 309 ms
コンパイル使用メモリ 82,348 KB
実行使用メモリ 94,620 KB
最終ジャッジ日時 2024-09-18 12:56:21
合計ジャッジ時間 22,803 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 45
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
from collections import deque, Counter
from math import gcd, log10, pi, sqrt, ceil
from bisect import bisect, bisect_left, bisect_right, insort
from typing import Iterable, TypeVar, Union, Tuple, Iterator, List
import copy
import heapq
import itertools
import math
import random
from fractions import Fraction
from functools import lru_cache, partial, cmp_to_key
import operator
input = sys.stdin.readline
sys.setrecursionlimit(10000000)
mod = 998244353
INF = 1 << 61
DIFF = 10 ** -9
DX = [1, 0, -1, 0, 1, 1, -1, -1]
DY = [0, 1, 0, -1, 1, -1, 1, -1]

def read_values(): return map(int, input().split())
def read_index(): return map(lambda x: int(x) - 1, input().split())
def read_list(): return list(read_values())
def read_lists(N): return [read_list() for _ in range(N)]


def main():
    N = int(input())

    T = {}
    for n in range(2, int(N ** 0.5) + 1):
        while N % n == 0:
            T[n] = T.get(n, 0) + 1
            N //= n

    if N > 1:
        T[N] = T.get(N, 0) + 1 

    T = list(T.items())
    X = {}
    V = set()

    def f(s, S):
        if len(S) == len(T):
            V.add(s)
            X[s] = (tuple(S))
            return
        
        k, v = T[len(S)] 
        S.append(0)
        for _ in range(v + 1):
            f(s, S)
            s *= k
            S[-1] += 1
        S.pop()
        s //= k ** v
    
    f(1, [])

    N = len(V)
    V = sorted(list(V))
    X = [X[v] for v in V]
    dp = [0] * N
    dp[0] = 1
    for i in range(N):
        for j in range(i + 1, N):
            r = 1
            for x1, x2 in zip(X[i], X[j]):
                if x2 > x1:
                    continue
                elif x2 == x1:
                    r *= x1 + 1
                    r %= mod
                else:
                    break
            else:
                dp[j] += dp[i] * r % mod
                dp[j] %= mod
    
    print(dp[-1])


if __name__ == "__main__":
    main()
0