結果

問題 No.271 next_permutation (2)
ユーザー gew1fw
提出日時 2025-06-12 20:59:47
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 1,646 bytes
コンパイル時間 433 ms
コンパイル使用メモリ 82,004 KB
実行使用メモリ 91,612 KB
最終ジャッジ日時 2025-06-12 21:03:18
合計ジャッジ時間 6,814 ms
ジャッジサーバーID
(参考情報)
judge2 / judge5
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other WA * 2 TLE * 1 -- * 18
権限があれば一括ダウンロードができます

ソースコード

diff #

import sys
import math
from itertools import permutations

MOD = 1000000007

def compute_inversion(p):
    n = len(p)
    inv = 0
    for i in range(n):
        for j in range(i+1, n):
            if p[i] > p[j]:
                inv += 1
    return inv

def main():
    N, K = map(int, sys.stdin.readline().split())
    p = list(map(int, sys.stdin.readline().split()))
    
    if K == 0:
        print(0)
        return
    
    if N < 20:
        m = math.factorial(N)
        perms = list(permutations(range(1, N+1)))
        inv_list = [compute_inversion(perm) for perm in perms]
        sum_cycle = sum(inv_list)
        prefix = [0] * (m + 1)
        for i in range(m):
            prefix[i+1] = (prefix[i] + inv_list[i]) % MOD
        
        target_perm = tuple(p)
        try:
            start_index = perms.index(target_perm)
        except ValueError:
            print(0)
            return
        
        sum_per = []
        current = start_index
        for i in range(m):
            sum_per.append(inv_list[current])
            current = (current + 1) % m
        
        full_cycles = K // m
        remainder = K % m
        total = (full_cycles * sum_cycle) % MOD
        total = (total + sum(sum_per[:remainder])) % MOD
        print(total)
    else:
        C = N * (N-1) // 2
        inv_initial = compute_inversion(p)
        if K == 0:
            print(0)
        else:
            inv_total = (K % MOD) * C
            inv_total = inv_total % MOD
            inv_total = inv_total * pow(2, MOD-2, MOD)
            inv_total = inv_total % MOD
            print(inv_total)
    
if __name__ == "__main__":
    main()
0