关闭

[后缀数组 后缀树] Codechef January Challenge 2018 #KILLKTH Killjee and k-th letter

86人阅读 评论(0) 收藏 举报
分类:

建出后缀树,记录每个子串的出现次数,然后二分下答案在哪个子串中就好了

退役选手不会写后缀自动机

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<vector>
#define pb push_back
using namespace std;
typedef long long ll;

inline char nc(){
  static char buf[100000],*p1=buf,*p2=buf;
  return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline void read(int &x){
  char c=nc(),b=1;
  for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
  for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc());
}
inline void read(ll &x){
  char c=nc(),b=1;
  for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
  for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc());
}
inline int read(char *s){
  char c=nc(); int len=0;
  for (;!(c>='a' && c<='z');c=nc());
  for (;c>='a' && c<='z';s[++len]=c,c=nc()); s[++len]=0; return len-1;
}

const int N=400005;

int n; char S[N];

int sa[N],t1[N],t2[N],c[N];  
int height[N],rnk[N];  


inline void SA(char *r,int m){  
  int *x=t1,*y=t2;  
  for (int i=0;i<=m;i++) c[i]=0;  
  for (int i=1;i<=n;i++) c[x[i]=r[i]-'a'+1]++;  
  for (int i=1;i<=m;i++) c[i]+=c[i-1];  
  for (int i=n;i;i--) sa[c[x[i]]--]=i;  
  for (int k=1;k<=n;k<<=1){  
    int p=0;  
    for (int i=n-k+1;i<=n;i++) y[++p]=i;  
    for (int i=1;i<=n;i++) if (sa[i]>k) y[++p]=sa[i]-k;  
    for (int i=0;i<=m;i++) c[i]=0;  
    for (int i=1;i<=n;i++) c[x[y[i]]]++;  
    for (int i=1;i<=m;i++) c[i]+=c[i-1];  
    for (int i=n;i;i--) sa[c[x[y[i]]]--]=y[i];  
    swap(x,y);  
    x[sa[1]]=1; p=1;  
    for (int i=2;i<=n;i++)  
      x[sa[i]]=(y[sa[i]]==y[sa[i-1]] && y[sa[i]+k]==y[sa[i-1]+k])?p:++p;  
    if (p>=n) break;  
    m=p;  
  }  
}   
inline void GetHeight(){  
  int j,k=0;  
  for (int i=1;i<=n;i++) rnk[sa[i]]=i;  
  for (int i=1;i<=n;height[rnk[i++]]=k)  
    for (k?k--:0,j=sa[rnk[i]-1];S[i+k]==S[j+k];)  
      k++;  
}

#define V G[p].v

struct edge{
  int u,v,next;
}G[N];
int head[N],inum;
inline void add(int u,int v,int p){
  G[p].u=u; G[p].v=v; G[p].next=head[u]; head[u]=p;
}

int depth[N],ed[N],cnt[N];
int fat[N];

int tot;
int sta[N],pnt;
inline void Build(){
  sta[++pnt]=++tot; sta[++pnt]=++tot;
  depth[tot]=n-sa[1]+1; ed[tot]=n; cnt[tot]++;
  for (int i=2;i<=n;i++){
    int h=height[i],last=0;
    while (depth[sta[pnt]]>h){
      if (last) add(sta[pnt],last,++inum);
      last=sta[pnt--];
    }
    if (depth[sta[pnt]]<h){
      sta[++pnt]=++tot;
      depth[tot]=h;
      ed[tot]=ed[last]-(depth[last]-h);
    }
    if (last) add(sta[pnt],last,++inum);
    if (n-sa[i]+1==h)
      cnt[tot]++;
    else{
      sta[++pnt]=++tot;
      depth[tot]=n-sa[i]+1;
      ed[tot]=n;
      cnt[tot]=1;
    }
  }
  while (pnt>1) add(sta[pnt-1],sta[pnt],++inum),pnt--;
}
int lst[N];
ll Cnt[N];
inline void dfs(int u,int fa){
  lst[++*lst]=u; fat[u]=fa;
  vector<int> v;
  for (int p=head[u];p;p=G[p].next)
    v.pb(V);
  for (int i=v.size()-1;~i;i--)
    if (v[i]!=fa)
      dfs(v[i],u),cnt[u]+=cnt[v[i]];
}

inline ll SS(int l,int r){
  return (ll)(l+r)*(r-l+1)/2;
}

int main(){
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  n=read(S); SA(S,26); GetHeight();
  Build();
  dfs(1,0);
  for (int i=2;i<=tot;i++)
    Cnt[i]=cnt[lst[i]]*(SS(depth[fat[lst[i]]]+1,depth[lst[i]]))+Cnt[i-1];
  int Q; read(Q); ll G=0,P,M,K;
  while (Q--){
    read(P); read(M); K=P*G%M+1;
    int it=lower_bound(Cnt+2,Cnt+tot+1,K)-Cnt,p=lst[it];
    K-=Cnt[it-1];
    int L=depth[fat[p]],R=depth[p],C=cnt[p],MID;
    while (L+1<R)
      if (K>SS(depth[fat[p]]+1,MID=(L+R)>>1)*C)
    L=MID;
      else
    R=MID;
    K-=SS(depth[fat[p]]+1,R-1)*C;
    K=(K-1)%R+1;
    char ans=S[ed[p]-depth[p]+K];
    G+=(int)ans;
    putchar(ans); putchar('\n');
  }
  return 0;
}
0
0
查看评论

