結果

問題 No.140 みんなで旅行
ユーザー 6soukiti296soukiti29
提出日時 2019-09-19 10:39:52
言語 Nim
(2.0.2)
結果
AC  
実行時間 90 ms / 5,000 ms
コード長 1,584 bytes
コンパイル時間 3,322 ms
コンパイル使用メモリ 65,152 KB
実行使用メモリ 36,608 KB
最終ジャッジ日時 2024-07-18 14:31:06
合計ジャッジ時間 5,173 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 29 ms
34,304 KB
testcase_01 AC 32 ms
34,304 KB
testcase_02 AC 32 ms
34,816 KB
testcase_03 AC 31 ms
34,304 KB
testcase_04 AC 31 ms
34,304 KB
testcase_05 AC 32 ms
34,432 KB
testcase_06 AC 29 ms
34,304 KB
testcase_07 AC 32 ms
34,432 KB
testcase_08 AC 32 ms
34,304 KB
testcase_09 AC 31 ms
34,432 KB
testcase_10 AC 32 ms
34,432 KB
testcase_11 AC 90 ms
36,608 KB
testcase_12 AC 36 ms
34,688 KB
testcase_13 AC 35 ms
34,432 KB
testcase_14 AC 88 ms
36,608 KB
testcase_15 AC 90 ms
36,608 KB
testcase_16 AC 52 ms
35,712 KB
testcase_17 AC 40 ms
35,456 KB
testcase_18 AC 76 ms
36,480 KB
testcase_19 AC 84 ms
36,480 KB
testcase_20 AC 39 ms
35,200 KB
testcase_21 AC 29 ms
34,432 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
/home/judge/data/code/Main.nim(1, 8) Warning: imported and not used: 'sequtils' [UnusedImport]

ソースコード

diff # 

import sequtils,strutils
const
    Max_num = 2_000_003
    modnum : int64 = 1_000_000_007

proc powInt(n : int64, m : int64, k = modnum):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

var
    kaizyou : array[Max_num + 1, int64]
    gyakuKaizyou : array[Max_num + 1, int64]
    
kaizyou[0] = 1
for i in 1..Max_num:
    kaizyou[i] = ((kaizyou[i - 1] * i.int64) mod modnum)
    
gyakuKaizyou[Max_num] = powInt(kaizyou[Max_num], modnum - 2)
for i in countdown(Max_num - 1,0):
    gyakuKaizyou[i] = (gyakuKaizyou[i + 1] * (i + 1)) mod modnum

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

var
    N = stdin.readline.parseInt
    memo : array[560, array[560, int64]]
    m = 1000_000_007
    ans : int64
    
proc F(x,y : int64):int64 =
    if x == 0 and y == 0:
        return 1
    if x < y or y <= 0:
        return 0
    if memo[x][y] > 0:
        return memo[x][y]
    
    memo[x][y] = (F(x - 1, y - 1) mod m + y * F(x - 1, y) mod m) mod m
    return memo[x][y]
        
for y in 1..(N.int64):
    for x in 1..(N.int64):
        var
            A = 1.int64
        A = combination(N, x) * F(x, y)
        A = A mod m
        A *= powInt((y * (y - 1) mod m), N.int64 - x)
        A = A mod m
        ans += A
        ans = ans mod m
echo ans
0