結果

問題 No.464 PPAP
ユーザー goodbaton
提出日時 2019-02-19 22:43:17
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 198 ms / 2,000 ms
コード長 2,151 bytes
コンパイル時間 939 ms
コンパイル使用メモリ 101,516 KB
実行使用メモリ 73,088 KB
最終ジャッジ日時 2024-10-12 23:07:08
合計ジャッジ時間 2,397 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 22
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>

#include <iostream>
#include <complex>
#include <string>
#include <algorithm>
#include <numeric>
#include <vector>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <unordered_map>
#include <unordered_set>

#include <functional>
#include <cassert>

typedef long long ll;
using namespace std;

#ifndef LOCAL
#define debug(x) ;
#else
#define debug(x) cerr << __LINE__ << " : " << #x << " = " << (x) << endl;

template <typename T1, typename T2>
ostream &operator<<(ostream &out, const pair<T1, T2> &p) {
  out << "{" << p.first << ", " << p.second << "}";
  return out;
}

template <typename T>
ostream &operator<<(ostream &out, const vector<T> &v) {
  out << '{';
  for (const T &item : v) out << item << ", ";
  out << "\b\b}";
  return out;
}
#endif

#define mod 1000000007 //1e9+7(prime number)
#define INF 1000000000 //1e9
#define LLINF 2000000000000000000LL //2e18
#define SIZE 200010

/* Manacher.cpp (Source: snukeさんブログ) */

// R[i] := 文字 i を中心とする最長の回文の半径( 全長+1 / 2 )
// 偶数長は文字間にダミーを入れる(例: a#b#a#c#b)
void Manacher(const char *S, int *R){
  int i = 0, j = 0, k, n = strlen(S);
  while (S[i]) {
    while (i-j >= 0 && i+j < n && S[i-j] == S[i+j]) j++;
    R[i] = j;
    for (k = 1; i-k >= 0 && i+k < n && k+R[i-k] < j; k++) R[i+k] = R[i-k];
    i += k; j -= k;
  }
}

char s[SIZE], t[SIZE];
vector<int> vec[SIZE];
ll dp[5][SIZE];

int main(){
  int n;

  scanf("%s", s);
  n = strlen(s);

  for(int i=0;i<n;i++) {
    t[i*2] = '#';
    t[i*2+1] = s[i];
  }
  t[n*2] = '#';

  int pal[SIZE];

  Manacher(t, pal);

  for(int i=0;i<n;i++)
    for(int j=0;j<pal[i*2+1]/2;j++)
      vec[i-j].push_back(i+j);

  for(int i=0;i<n-1;i++)
    for(int j=0;j<pal[i*2+2]/2;j++)
      vec[i-j].push_back(i+j+1);

  dp[0][0] = 1;

  ll sum = 0;

  for(int i=0;i<n;i++){
    dp[3][i] += sum;
    for(int to : vec[i]) {
      dp[1][to+1] += dp[0][i];
      dp[2][to+1] += dp[1][i];
      dp[4][to+1] += dp[3][i];
    }
    sum += dp[2][i];
  }

  cout << dp[4][n] << endl;

  return 0;
}
0