HDU5442

HDU5442

---

做法：把原串复制一份加在后边，中间插特殊入个特殊字符，再把翻转后的串加在后边，同样复制一份。然后做后缀数组，按题意处理细节即可。

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <map>
#include <cstring>
#include <cmath>
#include <queue>
#include <set>
#include <vector>
#include <iterator>
#include <string>
#include <deque>
#define rep(i,a,b) for(int i=a;i<=b;++i)
#define per(i,a,b) for(int i=a;i>=b;--i)
#define pb push_back
#define MP make_pair
#define fr first
#define sc second
#define PII pair<int,int>
#define VI vector<int>
typedef long long ll;
typedef unsigned long long ull;
const int N = 100005;
inline int readint() {
    char c=getchar();int x=0,f=1;
    while(!isdigit(c)){if(f=='-')f=-1;c=getchar();}
    while(isdigit(c)){x=x*10+c-'0';c=getchar();}
    return x*f;
}
using namespace std;
int n;
char c[N],str[N];
int rnk[N] , SA[N] , Height[N];
int X[N] , Y[N] , sum[N];
int f[101000][20] , fm[101000][20];
bool cmp(int *r,int a,int b,int l) {
    return ( r[a] == r[b] && r[a+l] == r[b+l] );
}
void calc() {
    int l , p , *x = X , *y = Y , m = 128;
    rep(i,0,m) sum[i] = 0;
    rep(i,1,n) sum[ x[i] = c[i] ] ++;
    rep(i,1,m) sum[i] += sum[i-1];
    per(i,n,1) SA[ sum[ x[i] ]-- ] = i;
    for ( l = 1 , p = 1 ; l <= n ; m = p , l *= 2 ) {
        p = 0;
        rep(i,n-l+1,n) y[++p] = i;
        rep(i,1,n) if ( SA[i] > l ) y[++p] = SA[i] - l;
        rep(i,0,m) sum[i] = 0;
        rep(i,1,n) sum[ x[y[i]] ] ++;
        rep(i,1,m) sum[i] += sum[i-1];
        per(i,n,1) SA[ sum[ x[y[i]] ]-- ] = y[i];
        swap( x , y );
        x[SA[1]] = 1; p = 2;
        rep(i,2,n)
            x[ SA[i] ] = cmp(y,SA[i-1],SA[i],l) ? p - 1 : p++;
    }
    rep(i,1,n) rnk[SA[i]] = i;
    p = 0;
    rep(i,1,n) {
        if ( rnk[i] == 1 ) continue;
        while ( c[i+p] == c[SA[rnk[i]-1]+p] ) p ++;
        Height[rnk[i]] = p;
        if ( p ) p --;
    }
}
void init() {
    n = strlen(str);
    int cc = 0;
    for(int i=0;i<n;++i) c[++cc] = str[i];
    for(int i=0;i<n;++i) c[++cc] = str[i];
    c[++cc] = '$';
    for(int i=n-1;i>=0;--i) c[++cc] = str[i];
    for(int i=n-1;i>=0;--i) c[++cc] = str[i];
    c[cc+1]=0;
    n = cc;
}
int Log[N],rmq[N][30];
void init_rmq() {
    Log[1] = 0;
    for(int i=2;i<=n;++i) Log[i] = Log[i>>1] + 1;
    for(int i=1;i<=n;++i) rmq[i][0] = Height[i];
    for(int j=1;j<=20;++j)
        for(int i=1;i+(1<<(j-1))<=n;++i)
            rmq[i][j] = min(rmq[i][j-1],rmq[i+(1<<j-1)][j-1]);
}
int RMQ_mn(int l,int r){
    int L=Log[r-l+1];
    return min(rmq[l][L],rmq[r-(1<<L)+1][L]);
}
int ask(int x,int y) {
    x=rnk[x],y=rnk[y];
    if(x>y)swap(x,y);
    return RMQ_mn(x+1,y);
}
int T, num;
int biao(int x) {
    if(x <= num) return x;
    return num - (x-num*2-1) + 1;
}
int cal_s(int x) {
    if(x<=num) return 0;
    return 1;
}
int main() {
    scanf("%d",&T);
    while(T--) {
        scanf("%d",&num);
        scanf(" %s",str);
        init();
        calc();
        init_rmq();
        int idx = n;
        int ansb = biao(SA[n]), anss = cal_s(SA[n]);
        for(int i=n;i>=1;--i) if(SA[i]!=num*2+1){
            if(SA[i]<=num) {
                idx = i;
                ansb = biao(SA[i]);
                anss = cal_s(SA[i]);
                break;
            }
            if(SA[i]>=2*num+1&&SA[i]<=3*num+1) {
                idx = i;
                ansb = biao(SA[i]);
                anss = cal_s(SA[i]);
                break;
            }
        }
        for(int i=idx-1;i>=1;--i) if(SA[i] != num*2+1) {
            if(SA[i]>num&&SA[i]<=num*2+1) continue;
            if(SA[i]>3*num+1) continue;
            int lcp = ask(SA[idx],SA[i]);
            if(lcp >= num) {
                int b = biao(SA[i]), s = cal_s(SA[i]);
                if(b<ansb) ansb = b,anss=s;
                else if(b==ansb&&s < anss) ansb=b,anss=s;
            }
        }
        printf("%d %d
",ansb,anss);
    }
    return 0;
}