[Contest] Codechef January Challenge 2018

发烧在家躺了一个星期呢… Rectangle 直接判一下是不是有两对相等的 Maximum Score 倒过来贪心,每次选最大的能选的 K-Concatenation 分三种情况讨论就完了,在一个串呢,在两个串呢,跨越n-2个串 Partition the number...
  • Coldef
  • Coldef
  • 2018-01-16 09:14
  • 96

[树形DP][后缀自动机][后缀树] CodeChef January Challenge 2018 KILLKTH

SolutionSolution 反串的后缀自动机的parent树就是原串的后缀树。 按dfs序记录子树大小前缀和,二分位置。 定位后就可以在节点内二分(可以直接算的啦) #include #define show(x) cerr using namespace std; typed...
  • Vectorxj
  • Vectorxj
  • 2018-01-23 21:04
  • 44

[Contest]CodeChef January Challenge 2018

RECTANGL:学习了ifif语句的用法。 MAXSC:贪心。 KCON:贪心。 PRTITION:乱搞题。 STRMRG:O(n2)O(n^2)DP。 MONSTER:整体二分+暴力。。 XYHUMOQ:奇妙的搜索。。设f0,f1,g0,g1f_0,f_1,g_0,g_1分别表示前(...
  • Vectorxj
  • Vectorxj
  • 2018-01-15 20:09
  • 71

[杜教筛] Codechef January Challenge 2018 #SQRGOOD Simplify the Square Root

二分转化为μ2\mu^2的前缀和。 然后转化为O(n13)O(n^{1\over 3})的运算,但是需要预处理μ\mu的前缀和,大力杜教筛求和。 然后感谢阿爷教我把二分改成了迭代,小范围内一个一个挪,用rho求μ(n)\mu(n),然后就能卡进去了。 复杂度似乎是萎的吧。 #include ...
  • u014609452
  • u014609452
  • 2018-01-19 14:16
  • 75

[整体二分] Codechef January Challenge 2018 #MONSTER Killing Monsters

整体二分,然后问题变成,子集加,单点查询,然后像CTSC吉夫特 可以用经典的二进制分高位低位的搞搞。调个参,大概是高5位低12位。 不知道在线怎么做。 #include #include #include #include #define pb push_back using namespa...
  • u014609452
  • u014609452
  • 2018-01-19 14:11
  • 104

USACO 2018 January Contest Platinum A: Lifeguards 题解

将所有的区间按左端点从小到大排序 我们处理那些被完全包含的区间,这些区间即使删除也不会使答案变坏 这样先删一波,如果发现这种小区间的个数多于k,就可以直接算答案了 否则我们要dp 设dp[i][j]为考虑到第i个区间,已经删除了j个区间,且第i个区间没有被删除的情况下最大的覆盖长度 显然有...
  • IcePrincess_1968
  • IcePrincess_1968
  • 2018-01-23 23:58
  • 82

实用算法实现-第8篇 后缀树和后缀数组 [1简介]

8.1    后缀树 一棵后缀树包含一个指定文本的所有后缀,对于在一个长度为N的文本中查找一个长度为M的子字符串,一个后缀树仅仅需要M次比较,而这个比较次数是查找该字符串所需要的最小比较次数。 后缀树有以下特征:一条边可以表示文本的任何非空子串,每个非终端节点,除了...
  • fsdev
  • fsdev
  • 2011-10-15 11:50
  • 3798

Codechef 2017 March Challenge 简要题解

在比赛进行到一半的时候才参加,,导致challenge没法搞得很好       写一个简要题解吧,,challenge太差了就不放了 Xenny and Alternating Tasks:直接按照题意模拟即可。 Bear and Extra Number:首先处理...
  • lych_cys
  • lych_cys
  • 2017-03-13 21:01
  • 752

利用后缀数组构造后缀树

由于蒟蒻azui前段时间忙着准备省选,并在省选中成功闷声滚大粗,博客停更了好久。。
  • u011542204
  • u011542204
  • 2016-04-24 19:06
  • 1676

CodeChef March Challenge 2017 题解

CC XENTASK由于要交替选,所以要不第一个人选位置为奇数的数,第二个人选位置为偶数的数,要不第一个人选位置为偶数的数,第二个人选位置为奇数的数。取最小值即可。#include <cstring> #include <cstdio> #include <algori...
  • YxuanwKeith
  • YxuanwKeith
  • 2017-03-13 20:26
  • 514
    个人资料
    • 访问:364252次
    • 积分:13016
    • 等级:
    • 排名:第1238名
    • 原创:975篇
    • 转载:3篇
    • 译文:0篇
    • 评论:56条
    文章分类
    最新评论