結果

問題 No.117 組み合わせの数
ユーザー 6soukiti296soukiti29
提出日時 2017-09-01 07:17:44
言語 Nim
(2.0.2)
結果
CE  
(最新)
AC  
(最初)
実行時間 -
コード長 1,500 bytes
コンパイル時間 844 ms
コンパイル使用メモリ 66,100 KB
最終ジャッジ日時 2024-11-14 20:14:36
合計ジャッジ時間 1,289 ms
ジャッジサーバーID
(参考情報)
judge2 / judge1
このコードへのチャレンジ
(要ログイン)
コンパイルエラー時のメッセージ・ソースコードは、提出者また管理者しか表示できないようにしております。(リジャッジ後のコンパイルエラーは公開されます)
ただし、clay言語の場合は開発者のデバッグのため、公開されます。

コンパイルメッセージ
/home/judge/data/code/Main.nim(24, 1) Error: index 3000002 not in 0 .. 3000001

ソースコード

diff #

import sequtils,strutils,strscans

proc powInt(n : int64, m : int64, k = 1_000_000_007):int64 =
    if m == 0:
        return 1
    elif m == 1:
        return (n mod k)
    if (m mod 2) == 0:
        return powInt((n * n) mod k,m div 2, k) mod k
    else:
        return (powInt((n * n) mod k,m div 2, k) * n) mod k

const
    Max_num = 3_000_002
    m : int64 = 1_000_000_007

var
    kaizyou : array[Max_num, int64]
    gyakuKaizyou : array[Max_num, int64]
kaizyou[0] = 1

for i in 1..Max_num:
    kaizyou[i] = ((kaizyou[i - 1] * i.int64) mod m)
gyakuKaizyou[Max_num] = powInt(kaizyou[Max_num], 1_000_000_005, 1_000_000_007)

for i in countdown(Max_num - 1,0):
    gyakuKaizyou[i] = (gyakuKaizyou[i + 1] * (i + 1)) mod m

proc combination(n : int64, k : int64): int64 =
    if n < k:
        return 0
    if n == k:
        return 1
    return (((kaizyou[n] * gyakuKaizyou[n - k]) mod m) * gyakuKaizyou[k]) mod m
    
proc permutation(n : int64, k : int64): int64 =
    if n < k:
        return 0
    if n == k:
        return kaizyou[n]
    return (kaizyou[n] * gyakuKaizyou[n - k]) mod m

var
    s : string
    T = stdin.readline.parseInt
    c : string
    a,b : int
for t in 1..T:
    s = stdin.readline
    discard scanf(s, "$w($i,$i)", c, a, b)
    if c == "C":
        echo combination(a.int64, b.int64)
    if c == "P":
        echo permutation(a.int64, b.int64)
    if c == "H":
        if b == 0:
            echo 1
        else:
            echo combination((a + b - 1).int64, b.int64)
0