結果
問題 |
No.2136 Dice Calendar?
|
ユーザー |
![]() |
提出日時 | 2022-10-14 12:40:41 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 4,427 ms / 5,000 ms |
コード長 | 1,882 bytes |
コンパイル時間 | 485 ms |
コンパイル使用メモリ | 82,432 KB |
実行使用メモリ | 505,512 KB |
最終ジャッジ日時 | 2024-06-26 12:35:11 |
合計ジャッジ時間 | 25,093 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 24 |
ソースコード
N = int(input()) S = [list(map(lambda x: int(x) - 1, input().split())) for _ in range(N)] factorial = [1] for i in range(N): factorial.append(factorial[i] * (i + 1)) def getPartition(diceSet, index): # 多重集合を表すbitでi番目に立っているbitの位置を求める return diceSet[1] >> 5 * index & 0b11111 def multichoose(diceSet): # この多重集合を並べてできる組合せ ret = factorial[getPartition(diceSet, 9) - 9] for i in range(9): ret //= factorial[getPartition(diceSet, i + 1) - getPartition(diceSet, i) - 1] return ret def nextSet(diceSet, dice, uniqueCheck, nextQueue): # diceを追加したときの多重集合をnextQueueに入れる for result in dice: # 出目がresultだった時 mask = (1 << getPartition(diceSet, result)) - 1 nextSet = (diceSet[0] & 0x1FFFFFFF - mask) << 1 | diceSet[0] & mask if (uniqueCheck[nextSet >> 6] >> (nextSet & 0x3F) & 1) == 0: # まだこの多重集合を計算対象にしていないなら uniqueCheck[nextSet >> 6] |= 1 << (nextSet & 0x3F) nextPartition = diceSet[1] + (0b00001_00001_00001_00001_00001_00001_00001_00001_00001_00001 & 0x3FFFFFFFFFFFF - ((0x20 << result * 5) - 1)) nextQueue.append((nextSet, nextPartition)) # (多重集合を、仕切りの考え方で見なした時のbit列, 上のbit列で立っているbitの位置)の二要素を状態とする nowQueue = [(0b111111111, 0b01001_01000_00111_00110_00101_00100_00011_00010_00001_00000)] # 初項M_0を求める uniqueCheck = [0] * (1 << N + 3) # 既に調べた多重集合を管理するためのBitSet、64MB程度 for dice in S: for i in nowQueue: uniqueCheck[i[0] >> 6] = 0 nextQueue = [] for diceSet in nowQueue: nextSet(diceSet, dice, uniqueCheck, nextQueue) # M_iからM_{i+1}を求める nowQueue = nextQueue ans = 0 for diceSet in nowQueue: ans += multichoose(diceSet) ans %= 998_244_353 print(ans)