結果
問題 | No.464 PPAP |
ユーザー |
|
提出日時 | 2016-12-15 23:44:12 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 157 ms / 2,000 ms |
コード長 | 1,306 bytes |
コンパイル時間 | 1,429 ms |
コンパイル使用メモリ | 169,440 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-11-30 09:00:44 |
合計ジャッジ時間 | 3,054 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 22 |
ソースコード
#include <bits/stdc++.h> using namespace std; int main() { string s; cin >> s; int n = s.size(); const int mod = 1e9 + 7; srand(time(NULL)); int64_t base = rand() + 2; vector<int64_t> pw(5050); pw[0] = 1; for (int i = 1; i < 5050; i++) { pw[i] = base * pw[i - 1] % mod; } vector<int64_t> P(n + 2); for (int i = 0; i < n; i++) { P[i + 1] = 1; for (int j = 0; j <= i; j++) { if (s[j] != s[i - j]) { P[i + 1] = 0; } } } vector<int64_t> PP(n + 2); for (int i = n - 1; i >= 0; i--) { int64_t a = 0; int64_t b = 0; for (int j = 0; i + j < n; j++) { a = a * base + s[i + j]; b = b + s[i + j] * pw[j]; a %= mod; b %= mod; if (a == b) { PP[i + j + 2] += P[i]; } } } for (int i = 0; i < n; i++) { PP[i + 1] += PP[i]; } int64_t ans = 0; for (int i = n - 1; i >= 0; i--) { bool pal = true; for (int j = 0; j <= n - 1 - i; j++) { if (s[i + j] != s[n - 1 - j]) { pal = false; } } if (pal) { ans += PP[i]; } } cout << ans << endl; }