#758208 (Lua) No.1330 Multiply or Divide

• 問題: No.1330 Multiply or Divide
• ユーザー: 👑 obakyan
• 提出日時: 2022-05-05 19:32:59
• 言語: Lua (LuaJit 2.1.1734355927)
• 結果: AC
• 実行時間: 1,209 ms / 2,000 ms
• コード長: 2,754 bytes
• コンパイル時間: 134 ms
• コンパイル使用メモリ: 5,120 KB
• 実行使用メモリ: 31,660 KB
• 最終ジャッジ日時: 2024-07-04 21:12:44
• 合計ジャッジ時間: 5,847 ms
• ジャッジサーバーID （参考情報）: judge2 / judge1

ファイルパターン	結果
other	AC * 46

ソースコード

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
local mmi, mma = math.min, math.max
local mfl, mce = math.floor, math.ceil
local bls, brs = bit.lshift, bit.rshift
local Heapq = {}
Heapq.create = function(self, lt)
  self.lt = lt
  self.cnt = 0
  self.t = {}
end
Heapq.empty = function(self)
  return self.cnt == 0
end
Heapq.push = function(self, v)
  local hqlt = self.lt
  local hqt = self.t
  local c = self.cnt + 1
  self.cnt = c
  hqt[c] = v
  while 1 < c do
    local p = brs(c, 1)
    if hqlt(hqt[c], hqt[p]) then
      hqt[c], hqt[p] = hqt[p], hqt[c]
      c = p
    else
      break
    end
  end
end
Heapq.pop = function(self)
  local hqlt = self.lt
  local hqt = self.t
  local ret = hqt[1]
  local c = self.cnt
  hqt[1] = hqt[c]
  c = c - 1
  self.cnt = c
  local p = 1
  while true do
    local d1, d2 = p * 2, p * 2 + 1
    if c < d1 then break
    elseif c < d2 then
      if hqlt(hqt[d1], hqt[p]) then
        hqt[d1], hqt[p] = hqt[p], hqt[d1]
      end
      break
    else
      if hqlt(hqt[d1], hqt[d2]) then
        if hqlt(hqt[d1], hqt[p]) then
          hqt[d1], hqt[p] = hqt[p], hqt[d1]
          p = d1
        else break
        end
      else
        if hqlt(hqt[d2], hqt[p]) then
          hqt[d2], hqt[p] = hqt[p], hqt[d2]
          p = d2
        else break
        end
      end
    end
  end
  return ret
end
Heapq.new = function(lt)
  local obj = {}
  setmetatable(obj, {__index = Heapq})
  obj:create(lt)
  return obj
end
local n, m, p = io.read("*n", "*n", "*n")
local a = {}
local amax = 0
for i = 1, n do
  a[i] = io.read("*n")
  amax = mma(amax, a[i])
end
if m < amax then print(1) os.exit() end
local c = {}
for i = 1, n do
  local ci = 1
  local ai = a[i]
  while ai % p == 0 do
    ci = ci + 1
    ai = mfl(ai / p)
  end
  if 1 < ai then
    if c[ai] then
      c[ai] = mmi(c[ai], ci)
    else
      c[ai] = ci
    end
  end
end
if not next(c) then print(-1) os.exit() end
local cm = {}
for k, v in pairs(c) do
  table.insert(cm, {k, v})
end
table.sort(cm, function(a, b) return a[1] < b[1] end)
local cn = #cm
local val, score, valmap = {}, {}, {}
local lim = mfl(m / amax) + 1
valmap[m] = 1
val[1] = 1
score[1] = 1
local hq = Heapq.new(function(a, b) return val[a] < val[b] end)
hq:push(1)
local ret = 1000000007
while not hq:empty() do
  local idx = hq:pop()
  local v, s = val[idx], score[idx]
  for i = 1, cn do
    local mul = cm[i][1]
    local cost = cm[i][2]
    if lim <= v * mul then
      ret = mmi(ret, s + cost)
    else
      local dst = v * mul
      local dstidx = valmap[dst]
      if dstidx then
        score[dstidx] = mmi(score[dstidx], s + cost)
      else
        table.insert(val, dst)
        table.insert(score, s + cost)
        valmap[dst] = #val
        hq:push(#val)
      end
    end
  end
end
print(ret)
-- print(os.clock())
 
 
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX