Codeforces Global Round 26 : D. "a" String Problem

かなりいまいちだった出来の回。
https://codeforces.com/contest/1984/problem/D

問題

文字列Sが与えられる。
以下を満たす、Tは何通りあるか。

• Tは"a"以外ではない。
• SをいくつかのTと"a"の連結で表現できる。ただしTが1個以上なければならない。

解法

Sが'a'だけで構成される場合は、コーナーケースとして除いていく。
'b'～'z'の登場頻度から、TがS中に現れる回数は、それらのGCDとして表現できる。
その回数sを総当たりしよう。
'b'～'z'の登場するindexを並べてs分割し、それぞれを含む部分文字列のhashを取ることで、s個が一致するか判定するとよい。

int T,N;
string S;
vector<int> C[26];

using VT = string;

struct RollingHash {
	static const ll mo0=1000000021,mo1=1000000009;
	static ll mul0,mul1;
	static const ll add0=1000010007, add1=1003333331;
	static vector<ll> pmo[2];
	VT s; int l; vector<ll> hash_[2];
	void init(VT s) {
		this->s=s; l=s.size(); int i,j;
		hash_[0]=hash_[1]=vector<ll>(1,0);
		if(!mul0) mul0=10009+(((ll)&mul0+time(NULL))>>5)%1259,mul1=10007+(time(NULL)+((ll)&mul1)>>5)%2257;
		if(pmo[0].empty()) pmo[0].push_back(1),pmo[1].push_back(1);
		FOR(i,l) hash_[0].push_back((hash_[0].back()*mul0+add0+s[i])%mo0);
		FOR(i,l) hash_[1].push_back((hash_[1].back()*mul1+add1+s[i])%mo1);
	}
	pair<ll,ll> hash(int l,int r) { // s[l..r]
		if(l>r) return make_pair(0,0);
		while(pmo[0].size()<r+2)
			pmo[0].push_back(pmo[0].back()*mul0%mo0), pmo[1].push_back(pmo[1].back()*mul1%mo1);
		return make_pair((hash_[0][r+1]+(mo0-hash_[0][l]*pmo[0][r+1-l]%mo0))%mo0,
			             (hash_[1][r+1]+(mo1-hash_[1][l]*pmo[1][r+1-l]%mo1))%mo1);
	}
	pair<ll,ll> hash(VT s) { init(s); return hash(0,s.size()-1); }
	static pair<ll,ll> concat(pair<ll,ll> L,pair<ll,ll> R,int RL) { // hash(L+R) RL=len-of-R
		while(pmo[0].size()<RL+2) pmo[0].push_back(pmo[0].back()*mul0%mo0), pmo[1].push_back(pmo[1].back()*mul1%mo1);
		return make_pair((R.first + L.first*pmo[0][RL])%mo0,(R.second + L.second*pmo[1][RL])%mo1);
	}
};
vector<ll> RollingHash::pmo[2]; ll RollingHash::mul0,RollingHash::mul1;
RollingHash rh;

void solve() {
	int i,j,k,l,r,x,y; string s;
	
	cin>>T;
	while(T--) {

		cin>>S;
		rh.init(S);
		N=S.size();
		FOR(i,26) C[i].clear();
		FOR(i,N) C[S[i]-'a'].push_back(i);
		
		if(C[0].size()==N) {
			cout<<N-1<<endl;
			continue;
		}
		int g=0;
		for(i=1;i<26;i++) g=__gcd(g,(int)C[i].size());
		ll ret=0;
		for(int s=1;s<=g;s++) if(g%s==0) {
			//s個に分割
			vector<pair<int,int>> cand;
			FOR(i,s) cand.push_back({N,0});
			int ok=1;
			for(i=1;i<26;i++) if(C[i].size()) {
				int span=C[i].size()/s;
				FOR(j,C[i].size()) {
					chmin(cand[j/span].first,C[i][j]);
					chmax(cand[j/span].second,C[i][j]);
				}
			}
			auto h1=rh.hash(cand[0].first,cand[0].second);
			FOR(i,s) {
				auto h2=rh.hash(cand[i].first,cand[i].second);
				if(h1!=h2) ok=0;
			}
			int TL=cand[0].first;
			int TR=N-1-cand.back().second;
			int mi=1<<20;
			FOR(i,s-1) {
				mi=min(mi,cand[i+1].first-cand[i].second-1);
				if(cand[i].second>=cand[i+1].first) ok=0;
			}
			if(ok==0) continue;
			if(TL>TR) swap(TL,TR);
			for(i=0;i<=min(TL,mi);i++) ret+=min(TR,mi-i)+1;
		}
		cout<<ret<<endl;
		
	}
}

まとめ

考え方は難しくないけど、ちょっとコード量多かったな。