bzoj 3790 神奇项链

这道题可以转化为用最少的回文子串覆盖这个字符串。比较神奇的是树状数组是可以维护后缀的信息的。具体方法是
void modify(int x,int d)
{
for (;x;x-=x&(-x)) c[x]=min(c[x],d);
}
int query(int x)
{
if (x==0) return 0;
int ans=inf;
for (;x<=n;x+=x&(-x)) ans=min(ans,c[x]);
return ans;
}
简单的说就是把modify和query的代码一交换

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>

#define md
#define ll long long
#define inf (int) 1e9
#define eps 1e-8
#define N 200010
using namespace std;
struct QQ {int l,r;} q[N];
int c[N],p[N];
int n,m;
char s[N],st[N];
bool cmp(QQ a,QQ b) { return a.r<b.r;}
void modify(int x,int d)
{
for (;x;x-=x&(-x)) c[x]=min(c[x],d);
}
int query(int x)
{
if (x==0) return 0;
int ans=inf;
for (;x<=n;x+=x&(-x)) ans=min(ans,c[x]);
return ans;
}

int manacher(int l)
{
m=2*l+1;
for (int i=1;i<=l;i++)
{
st[i<<1]=s[i]; st[(i<<1)|1]='#';
}
st[0]='+'; st[m+1]='-'; st[1]='#';
int mx=0,id=0;
memset(p,0,sizeof(p));
for (int i=1;i<=m;i++)
{
if (i<mx) p[i]=min(p[2*id-i],mx-i);
else p[i]=1;
while (st[i+p[i]]==st[i-p[i]]) p[i]++;
if (i+p[i]>mx) { mx=i+p[i]; id=i;}
}
int w=0;
for (int i=1;i<=m;i++)
{
w++;
q[w].l=(i-p[i])/2+1; q[w].r=(i+p[i])/2-1;
if (q[w].l>q[w].r) w--;
}
return w;
//for (int i=1;i<=m;i++) printf()
}



int main()
{
while (scanf("%s",s+1)!=EOF)
{
n=strlen(s+1);
m=manacher(n);
sort(q+1,q+m+1,cmp);
//for (int i=1;i<=m;i++) printf("%d %d\n",q[i].l,q[i].r);
for (int i=1;i<=n;i++) c[i]=inf;
for (int i=1;i<=m;i++)
{
int f=query(q[i].l-1)+1;
modify(q[i].r,f);
}
int ans=query(n);
printf("%d\n",ans-1);
}
return 0;
}



你可能感兴趣的:(bzoj 3790 神奇项链)