結果

問題 No.3188 K-th Lexmin
コンテスト
ユーザー kmjp
提出日時 2025-06-20 23:54:55
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 559 ms / 3,500 ms
コード長 2,838 bytes
コンパイル時間 2,571 ms
コンパイル使用メモリ 214,168 KB
実行使用メモリ 19,880 KB
最終ジャッジ日時 2025-06-20 23:55:16
合計ジャッジ時間 20,085 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 47
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
typedef signed long long ll;

#define _P(...) (void)printf(__VA_ARGS__)
#define FOR(x,to) for(x=0;x<(to);x++)
#define FORR(x,arr) for(auto& x:arr)
#define FORR2(x,y,arr) for(auto& [x,y]:arr)
#define ALL(a) (a.begin()),(a.end())
#define ZERO(a) memset(a,0,sizeof(a))
#define MINUS(a) memset(a,0xff,sizeof(a))
template<class T> bool chmax(T &a, const T &b) { if(a<b){a=b;return 1;}return 0;}
template<class T> bool chmin(T &a, const T &b) { if(a>b){a=b;return 1;}return 0;}
//-------------------------------------------------------

int T;
int N;
ll M,K;
vector<int> S[201010],G;
int TL[201010];
ll LS[201010];
pair<int,int> re[201010];

ll X[202020];
int ret[202020][3];
int nex;
ll did;

template<typename ST=vector<int>>
struct SuffixArray {
	int N; vector<int> rank,lcp,sa,rsa; ST S;
 
	void build(ST S){
		this->S=S;
		int i,h=0; vector<int> tmp;
		N=S.size(); rank.resize(N+1); sa.resize(N+1); tmp.resize(N+1);
		FOR(i,N+1) sa[i]=i, rank[i]=i==N?-1:S[i];
		
		for(int k=1; k<=N; k<<=1) {
			auto pred2 = [k,this](int& a,int &b)->bool{ return (((a+k<=N)?rank[a+k]:-1)<((b+k<=N)?rank[b+k]:-1));};
			auto pred = [pred2,k,this](int& a,int &b)->bool{ return (rank[a]!=rank[b])?(rank[a]<rank[b]):pred2(a,b);};
			int x=0;
			if(k!=1) for(i=1;i<N+1;i++) if(rank[sa[i]]!=rank[sa[x]]) sort(sa.begin()+x,sa.begin()+i,pred2), x=i;
			sort(sa.begin()+x,sa.end(),pred);
			FOR(i,N+1) tmp[sa[i]]=(i==0)?0:tmp[sa[i-1]]+pred(sa[i-1],sa[i]);
			swap(rank,tmp);
		}
		lcp.resize(N+1); rsa.resize(N+1);
		FOR(i,N+1) rsa[sa[i]]=i;
		FOR(i,N) {
			int j=sa[rsa[i]-1];
			for(h=max(h-1,0);i+h<N && j+h<N; h++) if(S[j+h]!=S[i+h]) break;
			lcp[rsa[i]-1]=h;
		}
	}
};
SuffixArray sa;


void solve() {
	int i,j,k,l,r,x,y; string s;
	
	cin>>T;
	while(T--) {
		cin>>N>>K;
		ll sum=0;
		G.clear();
		FOR(i,1) {
			G.resize(N);
			
			FOR(j,N) {
				cin>>G[j];
				TL[j]=N-j;
				re[j]={i+1,j+1};
				sum+=j+1;
			}
		}
	
		sa.build(G);
		int cur=0;
		vector<pair<int,int>> Rs={{0,0}};
		for(i=G.size();i>=0;i--) {
			while(Rs.back().first>sa.lcp[i]) {
				x=Rs.back().first;
				y=Rs.back().second;
				Rs.pop_back();
				k=max(Rs.back().first,sa.lcp[i]);
				sum-=1LL*(x-k)*(y-i);
				while(cur>=0&&K>sum) {
					ret[cur][0]=re[sa.sa[i+1]].first;
					ret[cur][1]=re[sa.sa[i+1]].second;
					ret[cur][2]=re[sa.sa[i+1]].second+k+(K-sum-1)/(y-i);
					cur--;
				}
				if(Rs.back().first<sa.lcp[i]) {
					Rs.push_back({sa.lcp[i],y});
				}
			}
			if(Rs.back().first<TL[sa.sa[i]]) Rs.push_back({TL[sa.sa[i]],i});
		}
		for(i=ret[0][1];i<=ret[0][2];i++) cout<<G[i-1]<<" ";
		cout<<endl;
	}
}


int main(int argc,char** argv){
	string s;int i;
	if(argc==1) ios::sync_with_stdio(false), cin.tie(0);
	FOR(i,argc-1) s+=argv[i+1],s+='\n'; FOR(i,s.size()) ungetc(s[s.size()-1-i],stdin);
	cout.tie(0); solve(); return 0;
}
0