結果

問題 No.852 連続部分文字列
ユーザー MarcusAureliusAntoninusMarcusAureliusAntoninus
提出日時 2019-07-26 21:42:57
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 307 ms / 3,153 ms
コード長 1,591 bytes
コンパイル時間 2,102 ms
コンパイル使用メモリ 198,284 KB
最終ジャッジ日時 2025-01-07 07:47:50
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 41
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

///////////////////////
// 部分列オートマトン //
//////////////////////

// 文字列の連続とは限らない部分列を重複なく列挙するオートマトンを行列の形式で作成する
// ret[i][j]=kはi番目の文字のあとに文字'a'+jがくる場合kに遷移することを意味する
// ただしインデックスは引数の文字列の前後に番兵を置いた文字列S'と一致している
// 例えばS="abcbca"ならば
// S' = $ a b c b c a $
// ret= 1 6 6 6 6 6 7 7
//      2 2 4 4 7 7 7 7
//      3 3 3 5 5 7 7 7
// となる(ここではabcの3文字のみで考えたが、実際は26文字のものが作成される)

std::vector<std::array<int, 26>> subsequenceAutomaton(const std::string& S, const char initial_char = 'a')
{
	std::vector<std::array<int, 26>> table(S.size() + 2);
	std::fill(table.back().begin(), table.back().end(), (int)table.size() - 1);
	table[table.size() - 2] = table.back();

	for (int s_i{(int)S.size() - 1}; s_i >= 0; s_i--)
	{
		std::copy(table[s_i + 1].begin(), table[s_i + 1].end(), table[s_i].begin());
		table[s_i][S[s_i] - initial_char] = s_i + 1;
	}
	return std::move(table);
}


int main()
{
	std::string S;
	std::cin >> S;
	auto table{subsequenceAutomaton(S)};
	int64_t count{};
	for (int c_i{}; c_i < 26; c_i++)
		for (int s_i{}; s_i < (int)S.size(); s_i++)
			count += (int)S.size() - table[s_i][c_i] + 1;
	long double ans{(long double)count};
	ans /= ((int64_t)S.size() + 1) * S.size() / 2;
	std::cout << std::setprecision(10) << std::fixed << ans << std::endl;

	return 0;
}
0