#303603 (C++11) No.765 ukuku 2

提出ソース

結果

問題	No.765 ukuku 2
ユーザー	chocorusk
提出日時	2018-12-13 01:23:20
言語	C++11 (gcc 11.4.0)
結果	WA
実行時間	-
コード長	3,674 bytes
コンパイル時間	944 ms
コンパイル使用メモリ	103,420 KB
実行使用メモリ	14,496 KB
最終ジャッジ日時	2023-10-25 08:50:27
合計ジャッジ時間	15,854 ms
ジャッジサーバーID （参考情報）	judge15 / judge14

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	3 ms 9,708 KB
testcase_01	AC	3 ms 9,708 KB
testcase_02	AC	3 ms 9,708 KB
testcase_03	AC	3 ms 9,708 KB
testcase_04	AC	3 ms 9,708 KB
testcase_05	AC	3 ms 9,708 KB
testcase_06	AC	3 ms 9,708 KB
testcase_07	WA	-
testcase_08	AC	3 ms 9,708 KB
testcase_09	AC	3 ms 9,708 KB
testcase_10	AC	3 ms 9,708 KB
testcase_11	AC	3 ms 9,708 KB
testcase_12	AC	3 ms 9,708 KB
testcase_13	AC	3 ms 9,708 KB
testcase_14	AC	3 ms 9,708 KB
testcase_15	AC	3 ms 9,708 KB
testcase_16	AC	3 ms 9,708 KB
testcase_17	AC	3 ms 9,708 KB
testcase_18	AC	3 ms 9,712 KB
testcase_19	AC	3 ms 9,712 KB
testcase_20	AC	3 ms 9,708 KB
testcase_21	AC	3 ms 9,712 KB
testcase_22	AC	3 ms 9,712 KB
testcase_23	AC	3 ms 9,712 KB
testcase_24	AC	3 ms 9,712 KB
testcase_25	AC	3 ms 9,708 KB
testcase_26	WA	-
testcase_27	AC	3 ms 9,712 KB
testcase_28	AC	3 ms 9,712 KB
testcase_29	AC	3 ms 9,712 KB
testcase_30	AC	598 ms 14,012 KB
testcase_31	AC	553 ms 13,964 KB
testcase_32	WA	-
testcase_33	AC	609 ms 14,496 KB
testcase_34	AC	609 ms 14,496 KB
testcase_35	AC	945 ms 14,496 KB
testcase_36	AC	942 ms 14,496 KB
testcase_37	AC	589 ms 14,496 KB
testcase_38	AC	621 ms 14,496 KB
testcase_39	AC	722 ms 14,496 KB
testcase_40	AC	568 ms 14,036 KB
testcase_41	AC	664 ms 14,168 KB
testcase_42	AC	693 ms 14,284 KB
testcase_43	WA	-
testcase_44	WA	-
testcase_45	AC	710 ms 14,308 KB
testcase_46	AC	587 ms 14,072 KB
testcase_47	AC	664 ms 14,372 KB
testcase_48	AC	669 ms 14,340 KB
testcase_49	AC	660 ms 14,344 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#include <cstdio>
#include <cstring>
#include <iostream>
#include <string>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
#include <complex>
#include <unordered_map>
#include <unordered_set>
#include <random>
using namespace std;
typedef long long int ll;
typedef pair<int, int> P;
string s;
int sa[400002], lcp[400002], rk[400002];
int tmp[400002];
int n, k;
bool compare_sa(int i, int j){
  if(rk[i]!=rk[j]) return rk[i]<rk[j];
  else{
    int ri, rj;
    if(i+k<=n) ri=rk[i+k];
    else ri=-1;
    if(j+k<=n) rj=rk[j+k];
    else rj=-1;
    return ri<rj;
  }
}
void construct_sa(string s, int *sa){
  n=s.size();
  for(int i=0; i<=n; i++){
    sa[i]=i;
    if(i<n){
      rk[i]=s[i]-'a';
    }else{
      rk[i]=-1;
    }
  }
  for(k=1; k<=n; k*=2){
    sort(sa, sa+n+1, compare_sa);
    tmp[sa[0]]=0;
    for(int i=1; i<=n; i++){
      tmp[sa[i]]=tmp[sa[i-1]];
      if(compare_sa(sa[i-1], sa[i])) tmp[sa[i]]++;
    }
    for(int i=0; i<=n; i++) rk[i]=tmp[i];
  }
}
void construct_lcp(string s, int *sa, int *lcp){
	int n=s.size();
	for(int i=0; i<n; i++) rk[sa[i]]=i;
	int h=0;
	lcp[0]=0;
	for(int i=0; i<n; i++){
		int j=sa[rk[i]-1];
		if(h>0) h--;
		for(; j+h<n && i+h<n; h++){
			if(s[j+h]!=s[i+h]) break;
		}
		lcp[rk[i]-1]=h;
	}
}
const int MAX_N=1<<19;
const int INF=1e9;
int m0, mn[2*MAX_N];
void init(int m_){
	m0=1;
	while(m0<m_) m0<<=1;
	for(int i=0; i<m_; i++){
		mn[m0+i]=lcp[i];
	}
	for(int i=m_+m0; i<=2*m0-1; i++) mn[i]=INF;
	for(int i=m0-1; i>=1; i--){
		mn[i]=min(mn[2*i], mn[2*i+1]);
	}
}
int query(int a, int b){
	a+=m0, b+=m0;
	int ans=INF;
	for(;a<b; a>>=1, b>>=1){
		if(b&1) ans=min(ans, mn[--b]);
		if(a&1) ans=min(ans, mn[a++]);
	}
	return ans;
}
int main()
{
	cin>>s;
	int m=s.size();
	string t=s;
	reverse(t.begin(), t.end());
	s+='$'+t;
	construct_sa(s, sa);
	construct_lcp(s, sa, lcp);
	for(int i=0; i<=s.size(); i++) rk[sa[i]]=i;
	init(s.size()+1);
	int ans=0;
	for(int i=0; i<m; i++){
		int j=2*m-i;
		int l=query(min(rk[i], rk[j]), max(rk[i], rk[j]));
		if(i-l>=0 || i+l<=m-1) ans=max(ans, 2*l-1);
		else ans=max(ans, m-2);
		if(i-l>0 && s[i-l-1]==s[i+l]){
			int j1=2*m-(i-l-1);
			int l1=query(min(rk[j1], rk[i+l]), max(rk[i+l], rk[j1]));
			ans=max(ans, 2*l-1+2*l1);
		}
		if(i+l<m-1 && s[i-l]==s[i+l+1]){
			int j1=2*m-(i-l);
			int l1=query(min(rk[i+l+1], rk[j1]), max(rk[i+l+1], rk[j1]));
			ans=max(ans, 2*l-1+2*l1);
		}
		if(i-l>=0 && i+l-1<=m-1 && s[i-l]==s[i+l-1]){
			int j1=2*m-(i-l);
			int l1=query(min(rk[i+l-1], rk[j1]), max(rk[i+l-1], rk[j1]));
			ans=max(ans, 2*l-3+2*l1);
		}
		if(i-l+1>=0 && i+l<=m-1 && s[i-l+1]==s[i+l]){
			int j1=2*m-(i-l+1);
			int l1=query(min(rk[i+l], rk[j1]), max(rk[i+l], rk[j1]));
			ans=max(ans, 2*l-3+2*l1);
		}
	}
	for(int i=1; i<m; i++){
		int j=2*m-i+1;
		int l;
		if(s[i-1]!=s[i]){
			l=0;
		}else{
			l=query(min(rk[i], rk[j]), max(rk[i], rk[j]));
		}
		if(i-1-l>=0 || i+l<=m-1) ans=max(ans, 2*l);
		else ans=max(ans, m-2);
		if(i-l-1>0 && s[i-l-2]==s[i+l]){
			int j1=2*m-(i-l-2);
			int l1=query(min(rk[j1], rk[i+l]), max(rk[i+l], rk[j1]));
			ans=max(ans, 2*l+2*l1);
		}
		if(i+l<m-1 && s[i-l-1]==s[i+l]){
			int j1=2*m-(i-l-1);
			int l1=query(min(rk[i+l+1], rk[j1]), max(rk[i+l+1], rk[j1]));
			ans=max(ans, 2*l+2*l1);
		}
		if(i-1-l>=0 && i+l-1<=m-1 && s[i-l-1]==s[i+l-1]){
			int j1=2*m-(i-l-1);
			int l1=query(min(rk[i+l-1], rk[j1]), max(rk[i+l-1], rk[j1]));
			ans=max(ans, 2*l-2+2*l1);
		}
		if(i-l>=0 && i+l<=m-1 && s[i-l]==s[i+l]){
			int j1=2*m-(i-l-1);
			int l1=query(min(rk[i+l], rk[j1]), max(rk[i+l], rk[j1]));
			ans=max(ans, 2*l-2+2*l1);
		}
	}
	cout<<ans<<endl;
	return 0;
}

yukicoder

結果

テストケース

ソースコード