結果

問題 No.187 中華風 (Hard)
ユーザー kotatsugamekotatsugame
提出日時 2020-09-17 04:11:22
言語 Ruby
(2.7.1p83)
結果
AC  
実行時間 52 ms / 3,000 ms
コード長 1,039 Byte
コンパイル時間 33 ms
使用メモリ 15,348 KB
最終ジャッジ日時 2020-09-17 04:11:25
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
使用メモリ
testcase_00 AC 13 ms
11,496 KB
testcase_01 AC 12 ms
11,572 KB
testcase_02 AC 25 ms
12,848 KB
testcase_03 AC 25 ms
12,744 KB
testcase_04 AC 44 ms
15,052 KB
testcase_05 AC 44 ms
15,124 KB
testcase_06 AC 47 ms
15,204 KB
testcase_07 AC 52 ms
15,160 KB
testcase_08 AC 52 ms
15,204 KB
testcase_09 AC 45 ms
15,300 KB
testcase_10 AC 47 ms
15,272 KB
testcase_11 AC 45 ms
15,164 KB
testcase_12 AC 45 ms
15,172 KB
testcase_13 AC 13 ms
11,560 KB
testcase_14 AC 13 ms
11,488 KB
testcase_15 AC 26 ms
12,860 KB
testcase_16 AC 26 ms
13,020 KB
testcase_17 AC 11 ms
11,456 KB
testcase_18 AC 13 ms
11,600 KB
testcase_19 AC 10 ms
11,572 KB
testcase_20 AC 38 ms
14,764 KB
testcase_21 AC 11 ms
11,552 KB
testcase_22 AC 45 ms
15,156 KB
testcase_23 AC 11 ms
11,452 KB
testcase_24 AC 11 ms
11,452 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
Main.rb:60: warning: ambiguous first argument; put parentheses or a space even after `-' operator
Syntax OK

ソースコード

diff #
# frozen_string_literal: true

# return [rem, mod] or [0, 0] (if no solution)
def crt(r, m)
  raise ArgumentError if r.size != m.size
  
  n = r.size
  r0, m0 = 0, 1
  n.times{ |i|
    raise ArgumentError if m[i] < 1

    r1, m1 = r[i]%m[i], m[i]
    if m0 < m1 then
      r0, r1 = r1, r0
      m0, m1 = m1, m0
    end

    if m0%m1 == 0 then
      return [0, 0] if r0%m1 != r1
      next
    end

    g, im = inv_gcd(m0, m1)
    u1 = m1/g
    return [0, 0] if (r1-r0)%g != 0

    x = (r1-r0)/g*im%u1
    r0 += x*m0
    m0 *= u1
    r0 += m0 if r0 < 0
  }

  return [r0, m0]
end

# internal method
# return [g, x] s.t. g = gcd(a, b), x*a = g (mod b), 0 <= x < b/g
def inv_gcd(a, b)
  a %= b
  return [b, 0] if a == 0

  s, t = b, a
  m0, m1 = 0, 1
  while t>0 do
    u, s = s.divmod(t)
    m0 -= m1*u
    s, t = t, s
    m0, m1 = m1, m0
  end
  m0 += b/s if m0 < 0
  return [s, m0]
end

gets
r, m = $<.map{ |e| e.split.map &:to_i }.transpose

r, m = crt(r, m)

if r == 0 && m == 0
  puts -1
else
  r = m if r == 0
  puts r % (10**9+7)
end
0