#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
実行時間	866 ms / 2,000 ms
コード長	3,499 bytes
コンパイル時間	3,627 ms
コンパイル使用メモリ	204,220 KB
実行使用メモリ	18,964 KB
最終ジャッジ日時	2024-07-03 04:40:42
合計ジャッジ時間	16,709 ms
ジャッジサーバーID （参考情報）	judge4 / judge1

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

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	3 ms 9,572 KB
testcase_01	AC	4 ms 10,464 KB
testcase_02	AC	4 ms 9,572 KB
testcase_03	AC	4 ms 8,300 KB
testcase_04	AC	3 ms 9,572 KB
testcase_05	AC	4 ms 9,572 KB
testcase_06	AC	3 ms 9,696 KB
testcase_07	AC	3 ms 9,568 KB
testcase_08	AC	3 ms 7,780 KB
testcase_09	AC	236 ms 10,008 KB
testcase_10	AC	244 ms 10,016 KB
testcase_11	AC	253 ms 9,896 KB
testcase_12	AC	238 ms 9,700 KB
testcase_13	AC	243 ms 8,724 KB
testcase_14	AC	392 ms 8,800 KB
testcase_15	AC	555 ms 10,176 KB
testcase_16	AC	804 ms 10,808 KB
testcase_17	AC	564 ms 9,540 KB
testcase_18	AC	387 ms 9,812 KB
testcase_19	AC	133 ms 8,424 KB
testcase_20	AC	138 ms 8,428 KB
testcase_21	AC	137 ms 9,828 KB
testcase_22	AC	138 ms 10,088 KB
testcase_23	AC	141 ms 9,956 KB
testcase_24	AC	800 ms 18,404 KB
testcase_25	AC	866 ms 18,964 KB
testcase_26	AC	747 ms 17,308 KB
testcase_27	AC	797 ms 17,600 KB
testcase_28	AC	778 ms 17,612 KB
testcase_29	AC	149 ms 9,916 KB
testcase_30	AC	150 ms 8,424 KB
testcase_31	AC	146 ms 8,424 KB
testcase_32	AC	148 ms 9,700 KB
testcase_33	AC	144 ms 9,832 KB
testcase_34	AC	3 ms 9,568 KB
testcase_35	AC	144 ms 10,212 KB
testcase_36	AC	144 ms 8,520 KB
testcase_37	AC	835 ms 18,892 KB

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

コンパイルメッセージ

main.cpp: In function 'std::string reduce(std::string)':
main.cpp:141:1: warning: control reaches end of non-void function [-Wreturn-type]
  141 | }
      | ^

ソースコード

raw source code

//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)); 
		}
	}
}

yukicoder

結果

テストケース

コンパイルメッセージ

ソースコード