bzoj 3790（manacher+树状数组）

传送门
问题：用manacher求出回文子串长度，每条回文串转化为一条线段，要用最少的线段，覆盖区间[1,n]

求最小覆盖用“dp+树状数组/线段树维护区间最小值”或者“贪心”都可以，下面的代码用的是dp+树状数组。

#include<bits/stdc++.h>
using namespace std;
const int maxn=50003,INF=0x3f3f3f3f;
char a[maxn],b[maxn<<1];
int n,m,p[maxn<<1];
int c[maxn],cnt;
struct NODE {
    int l,r;
    friend bool operator <(const NODE &a,const NODE &b) {
        return a.r<b.r;
    }
}q[maxn<<1];//由'#'扩展出来的区间也算,必须<<1 
inline void insert(int L,int R) {
    L=(L>>1)+1,R=(R>>1)-1;
    if (L>R) return ;//ESSENTIAL!!!
    q[++cnt].l=L,
    q[cnt].r=R;
}
inline void manacher() {
    m=n<<1|1;
    for (int i=1;i<=n;++i)
        b[i<<1]=a[i],b[i<<1|1]='#';
    b[0]='+',b[1]='#',b[m+1]='-';
    int mx=0,id=0;
    for (register int i=1;i<=m;++i) {
        p[i]=(mx>i)?min(mx-i,p[(id<<1)-i]):1;
        while (b[i+p[i]]==b[i-p[i]]) ++p[i];
        insert(i-p[i],i+p[i]);
        if (mx<i+p[i]) id=i,mx=i+p[i];
    }
}
inline int query(int x) {
    if (!x) return 0;//ESSENTIAL!!!
    int ans=INF;
    for (int i=x;i<=n;i+=(i&-i))
        ans=min(ans,c[i]);
    return ans;
}
inline void update(int x,int v) {
    for (int i=x;i;i-=(i&-i))
        c[i]=min(c[i],v);
}
int main() {
//  freopen("bzoj 3790.in","r",stdin);
    while (~scanf("%s",a+1)) {
        n=strlen(a+1),cnt=0;
        for (int i=1;i<=n;++i) c[i]=INF;
        manacher();
        sort(q+1,q+cnt+1);
        int ans=INF;
        for (int i=1;i<=cnt;++i) {
            int x=query(q[i].l-1)+1;
            update(q[i].r,x);
            if (q[i].r==n) ans=min(ans,x-1);//One individual string at first
        }
        printf("%d\n",ans);
    }
    return 0;

}