結果

問題 No.2136 Dice Calendar?
ユーザー CuriousFairy315CuriousFairy315
提出日時 2022-11-23 17:05:51
言語 C#(csc)
(csc 3.9.0)
結果
AC  
実行時間 1,769 ms / 5,000 ms
コード長 2,798 bytes
コンパイル時間 1,789 ms
コンパイル使用メモリ 116,440 KB
実行使用メモリ 110,760 KB
最終ジャッジ日時 2024-09-24 20:52:24
合計ジャッジ時間 12,886 ms
ジャッジサーバーID
(参考情報)
judge4 / judge2
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 31 ms
27,176 KB
testcase_01 AC 30 ms
27,216 KB
testcase_02 AC 38 ms
25,428 KB
testcase_03 AC 30 ms
25,008 KB
testcase_04 AC 29 ms
25,260 KB
testcase_05 AC 31 ms
25,272 KB
testcase_06 AC 31 ms
25,384 KB
testcase_07 AC 32 ms
25,392 KB
testcase_08 AC 33 ms
25,220 KB
testcase_09 AC 37 ms
27,092 KB
testcase_10 AC 41 ms
27,008 KB
testcase_11 AC 59 ms
25,264 KB
testcase_12 AC 69 ms
29,144 KB
testcase_13 AC 59 ms
25,316 KB
testcase_14 AC 84 ms
29,616 KB
testcase_15 AC 248 ms
37,648 KB
testcase_16 AC 313 ms
39,356 KB
testcase_17 AC 239 ms
51,872 KB
testcase_18 AC 702 ms
53,812 KB
testcase_19 AC 839 ms
76,144 KB
testcase_20 AC 1,005 ms
76,316 KB
testcase_21 AC 1,300 ms
104,356 KB
testcase_22 AC 1,567 ms
107,808 KB
testcase_23 AC 1,769 ms
110,384 KB
testcase_24 AC 77 ms
92,316 KB
testcase_25 AC 113 ms
92,460 KB
testcase_26 AC 1,582 ms
110,760 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
Microsoft (R) Visual C# Compiler version 3.9.0-6.21124.20 (db94f4cc)
Copyright (C) Microsoft Corporation. All rights reserved.

ソースコード

diff #

using System;
using System.Collections;
using System.Collections.Generic;
using System.Linq;

class Program
{
	static void Main(string[] args)
	{
		int N = int.Parse(Console.ReadLine());
		int[][] S = new int[N][];
		for (int i = 0;i < N;++ i) {
			S[i] = Console.ReadLine().Split(' ').Select(x => int.Parse(x) - 1).ToArray();
		}
		
		const int MOD = 998_244_353;
		long[] factorial = new long[N + 1];
		factorial[0] = 1;
		for (int i = 1;i < factorial.Length;++ i) factorial[i] = factorial[i - 1] * i;
		
		BitArray set = new BitArray(1 << N + 9); // 64MB程度
		Stack<int> now = new Stack<int>(3200000), next = new Stack<int>(3200000);
		now.Push(0b111111111); // 初項M_0を求める
		foreach (int[] dice in S) {
			foreach (int multiSet in now) nextSet(multiSet, dice, set, next); // M_iからM_{i+1}を求める
			Stack<int> swap = now;
			now = next;
			next = swap;
			next.Clear();
		}
		
		int ans = 0;
		foreach (int multiSet in now) ans = (int) ((ans + multichoose(factorial, multiSet)) % MOD);
		Console.WriteLine(ans);
	}
	
	static readonly int[] deBrujin32 = { 0, 1, 28, 2, 29, 14, 24, 3, 30, 22, 20, 15, 25, 17, 4, 8, 31, 27, 13, 23, 21, 19, 16, 7, 26, 12, 18, 6, 11, 5, 10, 9 };
	
	static int numberOfTrailingZeros(int i) { // 丁度1bit立っている値に対してその立っている位置を返す
		if (i == 0) return 32;
		return deBrujin32[(uint)(i * 0x077CB531) >> 27];
	}
	
	static long calcPartition(int multiSet) { // 与えられた多重集合に対して、立っているbitの位置を保持する数列Pを返す
		long partition = 0;
		for (int i = 5; i <= 45; i += 5) {
			int lob = multiSet & -multiSet;
			partition += 1L + numberOfTrailingZeros(lob) << i;
			multiSet -= lob;
		}
		return partition;
	}
	
	static int getPartition(long partition, int index) { // multiSetでindex番目に立っているbitの位置を求める
		return (int) (partition >> 5 * index & 0b11111);
	}
	
	static long multichoose(long[] factorial, int multiSet) { // multiSetで与えられた多重集合を並べてできる組合せ
		long partition = calcPartition(multiSet);
		long multichoose = factorial[getPartition(partition, 9) - 9];
		for (int i = 0; i < 9; ++i)
			multichoose /= factorial[getPartition(partition, i + 1) - getPartition(partition, i) - 1];
		return multichoose;
	}
	
	static void nextSet(int multiSet, int[] dice, BitArray set, Stack<int> stack) {// diceを追加したときにできる新たな多重集合のうち、新しく発見したものをstackに入れる
		long partition = calcPartition(multiSet);
		foreach (int result in dice) {
			int mask = (1 << getPartition(partition, result)) - 1;
			int next = (multiSet & ~mask) << 1 | multiSet & mask;
			if (!set.Get(next)) {
				set.Set(next, true);
				stack.Push(next);
			}
		}
	}
}
0