結果

問題 No.562 超高速一人かるた small
ユーザー fine
提出日時 2017-08-26 01:16:30
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 523 ms / 3,000 ms
コード長 1,522 bytes
コンパイル時間 2,287 ms
コンパイル使用メモリ 179,644 KB
実行使用メモリ 183,784 KB
最終ジャッジ日時 2024-10-15 17:00:51
合計ジャッジ時間 9,116 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 21
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

using ll = long long;

const ll MOD = 1e9 + 7;

ll cost[20][20], fact[20];
ll dp1[1 << 20][20], dp2[1 << 20], cnt[1 << 20], ans[21];

int main() {
	cin.tie(0);
	ios::sync_with_stdio(false);
	int n;
	cin >> n;
	vector<string> s(n);
	for (int i = 0; i < n; i++) {
		cin >> s[i];
	}
	sort(s.begin(), s.end());

	for (int i = 0; i < n; i++) {
		for (int j = 0; j < i; j++) {
			cost[i][j] = 1;
			for (int k = 0; k < min(s[i].length(), s[j].length()); k++) {
				if (s[i][k] != s[j][k]) break;
				cost[i][j]++;
			}
			cost[j][i] = cost[i][j];
		}
	}

	int hoge = (1 << n) - 1;

	for (int i = 0; i < n; i++) dp1[1 << i][i] = 1;
	for (int i = 1; i <= hoge; i++) {
		vector<int> tmp;
		for (int j = 0; j < n; j++) {
			if (i & (1 << j)) {
				tmp.push_back(j);
			}
		}

		for (int j = 0; j < tmp.size(); j++) {
			if (j > 0) {
				dp1[i][tmp[j]] = max(dp1[i][tmp[j]], cost[tmp[j - 1]][tmp[j]]);
			}
			if (j + 1 < tmp.size()) {
				dp1[i][tmp[j]] = max(dp1[i][tmp[j]], cost[tmp[j + 1]][tmp[j]]);
			}
		}
	}

	fact[0] = 1;
	for (int i = 1; i < n; i++) fact[i] = fact[i - 1] * i % MOD;;

	for (int i = 0; i <= hoge; i++) cnt[i] = __builtin_popcount(i);

	for (int i = hoge; i > 0; i--) {
		for (int j = 0; j < n; j++) {
			if (i & (1 << j)) {
				(dp2[i ^ (1 << j)] += dp2[i] + fact[n - cnt[i]] * dp1[i][j]) %= MOD;
			}
		}
	}

	for (int i = 0; i < hoge; i++) {
		(ans[n - cnt[i]] += dp2[i]) %= MOD;
	}

	for (int i = 1; i <= n; i++) {
		cout << ans[i] << endl;
	}
	return 0;
}
0