結果

問題 No.3285 Chorus with Friends
ユーザー highlighter
提出日時 2025-09-26 19:10:50
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 219 ms / 3,000 ms
コード長 1,337 bytes
コンパイル時間 6,852 ms
コンパイル使用メモリ 296,116 KB
実行使用メモリ 81,920 KB
最終ジャッジ日時 2025-09-26 19:11:01
合計ジャッジ時間 7,287 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 40
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
#include<atcoder/modint>
using namespace std;
using mint=atcoder::modint998244353;

int main(){
	int N,M;
	cin >> N >> M;
	vector<vector<int>> A(N,vector<int>(M+1,0));
	for(int i=0;i<N;i++) for(int j=0;j<M;j++) cin >> A[i][j];
	if(N<=M){
		vector<vector<mint>> dp(M+1,vector<mint>(1<<N));
		int val[301];
		memset(val,0,sizeof(val));
		for(int i=0;i<N;i++) val[A[i][0]]|=1<<i;
		for(int i=1;i<=300;i++) if(val[i]>0) dp[0][val[i]]++;
		for(int i=1;i<=M;i++){
			for(unsigned int j=1;j<1<<N;j++){
				int t=A[countr_zero(j)][0];
				for(int k=0;k<N;k++){
					if(j&(1<<k)){
						if(A[k][i]==t) dp[i][j]+=dp[i-1][j];
						else dp[i][j^(1<<k)]+=dp[i-1][j];
					}
				}
			}
		}
		cout << reduce(dp[M].begin(),dp[M].end()).val() << '\n';
	}
	else{
		mint ans;
		sort(A.begin(),A.end());
		mint dp[N][M+1][1<<M];
		for(int i=0;i<N;i++){
			if(i==0 || A[i][0]!=A[i-1][0]) dp[i][0][0]=1;
			else for(int j=0;j<=M;j++) for(int k=0;k<1<<M;k++) dp[i][0][k]+=dp[i-1][j][k];
			mint sum[1<<M];
			for(int j=0;j<1<<M;j++) sum[j]=dp[i][0][j];
			for(int j=1;j<=M;j++){
				for(int k=0;k<1<<M;k++){
					if(k&(1<<j-1)) dp[i][j][k]=sum[k^(1<<j-1)];
					if(A[i][j]==A[i][0]) sum[k]+=dp[i][j][k];
				}
			}
			if(i==N-1 || A[i][0]!=A[i+1][0]) for(int j=0;j<=M;j++) ans+=dp[i][j][(1<<M)-1];
		}
		cout << ans.val() << '\n';
	}
}
0