kmjp's blog

競技プログラミング参加記です

Codeforces ECR #113 : F. Palindromic Hamiltonian Path

コンテスト中のAC数がだいぶ少な目。
https://codeforces.com/contest/1569/problem/F

問題

N頂点(偶数)のグラフが与えられる。
各頂点に、アルファベット先頭K文字のいずれかを割り当てるとき、ハミルトンパスのうち訪問した頂点に割り当てた文字を並べると回文になるものが存在するものは何通りか。

解法

愚直に文字の割り当てを考えるとO(N^K)通りかかるので枝刈りする。
まず、各文字の登場回数は偶数でなければならない。
また、全文字を総当たりするのは無駄で、どの頂点同士の文字が等しい・異なるという条件だけわかれば、具体的な文字の列挙は最後に組み合わせ計算で行えばよい。

上記を踏まえて文字種を絞って枝刈りしよう。

int N,M,K;
int E[12][12];

static const int N_=1020;
static ll P_[N_][N_];
int V[12];
ll ret;

map<vector<pair<int,int>>,int> memo;
int valid(vector<pair<int,int>> P) {
	sort(ALL(P));
	if(memo.count(P)) return memo[P];
	int ret=0;
	do {
		int i;
		int ok=1;
		if(E[P[0].first][P[0].second]==0) ok=0;
		for(i=1;i<P.size();i++) {
			int can=(E[P[i].first][P[i-1].first]&&E[P[i].second][P[i-1].second])|(E[P[i].first][P[i-1].second]&&E[P[i].second][P[i-1].first]);
			ok&=can;
		}
		if(ok) return memo[P]=1;
	} while(next_permutation(ALL(P)));
	return memo[P]=0;
}


int ok(int mask,vector<pair<int,int>>& P) {
	if(mask==0) return valid(P);
	int i,j;
	FOR(i,N) if(mask&(1<<i)) break;
	mask^=1<<i;
	FOR(j,N) if((mask&(1<<j))&&V[i]==V[j]) {
		P.push_back({i,j});
		int v=ok(mask^(1<<j),P);
		P.pop_back();
		if(v) return 1;
	}
	return 0;
	
}


void dfs(int cur,int ma,int mask) {
	int i;
	if(cur==N) {
		if(mask==0) {
			vector<pair<int,int>> P;
			if(ok((1<<N)-1,P)) ret+=P_[K][ma];
			
		}
		return;
	}
	
	for(i=0;i<min(ma+1,N/2);i++) {
		V[cur]=i;
		dfs(cur+1,max(ma,i+1),mask^(1<<i));
	}
}


void solve() {
	int i,j,k,l,r,x,y; string s;
	
	FOR(i,15) {
		P_[i][0]=1;
		for(j=1;j<=i;j++) P_[i][j]=P_[i][j-1]*(i+1-j);
	}
	
	
	cin>>N>>M>>K;
	FOR(i,M) {
		cin>>x>>y;
		E[x-1][y-1]=E[y-1][x-1]=1;
	}
	dfs(0,0,0);
	cout<<ret<<endl;
}

まとめ

コード量がまぁまぁ多いな。