#457888 (C++14) No.1018 suffixsuffixsuffix

• 問題: No.1018 suffixsuffixsuffix
• ユーザー: HIR180
• 提出日時: 2020-04-03 22:32:24
• 言語: C++14 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 815 ms / 2,000 ms
• コード長: 3,499 bytes
• コンパイル時間: 2,025 ms
• コンパイル使用メモリ: 201,264 KB
• 実行使用メモリ: 18,384 KB
• 最終ジャッジ日時: 2023-09-16 03:21:38
• 合計ジャッジ時間: 15,282 ms
• ジャッジサーバーID （参考情報）: judge11 / judge14

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 7,548 KB
testcase_01	AC	3 ms; 8,244 KB
testcase_02	AC	3 ms; 8,224 KB
testcase_03	AC	3 ms; 8,236 KB
testcase_04	AC	2 ms; 7,532 KB
testcase_05	AC	2 ms; 7,540 KB
testcase_06	AC	2 ms; 7,548 KB
testcase_07	AC	2 ms; 7,644 KB
testcase_08	AC	2 ms; 7,544 KB
testcase_09	AC	228 ms; 8,380 KB
testcase_10	AC	234 ms; 8,364 KB
testcase_11	AC	247 ms; 8,528 KB
testcase_12	AC	228 ms; 8,300 KB
testcase_13	AC	229 ms; 8,432 KB
testcase_14	AC	368 ms; 9,128 KB
testcase_15	AC	552 ms; 9,804 KB
testcase_16	AC	764 ms; 10,548 KB
testcase_17	AC	537 ms; 9,632 KB
testcase_18	AC	380 ms; 8,980 KB
testcase_19	AC	131 ms; 8,568 KB
testcase_20	AC	132 ms; 8,300 KB
testcase_21	AC	134 ms; 8,276 KB
testcase_22	AC	138 ms; 8,464 KB
testcase_23	AC	138 ms; 8,220 KB
testcase_24	AC	771 ms; 17,352 KB
testcase_25	AC	815 ms; 18,384 KB
testcase_26	AC	734 ms; 17,348 KB
testcase_27	AC	762 ms; 17,288 KB
testcase_28	AC	751 ms; 17,604 KB
testcase_29	AC	144 ms; 8,484 KB
testcase_30	AC	147 ms; 8,520 KB
testcase_31	AC	145 ms; 8,440 KB
testcase_32	AC	145 ms; 8,660 KB
testcase_33	AC	138 ms; 8,436 KB
testcase_34	AC	2 ms; 7,532 KB
testcase_35	AC	142 ms; 8,320 KB
testcase_36	AC	139 ms; 8,420 KB
testcase_37	AC	784 ms; 18,336 KB

コンパイルメッセージ

main.cpp: 関数 ‘std::string reduce(std::string)’ 内:
main.cpp:141:1: 警告: 制御が非 void 関数の終りに到達しました [-Wreturn-type]
  141 | }
      | ^

ソースコード

//Let's join Kaede Takagaki Fan Club !!
#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
using namespace std;
typedef long long ll;
typedef pair<int,int> P;
typedef pair<int,P> P1;
typedef pair<P,P> P2;
#define pu push
#define pb push_back
#define mp make_pair
#define eps 1e-7
#define INF 1000000000
#define fi first
#define sc second
#define rep(i,x) for(int i=0;i<x;i++)
#define repn(i,x) for(int i=1;i<=x;i++)
#define SORT(x) sort(x.begin(),x.end())
#define ERASE(x) x.erase(unique(x.begin(),x.end()),x.end())
#define POSL(x,v) (lower_bound(x.begin(),x.end(),v)-x.begin())
#define POSU(x,v) (upper_bound(x.begin(),x.end(),v)-x.begin())
#define all(x) x.begin(),x.end()
template<class T>
void dmp(T a){
	rep(i,a.size()) cout << a[i] << " ";
	cout << endl;
}
template<class T>
bool chmax(T&a, T b){
	if(a < b){
		a = b;
		return 1;
	}
	return 0;
}
template<class T>
bool chmin(T&a, T b){
	if(a > b){
		a = b;
		return 1;
	}
	return 0;
}
template<class T>
void g(T &a){
	cin >> a;
}
template<class T>
void o(const T &a,bool space=false){
	cout << a << (space?' ':'\n');
}
//ios::sync_with_stdio(false);
const ll mod = 1000000007;//998244353
template<class T>
void add(T&a,T b){
	a+=b;
	if(a >= mod) a-=mod;
}


