manacher(无讲解)

BZOJ3325: [Scoi2013]密码

https://lydsy.com/JudgeOnline/problem.php?id=3325

分析:

根据前i个字符和一些不等和相等条件就可以确定每一位。
用manacher优化暴力的过程，发现就是manacher逆过来做。
相等的赋值，不等的打标记。

代码:

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
using namespace std;
#define N 200050
#define db(x) cerr<<#x<<" = "<<x<<endl
char ans[N];
int n,a[N],b[N],vis[30][N];
int rp[N];
int main() {
    scanf("%d",&n);
    int i,j;
    for(i=1;i<=n;i++) scanf("%d",&a[i]);
    for(i=1;i<n;i++) scanf("%d",&b[i]);
    int ln=2*n+1;
    for(i=1;i<=n;i++) {
        rp[2*i-1]=b[i-1]+1;
        rp[2*i]=a[i]+1;
    }
    rp[2*n+1]=1;
    
    int lst=1,mx=1;
    for(i=1;i<=ln;i++) {
        if(i&1) ans[i]='z'+1;
        else if(!ans[i]) {
            for(j=0;j<26;j++) if(!vis[j][i]) {
                ans[i]=j+'a'; break;
            }
        }
        int t=max(1,min(mx-i,rp[lst*2-i]));
        for(;t<=rp[i];t++) {
            if(i+t-1<=ln&&i-t+1>=1) ans[i+t-1]=ans[i-t+1];
        }
        if(i-t+1>=1&&i+t-1<=ln) vis[ans[i-t+1]-'a'][i+t-1]=1;
        if(i+rp[i]-1>mx) {
            mx=i+rp[i]-1;
            lst=i;
        }
    }
    for(i=2;i<=ln;i+=2) printf("%c",ans[i]);
}

BZOJ4166: 月宫的符卡序列

https://lydsy.com/JudgeOnline/problem.php?id=4166

分析:

用manacher+hash暴力建立回文树并在上面DP即可。

代码:

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <tr1/unordered_map>
#include <map>
using namespace std;
using namespace std::tr1;
#define db(x) cerr<<#x<<" = "<<x<<endl
#define N 1000050
#define M 2000050
#define base 131
typedef unsigned long long ull;
char w[N];
int n,ln,head[N],to[N],nxt[N],cnt,fa[N],tot,s[M],rp[M];
int val[N],ans;
ull mi[N],h[N];
unordered_map<ull,int>mp;
inline void add(int u,int v) {
    to[++cnt]=v; nxt[cnt]=head[u]; head[u]=cnt;
}
ull get_h(int l,int r) {
    return h[r]-h[l-1]*mi[r-l+1];
}
int work(int l,int r) {
    ull now=get_h(l,r);
    int p;
    if(mp.count(now)) p=mp[now];
    else p=mp[now]=++tot;


    if(r-l+1<=1||fa[p]) return p;
    
    fa[p]=work(l+1,r-1);
    add(fa[p],p);
    return p;
}
void dfs(int x) {
    int i;
    for(i=head[x];i;i=nxt[i]) {
        dfs(to[i]);
        val[x]^=val[to[i]];
    }
    if(x!=1&&val[x]>ans) ans=val[x];
}
bool check(int l,int r) {
    int i,j;
    for(i=l,j=r;i<=j;i++,j--) if(w[i]!=w[j]) return 0;
    return 1;
}
void solve() {
    ans=0;
    mp.clear();
    memset(rp,0,sizeof(rp));
    memset(val,0,sizeof(val));
    memset(head,0,sizeof(head)); cnt=0;
    memset(fa,0,sizeof(fa));

    scanf("%s",w+1);
    n=strlen(w+1);
    int i;
    for(mi[0]=i=1;i<=n;i++) {
        mi[i]=mi[i-1]*base;
        h[i]=h[i-1]*base+w[i];
    }
    mp[0]=1;
    for(i=2;i<=27;i++) {
        fa[i]=1; add(1,i); mp[i-2+'a']=i;
    }
    tot=27;

    int mx=1,lst=1;
    for(i=1;i<=n;i++) {
        s[(i<<1)-1]='{';
        s[(i<<1)]=w[i];
    }
    s[n*2+1]='{';
    ln=n*2+1;

    lst=mx=1;
    for(i=1;i<=ln;i++) {
        if(i<mx) rp[i]=min(mx-i+1,rp[lst*2-i]);
        else rp[i]=1;
        for(;i-rp[i]>=1&&i+rp[i]<=ln&&s[i-rp[i]]==s[i+rp[i]];rp[i]++);
        int t=(rp[i])>>1;
        if(i&1) {
            if(i!=1&&t) {
                int mid=i>>1;
                val[work(mid-t+1,mid+t)]^=(mid-1);
            }
        }else {
            if(t) {
                int mid=i>>1;
                val[work(mid-t+1,mid+t-1)]^=(mid-1);
            }
        }
        if(i+rp[i]-1>mx) {
            mx=i+rp[i]-1; lst=i;
        }
    }
    dfs(1);
    printf("%d\n",ans);
}
int main() {
    int T;
    scanf("%d",&T);
    while(T--) solve();
}

