結果

問題 No.389 ロジックパズルの組み合わせ
コンテスト
ユーザー むらため
提出日時 2019-01-27 18:46:56
言語 Nim
(2.2.0)
結果
WA  
(最新)
AC  
(最初)
実行時間 -
コード長 2,387 bytes
コンパイル時間 2,609 ms
コンパイル使用メモリ 65,348 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-07-01 11:03:04
合計ジャッジ時間 5,066 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 49 WA * 50
権限があれば一括ダウンロードができます
コンパイルメッセージ
/home/judge/data/code/Main.nim(1, 8) Warning: imported and not used: 'sequtils' [UnusedImport]
/home/judge/data/code/Main.nim(1, 32) Warning: imported and not used: 'strutils' [UnusedImport]

ソースコード

diff # 

import sequtils,algorithm,math,strutils
template times*(n:int,body) = (for _ in 0..<n: body)
template `max=`*(x,y) = x = max(x,y)
template `min=`*(x,y) = x = min(x,y)

template useModulo() =
  const MOD = 1_000_000_007
  type ModInt = object
    v:int # 0~MODに収まる
  proc toModInt*(a:int) : ModInt =
    if a < -MOD : result.v = ((a mod MOD) + MOD) mod MOD
    elif a < 0 : result.v = a + MOD
    elif a >= MOD: result.v = a mod MOD
    else: result.v = a
  proc `+`*(a,b:ModInt) : ModInt =
    result.v = a.v + b.v
    if result.v >= MOD : result.v = result.v mod MOD
  proc `*`*(a,b:ModInt) : ModInt =
    result.v = a.v * b.v
    if result.v >= MOD : result.v = result.v mod MOD
  proc `^`*(a:ModInt,b:int) : ModInt =
    if b == 0 : return 1.toModInt()
    if b == 1 : return a
    let pow = a^(b div 2)
    if b mod 2 == 0 : return pow * pow
    return pow * pow * a
  proc `+`*(a:int,b:ModInt) : ModInt = a.toModInt() + b
  proc `+`*(a:ModInt,b:int) : ModInt = a + b.toModInt()
  proc `-`*(a:ModInt,b:int) : ModInt = a + (-b)
  proc `-`*(a,b:ModInt) : ModInt = a + (-b.v)
  proc `-`*(a:int,b:ModInt) : ModInt = a.toModInt() + (-b.v)
  proc `*`*(a:int,b:ModInt) : ModInt = a.toModInt() * b
  proc `*`*(a:ModInt,b:int) : ModInt = a * b.toModInt()
  proc `/`*(a,b:ModInt) : ModInt = a * b^(MOD-2)
  proc `$`*(a:ModInt) : string = $a.v

  proc combination(n,k:int) : ModInt = # nCk を剰余ありで
    result = 1.toModInt()
    let x = k.max(n - k)
    let y = k.min(n - k)
    var fact = 1.toModInt()
    for i in 2..y: fact = fact * i
    for i in 1..y: result = result * (n+1-i)
    result = result / fact
useModulo()

proc getchar_unlocked():char {. importc:"getchar_unlocked",header: "<stdio.h>" .}
proc scan(): int =
  while true:
    let k = getchar_unlocked()
    if k < '0': break
    result = 10 * result + k.ord - '0'.ord
proc scanH(): tuple[sum,size:int] =
  result = (0,0)
  while true:
    var now = 0
    while true:
      let k = getchar_unlocked()
      if k < '0':
        result.size += 1
        result.sum += now
        if k == '\n': return
        break
      now = 10 * now + k.ord - '0'.ord

let m = scan()
let (hSum,hLen) = scanH()
if hLen == 1 :
  if hSum == 0: quit "1",0
let n = m - (hSum - hLen)
let h = hLen
let k = n - (h + h - 1)
if k < 0: quit "NA",0
if k == 0 : quit "1",0
# h箇所にk個入れる組み合わせ
echo (h+k).combination(h)
0