【HDU5442】 Favorite Donut(后缀数组)

题意,给你一个长度为n的字符串,让你选择从某一个地方断开(可以正着取也可以反着取)。问你最大的字典序是从哪里断开。

首先先将原串复制一遍,然后用后缀数组求出字典序最大的位置,如果字典序相同则求出最靠前的位置。

然后再将该串翻转,得到反向取字典序最大的串,再直接比较就好了。

code:

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N = 40005;
char s[N],ss[N];
int sa[N],t[N],t2[N];
int Rank[N],height[N],c[300];
void build_sa(int n,int m){
    int i,*x = t,*y = t2;
    for(i = 0 ; i < m ; i ++) c[i] = 0;
    for(i = 0 ; i < n ; i ++) c[ x[i] = s[i] ] ++;
    for(i = 1 ; i < m ; i ++) c[i] += c[i - 1];
    for(i = n - 1 ; i >= 0 ; i --) sa[--c[x[i]]] = i;
    for(int k = 1 ; k <= n ; k <<= 1){
        int p = 0;
        for(int i = n - k ; i < n ; i ++) y[p ++] = i;
        for(int i = 0 ; 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 = 0 ; i < n ; i ++) c[ x[i] ] ++;
        for(int i = 1 ; i < m ; i ++) c[i] += c[i - 1];
        for(int i = n - 1 ; i >= 0 ; i --) sa[--c[x[y[i]]]] = y[i];
        swap(x,y);
        p = 1;x[sa[0]] = 0;
        for(int i = 1 ; i < n ; i ++)
            x[sa[i]] = y[sa[i - 1]] == y[sa[i]] && y[sa[i - 1] + k] == y[sa[i] + k] ? p - 1 : p ++; 
        if(p >= n) break;
        m = p;
    }
}
void getheight(int n){
    int i,j,k = 0;
    for(int i = 1 ; i <= n ; i ++) Rank[sa[i]] = i;
    for(int i = 0 ; i < n ; i ++){
        if(k) k--;
        int j = sa[ Rank[i] - 1 ];
        while(s[i + k] == s[j + k]) k ++;
        height[Rank[i]] = k;
    }
}
int n;
void debug(){
    for(int i = 0 ; i <= n ; i ++) printf("%d ",sa[i]);
    printf("\n");
    for(int i = 0 ; i <= n ; i ++) printf("%d ",Rank[i]);
    printf("\n");
    for(int i = 0 ; i <= n ; i ++) printf("%d ",height[i]);
    printf("\n");
}
void cmp(int s1,int s2){
    int flag = 0;
    for(int i = 0 ; i < n ; i ++){
        if(ss[(s1 + i) % n] == ss[(s2 - i + n) % n]) continue;
        if(ss[(s1 + i) % n] > ss[(s2 - i + n) % n]) flag = 1;
            else flag = -1;
        if(flag != 0) break;
    }
    s1 ++,s2 ++;
    if(flag == 0){
        if(s1 <= s2) printf("%d 0\n",s1);
            else printf("%d 1\n",s2);
    }
    else if(flag > 0) printf("%d 0\n",s1);
    else printf("%d 1\n",s2);
}
void solve(){
    scanf("%d%s",&n,s);
    for(int i = 0 ; i < n ; i ++) ss[i] = s[n + i] = s[i];
    s[n * 2] = 0;
    build_sa(n * 2 + 1,127);
    getheight(n * 2);
    int flag = 1,st1,st2;
    st1 = st2 = sa[0];
    for(int i = 1 ; i <= n * 2 ; i ++){
        if(height[i] < n) flag = 0;
        if(sa[i] >= n) continue;
        if(flag) st1 = min(st1,sa[i]);
        else{
            flag = 1;
            st1 = sa[i];
        }
    }
    for(int i = 0 ; i < n ; i ++) swap(s[2 * n - 1 - i],s[i]);
    s[n * 2] = 0;
    build_sa(n * 2 + 1,127);
    getheight(n * 2);
    flag = 1;
    for(int i = 1 ; i <= n * 2 ; i ++){
        if(height[i] < n) flag = 0;
        if(sa[i] >= n) continue;
        if(flag) st2 = min(st2,n - 1 - sa[i]);
        else{
            flag = 1;
            st2 = n - 1 - sa[i];
        }
    }
    cmp(st1,st2);
}
int main(){
    int _;
    scanf("%d",&_);
    while(_--) solve();
    return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值