2342: [Shoi2011]双倍回文

https://lydsy.com/JudgeOnline/problem.php?id=2342

分析:

首先我们选择的一定是一个回文串。
这个限制去掉之后，只需令它的一半也是回文串。
假设对称轴左边位置为\(x\)，且\(y\)能成为其一半子串的右端点的条件是\(y-rp_y\le x\)且\(y\le x+rp_x/2\) 。
按\(y-rp_y\)排序然后用\(set\)维护前驱后继即可。

代码:

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <set>
using namespace std;
#define N 1000050
char w[N],s[N];
int n,rp[N],f[N],ans,t[N];
set<int>S;
inline bool cmp(const int &x,const int &y) {return x-f[x] < y-f[y];}
int main() {
    scanf("%d%s",&n,w+1);
    int i,ln=n*2+1;
    for(i=1;i<=n;i++) s[2*i-1]='{', s[2*i]=w[i];
    s[2*n+1]='{';
    int mx=1,lst=1;
    for(i=1;i<=ln;i++) {
        rp[i]=min(mx-i,rp[lst*2-i]);
        if(!rp[i]) rp[i]=1;
        for(;i-rp[i]>=1&&i+rp[i]<=ln&&s[i-rp[i]]==s[i+rp[i]];rp[i]++) ;
        if(i+rp[i]-1>mx) {
            mx=i+rp[i]-1, lst=i;
        }
    }
    for(i=1;i<n;i++) {
        f[i]=rp[i*2+1]>>1;
    }
    int j=1;
    for(i=1;i<n;i++) t[i]=i;
    sort(t+1,t+n,cmp);
    for(i=1;i<n;i++) {
        int x=i;
        for(;j<n&&t[j]-f[t[j]]<=x;j++) S.insert(t[j]);
        set<int>::iterator it=S.upper_bound((f[x]>>1)+x);
        if(it!=S.begin()) {
            --it;
            ans=max(ans,((*it)-x)*4);
        }
    }
    printf("%d\n",ans);
}

BZOJ5036: [Jsoi2014]回文串

https://lydsy.com/JudgeOnline/problem.php?id=5036

分析:

先用其他字符将每个字符隔开方便统计。
考虑\([l,r]\)这段区间以\(x\)为中心的回文串个数。
显然这玩意等于\(min(x-l+1,rp[x],r-x+1)\)
拆成两部分，每一部分用树状数组再拆开统计。

