結果
| 問題 |
No.1973 Divisor Sequence
|
| コンテスト | |
| ユーザー |
vwxyz
|
| 提出日時 | 2023-11-21 04:34:26 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 1,283 bytes |
| コンパイル時間 | 1,054 ms |
| コンパイル使用メモリ | 82,048 KB |
| 実行使用メモリ | 92,100 KB |
| 最終ジャッジ日時 | 2024-09-26 06:56:57 |
| 合計ジャッジ時間 | 14,494 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 21 TLE * 1 |
ソースコード
import bisect
import copy
import decimal
import fractions
import heapq
import itertools
import math
import random
import sys
import time
from collections import Counter,deque,defaultdict
from functools import lru_cache,reduce
from heapq import heappush,heappop,heapify,heappushpop,_heappop_max,_heapify_max
def _heappush_max(heap,item):
heap.append(item)
heapq._siftdown_max(heap, 0, len(heap)-1)
def _heappushpop_max(heap, item):
if heap and item < heap[0]:
item, heap[0] = heap[0], item
heapq._siftup_max(heap, 0)
return item
from math import gcd as GCD
read=sys.stdin.read
readline=sys.stdin.readline
readlines=sys.stdin.readlines
write=sys.stdout.write
#import pypyjit
#pypyjit.set_param('max_unroll_recursion=-1')
#sys.set_int_max_str_digits(10**9)
def Factorize(N):
assert N>=1
factors=defaultdict(int)
for p in range(2,N):
if p**2>N:
break
while N%p==0:
factors[p]+=1
N//=p
if N!=1:
factors[N]+=1
return factors
N,M=map(int,readline().split())
ans=1
mod=10**9+7
for p,e in Factorize(M).items():
dp=[0]*(e+1)
dp[0]=1
for i in range(N):
prev=dp
dp=[sum(prev[:e+1-i])%mod for i in range(e+1)]
ans*=sum(dp)%mod
ans%=mod
print(ans)
vwxyz