結果
問題 | No.3027 f-列とh-列 |
ユーザー |
![]() |
提出日時 | 2025-02-21 21:30:42 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 39 ms / 2,000 ms |
コード長 | 764 bytes |
コンパイル時間 | 551 ms |
コンパイル使用メモリ | 82,664 KB |
実行使用メモリ | 53,608 KB |
最終ジャッジ日時 | 2025-02-21 21:30:57 |
合計ジャッジ時間 | 2,523 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 21 |
ソースコード
def cmb(n, r, p):if (r < 0) or (n < r):return 0r = min(r, n - r)return fact[n] * factinv[r] * factinv[n-r] % pp=998244353N = 60 # N は必要分だけ用意するfact = [1, 1] # fact[n] = (n! mod p)factinv = [1, 1] # factinv[n] = ((n!)^(-1) mod p)inv = [0, 1] # factinv 計算用Sfor i in range(2, N + 1):fact.append((fact[-1] * i) % p)inv.append((-inv[p % i] * (p // i)) % p)factinv.append((factinv[-1] * inv[-1]) % p)N=int(input())F=list(map(int,input().split()))H=[0]*(N+1)for i in range(N+1):for j in range(i,N+1):if (j-i)%2==0:H[j]+=F[i]*cmb(N-i,j-i,p)else:H[j]-=F[i]*cmb(N-i,j-i,p)#print(i,j,H,f"cmb({N-i},{j-i},{p})",cmb(N-i,j-i,p))print(*H)