結果
| 問題 |
No.464 PPAP
|
| コンテスト | |
| ユーザー |
しらっ亭
|
| 提出日時 | 2016-11-25 18:10:29 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 13 ms / 2,000 ms |
| コード長 | 1,672 bytes |
| コンパイル時間 | 739 ms |
| コンパイル使用メモリ | 73,856 KB |
| 実行使用メモリ | 5,248 KB |
| 最終ジャッジ日時 | 2024-11-27 14:33:43 |
| 合計ジャッジ時間 | 1,731 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 4 |
| other | AC * 22 |
ソースコード
#include <iostream>
#include <vector>
#include <assert.h>
using namespace std;
struct Manachar {
vector<int> rads;
// O(|s|)
Manachar(const string &s) : rads(s.size() * 2 - 1){
int size = (int) rads.size();
int i = 0, j = 0;
while (i < size) {
while (i - j >= 0 && i + j < size && get(s, i - j) == get(s, i + j)) {
++j;
}
rads[i] = j;
int k = 1;
while (i - k >= 0 && i + k < size && k + rads[i - k] < j) {
rads[i + k] = rads[i - k], ++k;
}
i += k;
j -= k;
}
}
// s[l, r] is palindrome?
// O(1)
bool is_palindrome(int l, int r) {
assert(l >= 0);
assert(r >= l);
assert(r * 2 <= (int) rads.size());
return rads[l + r] >= r - l + 1;
}
private:
static char get(const string &s, int i) {
if (i & 1) return '^';
else return s[i >> 1];
}
};
int main() {
cin.tie(nullptr);
ios::sync_with_stdio(false);
string s;
cin >> s;
int n = (int) s.size();
assert(n >= 4);
assert(n <= 5000);
for (int i = 0; i < n; i++) {
assert('a' <= s[i] && s[i] <= 'z');
}
Manachar mana(s);
vector<int> p1;
for (int i = 0; i < n - 3; i++) {
if (mana.is_palindrome(0, i)) {
p1.push_back(i);
}
}
vector<int> p2(n - 1);
for (int i : p1) {
for (int j = i + 1; j < n - 2; j++) {
if (mana.is_palindrome(i + 1, j)) {
p2[j]++;
}
}
}
vector<long long> sum_p2(n);
for (int i = 0; i < n - 1; i++) sum_p2[i + 1] = sum_p2[i] + p2[i];
long long ans = 0;
for (int k = 2; k < n - 1; k++) {
if (mana.is_palindrome(k + 1, n - 1)) {
ans += sum_p2[k];
}
}
cout << ans << endl;
}
しらっ亭