HDU3948 The Number of Palindromes

传送门

题意：给出一个字符串，问你有多少个本质不同的回文子串。

题解：
理论上有3种做法，我只写了回文自动机和manacher。
这两个东西的坑我迟早会补上的

回文自动机：建完之后根据我抄的建树方法输出节点数-1即可。

#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
#define maxn 100005
using namespace std;
int t,n,la,ncnt,cnt[maxn];
struct node { int len,fail,ch[26];} tree[maxn];
char s[maxn];
void Init()
{
	memset(tree,0,sizeof(tree));
	tree[0].len=0,tree[1].len=-1,ncnt=1,la=0,tree[0].fail=1;
}
int Get(int i,int now)
{
	while(s[i-tree[now].len-1]!=s[i]) now=tree[now].fail;
	return now;
}
void Insert(int i)
{
	int t=Get(i,la),c=s[i]-'a';
	if(!tree[t].ch[c])
	{
		tree[++ncnt].len=tree[t].len+2;
		tree[ncnt].fail=tree[Get(i,tree[t].fail)].ch[c];
		tree[t].ch[c]=ncnt;
	}
	la=tree[t].ch[c];
	cnt[la]++;
}
int main() {
	scanf("%d",&t);
	for(int cas=1;cas<=t;cas++)
	{
		scanf("%s",s+1);
		n=strlen(s+1);
		Init();
		for(int i=1;i<=n;i++) Insert(i);
		printf("Case #%d: %d\n",cas,ncnt-1);
	}
}

manacher：每次进while时，当两端是特殊字符是hash一下丢到set里面去。好久没写hash了zz如我预处理幂的时候还没有写pw1[0]-pw2[0]=1wa了几次…

#include<cstdio>
#include<cstring>
#include<set>
#include<algorithm>
#define maxn 100005
#define maxm 26
#define P 131
#define LL unsigned long long
#define mod 998244353
using namespace std;
typedef pair<LL,LL> PLL; 
int T,p[maxn*2];
char s[maxn*2],t[maxn];
LL h1[maxn*2],h2[maxn*2],pw1[maxn*2],pw2[maxn*2];
set<PLL> st;
PLL HA(int l,int r)
{
	LL ha1=h1[r]-h1[l-1]*pw1[r-l+1];
	LL ha2=(h2[r]-h2[l-1]*pw2[r-l+1]%mod+mod)%mod;
	return PLL(ha1,ha2);
}
void manacher(int n)
{
	int id=0;
	for(int i=1;i<n;i++)
	{
		if(p[id]+id>i) p[i]=min(p[2*id-i],p[id]+id-i);
		else p[i]=1;
		while(s[i-p[i]]==s[i+p[i]])
		{ 
			p[i]++; 
			if((i+p[i]-1)%2) st.insert(HA(i-p[i]+1,i+p[i]-1));
		}
		if(id+p[id]<i+p[i]) id=i;
	}
}
void init()
{
	st.clear();
	memset(s,0,sizeof(s));
	memset(p,0,sizeof(p));
	memset(h1,0,sizeof(h1));
	memset(h2,0,sizeof(h2));
}
int main()
{
	scanf("%d",&T);
	pw1[0]=pw2[0]=1;
	for(int i=1;i<maxn*2;i++) pw1[i]=pw1[i-1]*P,pw2[i]=pw2[i-1]*P%mod;
	for(int cas=1;cas<=T;cas++)
	{
		scanf("%s",t);
		init();
		int len=strlen(t);
		int n=len*2+2;
		for(int i=0;i<len;i++) s[i*2+2]=t[i];
		s[0]=1,s[n]=2;
		h1[0]=h2[0]=1;
		for(int i=1;i<=n;i++) h1[i]=h1[i-1]*P+s[i],h2[i]=(h2[i-1]*P+s[i])%mod;
		manacher(n);
		printf("Case #%d: %d\n",cas,st.size());
	}
}

说句闲话：补文化课真是刺激惨了。