結果
| 問題 | No.1521 Playing Musical Chairs Alone |
| コンテスト | |
| ユーザー |
titia
|
| 提出日時 | 2026-02-19 03:06:16 |
| 言語 | PyPy3 (7.3.17) |
| 結果 |
AC
|
| 実行時間 | 1,480 ms / 2,000 ms |
| コード長 | 1,030 bytes |
| 記録 | |
| コンパイル時間 | 269 ms |
| コンパイル使用メモリ | 82,500 KB |
| 実行使用メモリ | 89,388 KB |
| 最終ジャッジ日時 | 2026-02-19 03:06:38 |
| 合計ジャッジ時間 | 19,340 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 23 |
ソースコード
import sys
input = sys.stdin.readline
mod=10**9+7
N,K,L=list(map(int,input().split()))
# 行列の計算(numpyを使えないとき,modを使用)
def prod(A,B,k,l,m):# A:k*l,B:l*m
C=[[None for i in range(m)] for j in range(k)]
for i in range(k):
for j in range(m):
ANS=0
for pl in range(l):
ANS=(ANS+A[i][pl]*B[pl][j])%mod
C[i][j]=ANS
return C
def plus(A,B,k,l):# a,B:k*l
C=[[None for i in range(l)] for j in range(k)]
for i in range(k):
for j in range(l):
C[i][j]=(A[i][j]+B[i][j])%mod
return C
A=[[0]*N for i in range(N)]
for i in range(N):
for j in range(1,L+1):
A[i][(i+j)%N]+=1
POWA=[A]
# 漸化式を行列累乗で求める(ダブリング)
for i in range(60):
POWA.append(prod(POWA[-1],POWA[-1],N,N,N)) # ベキを求めて
X=[[0]*N]
X[0][0]=1
n=K
while n:
X=prod(X,POWA[n.bit_length()-1],1,N,N) # n乗の場合
n-=1<<(n.bit_length()-1)
for a in X[0]:
print(a%mod)
titia