【CF932G】Palindrome Partition(回文树,动态规划)
【CF932G】Palindrome Partition(回文树,动态规划)
题面
CF
翻译:
给定一个串,把串分为偶数段
假设分为了\(s1,s2,s3....sk\)
求,满足\(s_1=s_k,s_2=s_{k-1}......\)的方案数
题解
反正我是不会做
基本就是照着\(laofu\)的打了一遍(laofu太强啦)
这题分成了两个步骤
如果直接分\(k\)段我们是没法直接判断的
假设两段\(s_i,s_{k-i+1}\)
因为\(s_i=s_{k-i+1}=x_1x_2.....x_j\)
假设\(s_i\)的开始位置为\(p\)
假设原串\(S\)的长度为\(n\)
\(s_i=S[p]S[p+1]....S[p+j-1]\)
\(s_{k-i+1}=S[n-j-p+1]S[n-j-p+2]...S[n-p+1]\)
对应相等的关系是
\(S[p]=S[n-p-j+1]\)
如果\(p=1\)考虑一下特殊情况
那么就是\(S[1]=S[n-j]\)
那么,如果我们有一个串\(S'\)
\(S'=S[1]S[n]S[2]S[n-2].....\)
那么,对应到上面的相等关系里面,
就是\(S'[1..j]\) 是一个回文串
其他的每一组对应关系也是如此
所以题目转换成了:
告诉你了\(S'\)回答把\(S'\)分为若干个长度为偶数的回文串的方案数
(如果没有搞清楚可以手玩一下)
这个就是一个裸的\(dp\)了
设\(f[i]\)表示\(S'[1..i]\)划分为若干长度为偶数的回文串的方案数
得到转移方程:
\]
其中,\(S'[j,i]\)是回文串
那么,每一个\(j\)对应的位置相当于是\(S'[1..i]\)的回文后缀
构建出回文树之后沿着\(last\)跳\(fail\)就可以得到所有的\(j\)的位置
然后我们就写出来了一份代码:
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define ll long long
#define RG register
#define MAX 1000100
#define MOD 1000000007
inline int read()
{
RG int x=0,t=1;RG char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=-1,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return x*t;
}
char ch[MAX],s[MAX];
int n,f[MAX];
struct Palindromic_Tree
{
struct Node
{
int son[26];
int ff,len;
}t[MAX];
int last,tot;
void init()
{
t[tot=1].len=-1;
t[0].ff=t[1].ff=1;
}
void extend(int c,int n,char *s)
{
int p=last;
while(s[n-t[p].len-1]!=s[n])p=t[p].ff;
if(!t[p].son[c])
{
int v=++tot,k=t[p].ff;
t[v].len=t[p].len+2;
while(s[n-t[k].len-1]!=s[n])k=t[k].ff;
t[v].ff=t[k].son[c];
t[p].son[c]=v;
}
last=t[p].son[c];
}
}PT;
int main()
{
PT.init();
scanf("%s",ch+1);
n=strlen(ch+1);
if(n&1){puts("0");return 0;}
for(int i=1;i<=n;i+=2)s[i]=ch[(i+1)/2];
reverse(&ch[1],&ch[n+1]);
for(int i=2;i<=n;i+=2)s[i]=ch[(i+1)/2];
f[0]=1;
for(int i=1;i<=n;++i)
{
PT.extend(s[i]-97,i,s);
int p=PT.last;
while(p!=1)
{
if(PT.t[p].len%2==0&&PT.t[p].len>0)f[i]=(f[i]+f[i-PT.t[p].len])%MOD;
p=PT.t[p].ff;
}
}
printf("%d\n",f[n]);
return 0;
}
然后在\(CF\)上交一发:
看一看数据是啥呢?
如果我们像这份代码一样,沿着\(fail\)一路上跳
因为所有字符都相同
所以要跳\(O(n)\)次
复杂度变为了\(O(n^2)\)
完美\(TLE\)
所以肯定就不能这么做。。。
然后我就不会了
假设对于某个位置,它对应的若干回文后缀
如果他们的长度\(>len/2\)
可以证明他们的长度等差
证明?根据回文串的性质,上面的那些圈圈都是相等的串
证毕
如果把这一些看成一组
每一组都至少要\(÷ 2\)
所以至多只会有\(log\)组
这样能够保证复杂度,所以我们考虑能否一组一组转移
设\(g[x]\)表示这一组东西的和
因为一组不能直接在串中表示
所以用回文树上的节点来表示
所以,\(g[x]\)表示从节点\(x\)开始,一直到缩短的值不再是当前这个等差的位置产生的贡献
比如对于这一组,它产生的贡献就是
\(g[x]=f[j1]+f[j2]+f[j3]\)
(最底下那个是原串)
好的,假设我们当前位置是\(i\)
也就是最底下的原串是\(S'[1..i]\)
最长的回文后缀,也就是当前回文树上的\(last\)
令\(p=last\)
因为是最长的一个,所以之前肯定不会计算过这个位置的值
所以\(g[p]=f[j_1]\)
然后往上面的那些回文后缀看
这些红色的代表着关于下面的父亲对称的回文前缀
我们惊奇的发现:
它们的开始位置怎么就这么巧呢?
正好就是我要算的东西
也就是说\(g[p]+=g[p.fail]\)
当然,前提是\(p.fail\)还在这一组等差的串里面
然后\(p\)就跳到第一个不是它所在的等差的回文后缀的位置
所以\(f[i]=\sum_{p}g[p]\),前提是\(i\%2=0\)
最后,总结一下算法
对于给定的串\(S\)
把它变成\(S'=S[1]S[n]S[2]S[n-1].....\)
然后依次构建回文树
每个节点要记录一下它所在的这个等差的回文后缀在哪个节点结束
然后对于每次的\(last\)节点向上跳
跳到第一个不是自己这一段回文的位置
因为每次至少要\(÷2\),所以至多跳\(log\)次
综上,时间复杂度\(O(nlog)\)
空间复杂度\(O(n\sum)\),\(\sum\)代表字符集大小
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define ll long long
#define RG register
#define MAX 1000100
#define MOD 1000000007
inline int read()
{
RG int x=0,t=1;RG char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=-1,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return x*t;
}
char ch[MAX],s[MAX];
int n,anc[MAX],diff[MAX];
int ans[MAX],f[MAX];
struct Palindromic_Tree
{
struct Node
{
int son[26];
int ff,len;
}t[MAX];
int last,tot;
void init()
{
t[tot=1].len=-1;
t[0].ff=t[1].ff=1;
anc[0]=1;
}
void extend(int c,int n,char *s)
{
int p=last;
while(s[n-t[p].len-1]!=s[n])p=t[p].ff;
if(!t[p].son[c])
{
int v=++tot,k=t[p].ff;
t[v].len=t[p].len+2;
while(s[n-t[k].len-1]!=s[n])k=t[k].ff;
t[v].ff=t[k].son[c];
t[p].son[c]=v;
diff[v]=t[v].len-t[t[v].ff].len;
anc[v]=(diff[v]==diff[t[v].ff])?anc[t[v].ff]:t[v].ff;
}
last=t[p].son[c];
}
}PT;
int main()
{
PT.init();
scanf("%s",ch+1);
n=strlen(ch+1);
if(n&1){puts("0");return 0;}
for(int i=1;i<=n;i+=2)s[i]=ch[(i+1)/2];
reverse(&ch[1],&ch[n+1]);
for(int i=2;i<=n;i+=2)s[i]=ch[(i+1)/2];
ans[0]=1;
for(int i=1;i<=n;++i)
{
PT.extend(s[i]-97,i,s);
for(int k=PT.last;k;k=anc[k])
{
f[k]=ans[i-PT.t[anc[k]].len-diff[k]];
if(anc[k]!=PT.t[k].ff)
f[k]=(f[k]+f[PT.t[k].ff])%MOD;
if(!(i&1))ans[i]=(ans[i]+f[k])%MOD;
}
}
printf("%d\n",ans[n]);
return 0;
}
【CF932G】Palindrome Partition(回文树,动态规划)的更多相关文章
- Codeforces 932G Palindrome Partition - 回文树 - 动态规划
题目传送门 通往???的传送点 通往神秘地带的传送点 通往未知地带的传送点 题目大意 给定一个串$s$,要求将$s$划分为$t_{1}t_{2}\cdots t_{k}$,其中$2\mid k$,且$ ...
- CF932G Palindrome Partition(回文自动机)
CF932G Palindrome Partition(回文自动机) Luogu 题解时间 首先将字符串 $ s[1...n] $ 变成 $ s[1]s[n]s[2]s[n-1]... $ 就变成了求 ...
- Codeforces 932G Palindrome Partition 回文树+DP
题意:给定一个串,把串分为偶数段 假设分为$s_1,s_2,s_3....s_k$ 求满足$ s_1=s_k,s_2=s_{ k-1 }... $的方案数模$10^9+7$ $|S|\leq 10^6 ...
- HDU 6599 I Love Palindrome String (回文树+hash)
题意 找如下子串的个数: (l,r)是回文串,并且(l,(l+r)/2)也是回文串 思路 本来写了个回文树+dfs+hash,由于用了map所以T了 后来发现既然该子串和该子串的前半部分都是回文串,所 ...
- 【CF932G】Palindrome Partition 回文自动机
[CF932G]Palindrome Partition 题意:给你一个字符串s,问你有多少种方式,可以将s分割成k个子串,设k个子串是$x_1x_2...x_k$,满足$x_1=x_k,x_2=x_ ...
- 2019牛客暑期多校训练营(第六场)Palindrome Mouse 回文树+dfs
题目传送门 题意:给出一个字符串,将字符串中所有的回文子串全部放入一个集合里,去重后.问这个集合里有几对<a,b>,使得a是b的子串. 思路:一开始想偏了,以为只要求每个回文串的回文后缀的 ...
- Palindrome Partition CodeForces - 932G 回文树+DP+(回文后缀的等差性质)
题意: 给出一个长度为偶数的字符串S,要求把S分成k部分,其中k为任意偶数,设为a[1..k],且满足对于任意的i,有a[i]=a[k-i+1].问划分的方案数. n<=1000000 题解: ...
- Palindrome Mouse(2019年牛客多校第六场C题+回文树+树状数组)
目录 题目链接 题意 思路 代码 题目链接 传送门 题意 问\(s\)串中所有本质不同的回文子串中有多少对回文子串满足\(a\)是\(b\)的子串. 思路 参考代码:传送门 本质不同的回文子串肯定是要 ...
- 2019牛客暑期多校训练营(第六场)C:Palindrome Mouse(回文树+树剖)
题意:给定字符串Str,求出回文串集合为S,问S中的(a,b)满足a是b的子串的对数. 思路:开始和题解的思路差不多,维护当前后缀的每个串的最后出现位置,但是不知道怎么套“最小回文分割”,所以想到了树 ...
随机推荐
- 在Arrays.asList()引发的问题中进一步学习集合与泛型等内容
前言 最近在网上看到一个问题,情况类似如下(记为问题1): public class Demo { public static void main(String[] args) { System.ou ...
- Godep的基本使用
[http://studygolang.com/articles/4385] 关于Godep 发现好多golang项目都使用到godep作为包管理的工具,像比较大型的项目,比如kubernetes这种 ...
- html拨打电话、发送短信、发送邮件的链接写法
拨打电话 <a href="tel:88888888">呼叫</a> 发送短信 <a href="sms:88888888"> ...
- C# 使用AngleSharp 爬虫图片
AngleSharp 简介 AngleSharp是基于.NET(C#)开发的专门解析HTML源码的DLL组件.根据HTML的DOM结构操作HTML,整个DOM已传输到逻辑类结构中.这种结构可以更好的操 ...
- Servlet3.0+springmvc5+log4j2正确的开启姿势(WebLookUp)
前言 java社区占据市场份额比较大的日志组件由log4j 1.×,到logback,再到整合后的升级版 log4j 2.×,有网友测试后据说log4j2的性能最NB.于是开始往自己的springmv ...
- Vscode 插件
HTML Snippets Markdown All in One Markdown PDF Markdown Priview Enhanced Markdown TOC Open HTML in D ...
- 如何使用 VS生成动态库?
如何使用 VS生成动态库? //.cpp 文件默认定义了 __cplusplus 宏 #ifdef __cplusplus extern "C"{ #endif xxx #ifde ...
- 搭建多系统yum服务器
一.多系统服务器搭建 1.首先挂载光盘 2.安装vsftp 3.使用rpm -ql vsftpd查看vsftpd安装时都产生了哪些文件,找到以.server结尾的文件路径.此文件的文件名就是vsftp ...
- tomcat证书配置
第一步:为服务器生成证书 1.进入%JAVA_HOME%/bin目录 2.使用keytool为Tomcat生成证书,假定目标机器的域名是"localhost",keystore文件 ...
- zTree实现地市县三级级联DAO接口
zTree实现地市县三级级联DAO接口 ProvinceDao.java: /** * @Title:ProvinceDao.java * @Package:com.gwtjs.dao * @Desc ...