結果

問題 No.573 a^2[i] = a[i]
コンテスト
ユーザー 6soukiti29
提出日時 2017-10-06 23:50:22
言語 Nim
(2.2.0)
結果
AC  
実行時間 42 ms / 2,000 ms
コード長 1,065 bytes
コンパイル時間 3,166 ms
コンパイル使用メモリ 66,816 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-06-30 03:15:58
合計ジャッジ時間 4,794 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 47
権限があれば一括ダウンロードができます
コンパイルメッセージ
/home/judge/data/code/Main.nim(1, 8) Warning: imported and not used: 'sequtils' [UnusedImport]
/home/judge/data/code/Main.nim(1, 26) Warning: imported and not used: 'strscans' [UnusedImport]

ソースコード

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 = 100_003
    m : int64 = 1_000_000_007

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 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

var
    N = stdin.readline.parseInt
    ans : int64
    
for p in 1..N:
    ans = (ans + combination(N, p) * powInt(p, N - p)) mod m

echo ans
0