結果
問題 | No.1659 Product of Divisors |
ユーザー | mkawa2 |
提出日時 | 2021-08-27 23:16:34 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 67 ms / 2,000 ms |
コード長 | 1,525 bytes |
コンパイル時間 | 621 ms |
コンパイル使用メモリ | 82,004 KB |
実行使用メモリ | 77,892 KB |
最終ジャッジ日時 | 2024-05-01 04:11:03 |
合計ジャッジ時間 | 3,016 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 54 ms
75,776 KB |
testcase_01 | AC | 57 ms
76,468 KB |
testcase_02 | AC | 57 ms
76,396 KB |
testcase_03 | AC | 57 ms
76,092 KB |
testcase_04 | AC | 56 ms
76,004 KB |
testcase_05 | AC | 57 ms
75,832 KB |
testcase_06 | AC | 55 ms
77,184 KB |
testcase_07 | AC | 56 ms
76,964 KB |
testcase_08 | AC | 55 ms
76,164 KB |
testcase_09 | AC | 55 ms
76,504 KB |
testcase_10 | AC | 60 ms
76,344 KB |
testcase_11 | AC | 54 ms
75,548 KB |
testcase_12 | AC | 54 ms
75,652 KB |
testcase_13 | AC | 54 ms
75,864 KB |
testcase_14 | AC | 56 ms
75,668 KB |
testcase_15 | AC | 57 ms
76,528 KB |
testcase_16 | AC | 56 ms
76,400 KB |
testcase_17 | AC | 56 ms
75,832 KB |
testcase_18 | AC | 62 ms
75,704 KB |
testcase_19 | AC | 54 ms
75,720 KB |
testcase_20 | AC | 57 ms
76,836 KB |
testcase_21 | AC | 57 ms
77,892 KB |
testcase_22 | AC | 67 ms
76,356 KB |
testcase_23 | AC | 55 ms
76,980 KB |
testcase_24 | AC | 58 ms
77,300 KB |
ソースコード
import sys sys.setrecursionlimit(200005) int1 = lambda x: int(x)-1 p2D = lambda x: print(*x, sep="\n") def II(): return int(sys.stdin.readline()) def LI(): return list(map(int, sys.stdin.readline().split())) def LI1(): return list(map(int1, sys.stdin.readline().split())) def LLI(rows_number): return [LI() for _ in range(rows_number)] def LLI1(rows_number): return [LI1() for _ in range(rows_number)] def SI(): return sys.stdin.readline().rstrip() # dij = [(0, 1), (-1, 0), (0, -1), (1, 0)] dij = [(0, 1), (-1, 0), (0, -1), (1, 0), (1, 1), (1, -1), (-1, 1), (-1, -1)] inf = 10**16 # md = 998244353 md = 10**9+7 def nCr(com_n, com_r): if com_r < 0: return 0 if com_n < com_r: return 0 res = 1 for i in range(com_n, com_n-com_r, -1): res = res*i%md return res*ifac[com_r]%md # 準備 n_max = 200005 fac = [1] for i in range(1, n_max+1): fac.append(fac[-1]*i%md) ifac = [1]*(n_max+1) ifac[n_max] = pow(fac[n_max], md-2, md) for i in range(n_max-1, 1, -1): ifac[i] = ifac[i+1]*(i+1)%md def PrimeFactorization(x): def plist(x): if x < 2: return [] if x & 1 == 0: return [2]+plist(x >> 1) for p in range(3, x+1, 2): if x%p == 0: return [p]+plist(x//p) if p**2 > x: return [x] pl = plist(x) pp, ee = [], [] for p in pl: if not pp or p != pp[-1]: pp += [p] ee += [0] ee[-1] += 1 return ee n, k = LI() ee = PrimeFactorization(n) ans = 1 for e in ee: ans = ans*nCr(k+e, e)%md print(ans)