結果
問題 |
No.1127 変形パスカルの三角形
|
ユーザー |
|
提出日時 | 2021-03-25 22:31:16 |
言語 | Nim (2.2.0) |
結果 |
AC
|
実行時間 | 106 ms / 1,500 ms |
コード長 | 1,526 bytes |
コンパイル時間 | 4,413 ms |
コンパイル使用メモリ | 87,032 KB |
実行使用メモリ | 26,564 KB |
最終ジャッジ日時 | 2024-11-27 07:52:41 |
合計ジャッジ時間 | 8,150 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 30 |
コンパイルメッセージ
/home/judge/data/code/Main.nim(2, 18) Warning: Use the new 'sugar' module instead; future is deprecated [Deprecated] /home/judge/data/code/Main.nim(1, 8) Warning: imported and not used: 'times' [UnusedImport] /home/judge/data/code/Main.nim(2, 26) Warning: imported and not used: 'strformat' [UnusedImport] /home/judge/data/code/Main.nim(2, 18) Warning: imported and not used: 'future' [UnusedImport] /home/judge/data/code/Main.nim(2, 37) Warning: imported and not used: 'deques' [UnusedImport] /home/judge/data/code/Main.nim(1, 41) Warning: imported and not used: 'algorithm' [UnusedImport] /home/judge/data/code/Main.nim(1, 52) Warning: imported and not used: 'tables' [UnusedImport] /home/judge/data/code/Main.nim(1, 60) Warning: imported and not used: 'sets' [UnusedImport] /home/judge/data/code/Main.nim(1, 66) Warning: imported and not used: 'lists' [UnusedImport] /home/judge/data/code/Main.nim(1, 73) Warning: imported and not used: 'intsets' [UnusedImport] /home/judge/data/code/Main.nim(2, 8) Warning: imported and not used: 'critbits' [UnusedImport]
ソースコード
import times, strutils, sequtils, math, algorithm, tables, sets, lists, intsets import critbits, future, strformat, deques template `max=`(x,y) = x = max(x,y) template `min=`(x,y) = x = min(x,y) template `mod=`(x,y) = x = x mod y template scan2 = (scan(), scan()) template scan3 = (scan(), scan()) let read* = iterator: string {.closure.} = while true: (for s in stdin.readLine.split: yield s) proc scan(): int = read().parseInt proc scanf(): float = read().parseFloat proc toInt(c:char): int = return int(c) - int('0') var D = 1E9.int+7 m:int = 1E6.int fac:seq[int] finv:seq[int] inv:seq[int] proc comInit()= fac = newseqwith(m,1) finv = newseqwith(m,1) inv = newseqwith(m,1) for i in 2..<m: fac[i] = fac[i-1]*i mod D inv[i] = D - inv[D mod i] * (D div i) mod D finv[i] = finv[i-1] * inv[i] mod D proc com(n,k:int):int= if n<k:return 0 elif n<0 or k<0: return 0 else: return fac[n] * (finv[k]*finv[n-k] mod D) mod D proc H(n,r:int):int= return com(n+r-1,r) comInit() proc modPow(n,x,md:int):int= var n = n x = x result = 1 while n>0: if (n and 1)>0: result*=x result.mod=md n = n shr 1 x = x*x x.mod=md proc solve():int= var md = 10^9+7 (a,b)=(scan().mod(md),scan().mod(md)) (n,k)=(scan(),scan()) sm = 0 ca = com(n-1,k-1) cb = com(n-1,k-2) echo (ca*a+cb*b).mod(md) for i in 0..n+1: sm += ((com(n-1,i-1)*a+com(n-1,i-2)*b).mod(md))^2 sm.mod=md return sm echo solve()