結果

問題 No.2606 Mirror Relay
ユーザー 沙耶花沙耶花
提出日時 2024-01-12 23:15:33
言語 C++17(gcc12)
(gcc 12.3.0 + boost 1.87.0)
結果
AC  
実行時間 497 ms / 2,000 ms
コード長 3,678 bytes
コンパイル時間 5,268 ms
コンパイル使用メモリ 280,292 KB
実行使用メモリ 32,004 KB
最終ジャッジ日時 2024-09-27 23:58:06
合計ジャッジ時間 13,089 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 69
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <stdio.h>
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace atcoder;
using mint = modint998244353;
using namespace std;
#define rep(i,n) for (int i = 0; i < (n); ++i)
#define Inf32 1000000001
#define Inf64 1000000000000000001

int op(int a,int b){
	return min(a,b);
}
int e(){
	return Inf32;
}
int X = 0;
bool f(int x){
	return x >= X;
}
struct rolling_hash{

	long long t_hash;
	static inline vector<long long> power;
	static const long long MOD = (1LL<<61)-1;
	static const long long b = 123456;
	int sz;
	
	rolling_hash(){		
		sz = 0;
		t_hash = 0;
	}
	
	rolling_hash(char c){
		sz = 1;
		t_hash = b*c;
	}

	long long mul(__int128 x,__int128 y){
		__int128 t = x*y;
		t = (t>>61) + (t&MOD);
		
		if(t>=MOD)t -= MOD;
		return t;
	}
	
	long long get_pow(int sz){
		if(power.size()>sz)return power[sz];
		
		while(power.size()<=sz){
			if(power.size()==0)power.push_back(1);
			else power.push_back(mul(power.back(),b));
		}
		return power.back();
		
	}
	
	rolling_hash &operator+=(const rolling_hash &another){
		
		(*this).t_hash = mul((*this).t_hash,get_pow(another.sz));
		(*this).t_hash += another.t_hash;
		if((*this).t_hash>=MOD)(*this).t_hash -= MOD;
			
		(*this).sz += another.sz;
		
		return (*this);
	}
	
	rolling_hash operator+(const rolling_hash &another)const{
		return (rolling_hash(*this)+=another);
	}
	
	rolling_hash &operator-=(const rolling_hash &another){

		(*this).t_hash += MOD - mul(another.t_hash,get_pow((*this).sz-another.sz));
		if((*this).t_hash>=MOD)(*this).t_hash -= MOD;
			
		(*this).sz -= another.sz;

		return (*this);
	}
	
	rolling_hash operator-(const rolling_hash &another)const{
		return (rolling_hash(*this)-=another);
	}
	
	bool operator<(const rolling_hash &another)const{
		if((*this).t_hash!=another.t_hash)return (*this).t_hash<another.t_hash;
		return (*this).sz<another.sz;
	}
	
	bool operator==(const rolling_hash &another)const{
		return ((*this).t_hash==another.t_hash && (*this).sz==another.sz);
	}

	
};

int main(){
	
	string s;
	cin>>s;
	int n = s.size();
	vector<rolling_hash> R(n+1),rR(n+1);
	rep(i,n){
		R[i+1] = R[i] + rolling_hash(s[i]);
	}
	reverse(s.begin(),s.end());
	rep(i,n){
		rR[i+1] = rR[i] + rolling_hash(s[i]);
	}
	
	reverse(s.begin(),s.end());
	set<rolling_hash> ss;
	vector<pair<int,int>> p;
	rep(i,n){
		int ok = 0,ng = n;
		while(ng-ok>1){
			int mid = (ok+ng)/2;
			int l = i-mid,r = i+1+mid;
			if(l<0||r>n){
				ng = mid;
				continue;
			}
			if((R[r]-R[l]) == (rR[n-l] - rR[n-r]))ok = mid;
			else ng = mid;
		}
		int l = i-ok,r = i+1+ok;
		while(l<r){
			if(ss.count(R[r]-R[l]))break;
			ss.insert(R[r]-R[l]);
			p.emplace_back(l,r);
			l++;
			r--;
		}
	}
	
	rep(i,n-1){
		if(s[i]!=s[i+1])continue;
		int ok = 0,ng = n;
		while(ng-ok>1){
			int mid = (ok+ng)/2;
			int l = i-mid,r = i+2+mid;
			if(l<0||r>n){
				ng = mid;
				continue;
			}
			if((R[r]-R[l]) == (rR[n-l] - rR[n-r]))ok = mid;
			else ng = mid;
		}
		int l = i-ok,r = i+2+ok;
		while(l<r){
			if(ss.count(R[r]-R[l]))break;
			ss.insert(R[r]-R[l]);
			p.emplace_back(l,r);
			l++;
			r--;
		}
	}
	sort(p.begin(),p.end());
	auto sa = suffix_array(s);
	auto la = lcp_array(s,sa);
	segtree<int,op,e> seg(la);
	vector<long long> imos(n+1);
	rep(i,sa.size()){
		int d = distance(p.begin(),lower_bound(p.begin(),p.end(),make_pair(sa[i],-1)));
		for(int j=d;j<p.size();j++){
			if(p[j].first!=sa[i])break;
			int len = p[j].second - p[j].first;
			X = len;
			int ll = seg.min_left<f>(i);
			int rr = seg.max_right<f>(i);
			rr++;
			long long t = rr-ll;
			t *= len;
			imos[ll] += t;
			imos[rr] -= t;
			
		}
	}
	rep(i,n)imos[i+1] += imos[i];
	cout<<(*max_element(imos.begin(),imos.end()))<<endl;
	
	return 0;
}
0