横一列に並んだN個のマスを、N-1台のマシンで塗っていく。i番目のマシンはi番目とi+1番目のマスを黒く塗る。
ある順列Pを与え、その番号の順番でマシンを動かしていき、N個のマスが初めて黒く濡れた時点のPのindex(1-indexed)がスコアとなる。
全ての順列についてこのスコアを求め、その総和を計算する問題。
解けなかった。解法メモ。
この問題の方針は、各スコアKについて、全ての順列の内スコアKになるような通り数を求め、その合計を計算することで解を求める。
コンテスト中もこの方針で計算したが、となるを計算できず積んだ。
このでは、通り存在するスコアKの各順列について、最後の取れるマシンの番号の通り数を求め、それらの和を取る関数を求めたかった。(最後に取れないマシンの番号iは、番号i以外のK-1台のマシンでi番目とi+1番目が既に塗られている場合である)
想定解法は、K以下のスコアを取る順列の通り数()を求め、
(mod )
を解として計算する。
はで計算できる。
この通り数は、稼働するK台のマシンの間の(K-1)箇所に、稼働しないN-K-1台のマシンが高々1台ずつ入り、稼働しないマシンと稼働しないマシンの並べ方を考慮に入れることで計算できる。(どの番号のマシンが稼働するかは、この時点で考慮されている)
簡単な計算だが、逆元の計算含めで解を求めるとPythonではTLEするため、工夫して逆元の計算回数を減らしに落とす必要がある。
下記のように考えればで計算できる。下線部もループ処理で逆元なしで計算できる。
提出コード (Python3): Submission #2434722 - AtCoder Grand Contest 023
N = int(input()) MOD = 10**9 + 7 fact = [1]*(N+1) for i in range(1, N+1): fact[i] = r = i*fact[i-1] % MOD cnts = [0]*(N+1) rev = 1 # rev: 上記の下線部部分 for K in range(N-1, (N+1)//2-1, -1): cnts[K] = fact[K]*fact[K-1]*rev % MOD rev = rev * (2*K-N) * (2*K-N-1) % MOD ans = 0 for K in range((N+1)//2, N): ans += (cnts[K] - cnts[K-1]) * K % MOD ans %= MOD # 最後に一回逆元計算 ans = (ans * pow(fact[N-2], MOD-2, MOD)) % MOD print(ans)