代码:

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cstdlib>
using namespace std;
typedef long long ll;
#define N 200050
#define M 300050
char w[N],s[N];
int n,rp[N],pi[N],Q,ln;
ll ans[M];
struct A {
    int id,l,r,mid;
}q[M];
struct B {
    ll c[N];
    void fix(int x,int v) {for(;x<=ln;x+=x&(-x)) c[x]+=v;}
    ll inq(int x) {ll re=0;for(;x;x-=x&(-x)) re+=c[x]; return re;}
}t1,t2,t3;
bool cmp1(const A &x,const A &y) {return x.l<y.l;}
bool cmp2(const A &x,const A &y) {return x.r>y.r;}
bool cmp3(const int &x,const int &y) {return x-rp[x]<y-rp[y];}
bool cmp4(const int &x,const int &y) {return x+rp[x]>y+rp[y];}
int main() {
    scanf("%s%d",w+1,&Q);
    n=strlen(w+1);
    int i; ln=2*n+1;
    for(i=1;i<=n;i++) s[2*i-1]='{',s[2*i]=w[i];
    s[2*n+1]='{';
    int mx=1,lst=1;
    for(i=1;i<=ln;i++) {
        rp[i]=min(mx-i,rp[2*lst-i]);
        if(!rp[i]) rp[i]=1;
        for(;i-rp[i]>=1&&i+rp[i]<=ln&&s[i-rp[i]]==s[i+rp[i]];rp[i]++);
        if(i+rp[i]-1>mx) {
            mx=i+rp[i]-1; lst=i;
        }
    }
    // for(i=1;i<=ln;i++) printf("%d ",rp[i]); puts("");
    for(i=1;i<=Q;i++) {
        scanf("%d%d",&q[i].l,&q[i].r);
        q[i].l=2*q[i].l-1, q[i].r=2*q[i].r+1;
        q[i].id=i; q[i].mid=(q[i].l+q[i].r)>>1;
    }
    sort(q+1,q+Q+1,cmp1);
    for(i=1;i<=ln;i++) {
        t1.fix(i,rp[i]);
        pi[i]=i;
    }
    sort(pi+1,pi+ln+1,cmp3);
    int j=1;
    for(i=1;i<=Q;i++) {
        for(;j<=ln&&pi[j]-rp[pi[j]]+1<q[i].l;j++) {
            t1.fix(pi[j],-rp[pi[j]]);
            t2.fix(pi[j],1);
            t3.fix(pi[j],pi[j]);
        }
        ans[q[i].id]+=t1.inq(q[i].mid)-t1.inq(q[i].l-1)
                    - (t2.inq(q[i].mid)-t2.inq(q[i].l-1))*(q[i].l-1)
                    + (t3.inq(q[i].mid)-t3.inq(q[i].l-1));
    }
    memset(t1.c,0,sizeof(t1.c));
    memset(t2.c,0,sizeof(t2.c));
    memset(t3.c,0,sizeof(t3.c));
    sort(q+1,q+Q+1,cmp2);
    for(i=1;i<=ln;i++) t1.fix(i,rp[i]);
    sort(pi+1,pi+ln+1,cmp4);
    j=1;
    for(i=1;i<=Q;i++) {
        for(;j<=ln&&pi[j]+rp[pi[j]]-1>q[i].r;j++) {
            t1.fix(pi[j],-rp[pi[j]]);
            t2.fix(pi[j],1);
            t3.fix(pi[j],pi[j]);
        }
        ans[q[i].id]+=t1.inq(q[i].r)-t1.inq(q[i].mid)
                    + (t2.inq(q[i].r)-t2.inq(q[i].mid))*(q[i].r+1)
                    - (t3.inq(q[i].r)-t3.inq(q[i].mid));
    }

    for(i=1;i<=Q;i++) {
        ans[q[i].id]=(ans[q[i].id]-((q[i].r-q[i].l)>>1)-1)>>1;
    }
    for(i=1;i<=Q;i++) printf("%lld\n",ans[i]);
}

BZOJ3325: [Scoi2013]密码

BZOJ4166: 月宫的符卡序列

2342: [Shoi2011]双倍回文

BZOJ5036: [Jsoi2014]回文串

猜你喜欢