bzoj 3790 神奇项链

最新推荐文章于 2018-02-01 16:17:00 发布

heheda_is_an_OIer

最新推荐文章于 2018-02-01 16:17:00 发布

阅读量390

点赞数

分类专栏：字符串

本文链接：https://blog.csdn.net/heheda_is_an_OIer/article/details/49862599

版权

字符串专栏收录该内容

7 篇文章 0 订阅

订阅专栏

这道题可以转化为用最少的回文子串覆盖这个字符串。比较神奇的是树状数组是可以维护后缀的信息的。具体方法是

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;
}