結果

問題 No.1973 Divisor Sequence
コンテスト
ユーザー TKG1
提出日時 2022-06-10 23:47:25
言語 Crystal
(1.14.0)
結果
TLE  
実行時間 -
コード長 865 bytes
コンパイル時間 12,500 ms
コンパイル使用メモリ 298,400 KB
実行使用メモリ 249,216 KB
最終ジャッジ日時 2024-09-21 07:02:44
合計ジャッジ時間 16,072 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2 TLE * 1
other -- * 22
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

alias LL = Int64
macro chmin(x, y); {{x}} = {{y}} if {{x}} > {{y}}; end
macro chmax(x, y); {{x}} = {{y}} if {{x}} < {{y}}; end
macro yn(x); {{x}} ? "Yes" : "No"; end

# 約数列挙
def divsors(n : Int64)
  res = Set(Int64).new
  i = 1i64
  while i * i <= n
    if n % i == 0
      res << i
      res << (n // i)
    end
    i += 1
  end
  return res
end

# ---------------------------------------------------- :)

n, m = read_line.split.map(&.to_i64)
mod = 10**9+7
div = divsors(m)
dp = Array.new(n){Hash(LL, LL).new(0)}
sum = Hash(LL, Array(LL)).new

div.each do |i|
  dp[0][i] = 1
  div.each do |j|
    if (1i128*i*j) % m == 0
      sum[i] ||= [] of Int64
      sum[i] << j
    end
  end
end

1.upto(n-1) do |i|
  div.each do |j|
    sum[j].each do |k|
      dp[i][j] += dp[i-1][k]
      dp[i][j] %= mod
    end
  end
end

puts dp[n-1].values.sum(&.% mod) % mod
0