int NN,k;
int ran[400005];
int tmp[400005];
int sa[400005];

bool compare_sa(int i,int j)
{
	if(ran[i] != ran[j]) return ran[i] < ran[j];
	else
	{
		int ri = i+k<=NN ? ran[i+k]: -1;
		int rj = j+k<=NN ? ran[j+k]: -1;
		
		return ri < rj;
	}
}

void construct_sa(string S)
{
	NN = S.size();
	for(int i=0;i<=NN;i++)
	{
		sa[i] = i;
		ran[i] = i<NN?S[i]:-1;
	}
	
	for(k=1;k<=NN;k*=2)
	{
		sort(sa,sa+NN+1,compare_sa);
		
		tmp[sa[0]] = 0;
		for(int i=1;i<=NN;i++)
		{
			tmp[sa[i]] = tmp[sa[i-1]] + compare_sa(sa[i-1],sa[i]);
		}
		for(int i=0;i<=NN;i++)
		{
			ran[i] = tmp[i];
		}
	}
}
int lcp[400005];
void construct_lcp(string S)
{
	int n = S.size();
	for(int i=0;i<=n;i++) ran[sa[i]] = i;
	
	int h = 0;
	lcp[0] = 0;
	
	for(int i=0;i<n;i++)
	{
		int j = sa[ran[i]-1];
		
		if(h) h--;
		for(;j+h<n && i+h<n;h++)
		{
			if(S[j+h] != S[i+h]) break;
		}
		lcp[ran[i]-1] = h;
	}
}

string s;
ll n,m,q;
ll ans[100005];
ll ans2[100005];

string reduce(string str){
	int sz = str.size();
	for(int i=1;i<=sz;i++){
		if(sz%i != 0) continue;
		string ret = str.substr(0, i);
		for(int j=i;j<sz;j++){
			if(str[j-i] != str[j]) goto fail;
		}
		return ret;
		fail:;
	}
}
int main(){
	ios::sync_with_stdio(false);
	cin >> n >> m >> q;
	cin >> s;
	s = reduce(s);
	m *= n/s.size();
	n = s.size();
	//cout << n << " " << m << endl;
	if(m > 4){
		int n = s.size();
		s = s+s+s+s;
		construct_sa(s);
		construct_lcp(s);
		ll cur = 0;
		rep(i,100005) ans2[i] = -1;
		
		map<ll,int>M;
		vector<pair<ll,int>>vec;
		for(int i=1;i<=4*n;i++){
		    //cout << sa[i]%n << " " << (sa[i] < 3*n?1:0) << endl;
			if(sa[i] >= 3*n){
				M[++cur] = sa[i]%n;
			}
			else if(ans2[sa[i]%n] == -1){
				//answer for n*(m-2)+sa[i]%n
				ans2[sa[i]%n] = 1;
				vec.pb(mp(++cur, sa[i]%n));
				cur += (m-2);
			}
		}
		rep(i,q) {
			ll a; g(a);
			if(M.find(a) != M.end()){
				o(1LL*(m-1)*n + M[a] + 1, (i==q-1?0:1));
			}
			else{
				int x = POSL(vec, mp(a, INF));
				x--;
				assert(vec[x].fi <= a && a <= vec[x].fi+m-2);
				int id = vec[x].sc;
				o(1LL*(m-2)*n+id - 1LL*n*(a-vec[x].fi) + 1, (i==q-1?0:1));
			}
		}
	}
	else{
		string t = "";
		rep(i,m) t += s;
		construct_sa(t);
		construct_lcp(t);
		rep(i,q) {
			int a; g(a); o(1+sa[a], (i==q-1?0:1)); 
		}
	}
}