結果

問題 No.2975 単調増加部分積
ユーザー katonyonkokatonyonko
提出日時 2024-11-29 23:16:17
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 1,834 bytes
コンパイル時間 247 ms
コンパイル使用メモリ 82,304 KB
実行使用メモリ 852,224 KB
最終ジャッジ日時 2024-11-29 23:16:31
合計ジャッジ時間 13,718 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 52 ms
55,040 KB
testcase_01 AC 48 ms
54,400 KB
testcase_02 AC 46 ms
54,144 KB
testcase_03 AC 46 ms
54,784 KB
testcase_04 WA -
testcase_05 WA -
testcase_06 AC 47 ms
54,400 KB
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 AC 62 ms
65,280 KB
testcase_11 WA -
testcase_12 WA -
testcase_13 AC 137 ms
88,960 KB
testcase_14 AC 81 ms
77,696 KB
testcase_15 AC 179 ms
99,072 KB
testcase_16 AC 60 ms
64,896 KB
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
testcase_20 MLE -
testcase_21 MLE -
testcase_22 MLE -
testcase_23 MLE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

import io
import sys
from collections import defaultdict, deque, Counter
from itertools import permutations, combinations, accumulate
from heapq import heappush, heappop
from bisect import bisect_right, bisect_left
from math import gcd
import math

_INPUT = """\
6
1 1 2
2 1 3
2 2 5
"""

def input():
  return sys.stdin.readline()[:-1]

def solve(test):
  N, M, P = map(int, input().split())
  flg = 0
  for i in range(N-M+1, N+1):
    if gcd(i, P)!=1:
      flg = 1
      break
  if flg==1:
    print(-1)
  else:
    dp=[None]*M*N
    def idx(i,j):
      return i*N+j
    for j in reversed(range(N)):
      dp[idx(M-1,j)]=(j+1)%P
    adp=list(accumulate([dp[idx(M-1,j)] for j in reversed(range(N))]))[::-1]
    adp=[adp[i]%P for i in range(N)]
    for i in reversed(range(M-1)):
      for j in reversed(range(N)):
        dp[idx(i,j)]=(j+1)*((adp[j+1] if j!=N-1 else 0) +1)%P
      tmp=list(accumulate([dp[idx(i,j)] for j in reversed(range(N))]))[::-1]
      adp=[(tmp[j]+adp[j])%P for j in range(N)]
    ans=sum(dp)%P
    for i in range(N-M+1,N+1):
      ans=(ans*pow(i,P-2,P))%P
    print(ans)

def random_input():
  from random import randint,shuffle
  N=randint(1,10)
  M=randint(1,N)
  A=list(range(1,M+1))+[randint(1,M) for _ in range(N-M)]
  shuffle(A)
  return (" ".join(map(str, [N,M]))+"\n"+" ".join(map(str, A))+"\n")*3

def simple_solve():
  return []

def main(test):
  if test==0:
    solve(0)
  elif test==1:
    sys.stdin = io.StringIO(_INPUT)
    case_no=int(input())
    for _ in range(case_no):
      solve(0)
  else:
    for i in range(1000):
      sys.stdin = io.StringIO(random_input())
      x=solve(1)
      y=simple_solve()
      if x!=y:
        print(i,x,y)
        print(*[line for line in sys.stdin],sep='')
        break

#0:提出用、1:与えられたテスト用、2:ストレステスト用
main(0)
0