题目传送门

题目大意:给你一个字符串,让你求出有多少对相交的回文子串

啊啊啊啊降智了,我怎么又忘了正难则反

求相交会很难搞。把问题转化成求互不相交的回文子串再减一下就行了

先利用$PAM$求出以每个位置为末尾的回文子串数量,这个数量就是此时构造末尾节点在$fail$树中的深度

再把串翻过来,用同样的方法求出每个位置为开头的回文子串数量

对其中一个数组求前缀和,用乘法原理算一下就行了

空间开不下怎么办?以时间换空间,用邻接表存儿子!每次跳儿子都暴力遍历一次邻接表

 #include <cstdio>
#include <cstring>
#include <algorithm>
#define ll long long
#define N1 2000010
using namespace std;
const int p=; template <typename _T> void read(_T &ret)
{
ret=; _T fh=; char c=getchar();
while(c<''&&c>''){ if(c=='-') fh=-; c=getchar(); }
while(c>=''&&c<=''){ ret=ret*+c-''; c=getchar(); }
ret=ret*fh;
}
/*void exgcd(ll a,ll b,ll &x,ll &y)
{
if(!b){ x=1; y=0; return; }
exgcd(b,a%b,x,y); ll t=x; x=y; y=t-a/b*y;
}*/ int idx(char c){ return c-'a'; } struct Edge{
int to[N1],nxt[N1],val[N1],head[N1],cte;
void ae(int u,int v,int w)
{ cte++; to[cte]=v; nxt[cte]=head[u]; val[cte]=w; head[u]=cte; }
}e; namespace PAM{
int pre[N1],dep[N1],sum[N1],sz[N1],la,tot;
void clr()
{
memset(pre,,sizeof(pre)); memset(dep,,sizeof(dep));
memset(sum,,sizeof(sum)); memset(&e,,sizeof(e));
}
void init(){ la=tot=; pre[]=pre[]=; dep[]=-; }
int same(char *str,int p,int i){ return str[i-dep[p]-]==str[i]; }
int trs(int x,int c)
{
for(int j=e.head[x];j;j=e.nxt[j])
if(e.val[j]==c) return e.to[j];
return ;
}
int insert(char *str,int i)
{
int p=la,np,fp,tp,c=idx(str[i]);
while(!same(str,p,i)) p=pre[p];
if(!(tp=trs(p,c))) //!trs[p][c]
{
np=++tot;
dep[np]=dep[p]+;
fp=pre[p];
while(!same(str,fp,i)) fp=pre[fp];
pre[np]=trs(fp,c); //trs[fp][c]
e.ae(p,np,c); //trs[p][c]=np;
p=np;
}else p=tp;
la=p;
sum[p]=sum[pre[p]]+;
return p;
}
}; int n;
char str[N1];
ll lsum[N1],rsum[N1]; int main()
{
int i,j,x; ll ans=;
scanf("%d",&n);
scanf("%s",str+);
PAM::init();
for(i=;i<=n;i++) x=PAM::insert(str,i), lsum[i]=PAM::sum[x];// ans+=lsum[i];
PAM::clr(); PAM::init();
reverse(str+,str+n+);
for(i=;i<=n;i++) x=PAM::insert(str,i), rsum[i]=PAM::sum[x], rsum[i]+=rsum[i-]; //rsum[i]+=rsum[i-1];
if(rsum[n]&) ans=1ll*((rsum[n]-)/%p)*(rsum[n]%p)%p;
else ans=1ll*(rsum[n]/%p)*((rsum[n]-)%p)%p;
for(i=;i<=n;i++) ans=(ans-1ll*lsum[i-]*rsum[n-i+]%p+p)%p;
printf("%I64d\n",(ans%p+p)%p);
return ;
}

CF17E Palisection (回文自动机+DP)的更多相关文章

  1. bzoj2084/luoguP3501 [Poi2010]Antisymmetry(回文自动机+dp)

    bzoj2084/luoguP3501 [Poi2010]Antisymmetry(回文自动机+dp) bzoj Luogu 对于一个01字符串,如果将这个字符串0和1取反后,再将整个串反过来和原串一 ...

  2. bzoj4044/luoguP4762 [Cerc2014]Virus synthesis(回文自动机+dp)

    bzoj4044/luoguP4762 [Cerc2014]Virus synthesis(回文自动机+dp) bzoj Luogu 你要用ATGC四个字母用两种操作拼出给定的串: 1.将其中一个字符 ...

  3. 洛谷P4762 [CERC2014]Virus synthesis(回文自动机+dp)

    传送门 回文自动机的好题啊 先建一个回文自动机,然后记$dp[i]$表示转移到$i$节点代表的回文串的最少的需要次数 首先肯定2操作越多越好,经过2操作之后的串必定是一个回文串,所以最后的答案肯定是由 ...

  4. bzoj 4044: [Cerc2014] Virus synthesis【回文自动机+dp】

    建回文自动机,注意到一个回文串是可以通过一个长度小于等于这个串长度的一半的回文串添上一些字符然后复制得到的,也就是在自动机上向fa走,相当于treedp 每次都走显然会T,记录一个up,指向祖先中最下 ...

  5. [CERC2014]Virus synthesis【回文自动机+DP】

    [CERC2014]Virus synthesis 初始有一个空串,利用下面的操作构造给定串 SS . 1.串开头或末尾加一个字符 2.串开头或末尾加一个该串的逆串 求最小化操作数, \(|S| \l ...

  6. BZOJ 4044 Virus synthesis (回文自动机+dp)

    题目大意: 你可以在一个串的开头或者末尾加入一个字符,或者把当前整个串$reverse$,然后接在前面或者后面,求达到目标串需要的最少操作次数 对目标串建出$PAM$ 定义$dp[x]$表示当前在回文 ...

  7. [BZOJ4044]Virus synthesis 回文自动机的DP

    4044: [Cerc2014] Virus synthesis Time Limit: 20 Sec  Memory Limit: 128 MB Description Viruses are us ...

  8. BZOJ 4044 Luogu P4762 [CERC2014]Virus Synthesis (回文自动机、DP)

    好难啊..根本不会做..基本上是抄Claris... 题目链接: (bzoj)https://www.lydsy.com/JudgeOnline/problem.php?id=4044 (luogu) ...

  9. BZOJ 3790 神奇项链(回文自动机+线段树优化DP)

    我们预处理出来以i为结尾的最长回文后缀(回文自动机的构建过程中就可以求出)然后就是一个区间覆盖,因为我懒得写贪心,就写了线段树优化的DP. #include<iostream> #incl ...

随机推荐

  1. 【Android自己定义View实战】之自己定义超简单SearchView搜索框

    [Android自己定义View实战]之自己定义超简单SearchView搜索框 这篇文章是对之前文章的翻新,至于为什么我要又一次改动这篇文章?原因例如以下 1.有人举报我抄袭,原文链接:http:/ ...

  2. Android 自己定义Activity基类

    我们在开发App的时候有时候碰到多个界面有一个共同点的时候.比方,都有同样的TitleBar.而且TitleBar能够设置显示的文字.TitleBar上的点击事件,假设给每个Activity都写一遍T ...

  3. 箭头函数的 this

    在非严格模式下,匿名函数和定时器中的 this 由于没有默认的宿主对象,因此指向 window: 而在严格模式下,匿名函数和定时器中的 this 由于没有默认的宿主对象,因此为 undefined. ...

  4. rel='canonical'

    rel='canonical' <!DOCTYPE html><html><head> <meta charset="utf-8"/> ...

  5. 20170410 --- Linux备课资料 --- 压缩与解压缩

    这节课我们来学习一下压缩与解压缩,那什么是压缩与解压缩呢? 联想一下Windows系统: 选中文件,右键选择即可 如果压缩,可以选择要压缩的格式,而解压缩直接选择就可以完成了 Linux是通过命令的方 ...

  6. ios9--UIImageView的帧动画

    // // ViewController.m // 05-UIImageView的帧动画 // #import "ViewController.h" @interface View ...

  7. tiny4412 裸机程序 七、重定位代码到DRAM【转】

    本文转载自:http://blog.csdn.net/eshing/article/details/37116637 一.关于DRAM 上一章我们讲解了如何对代码进行重定位,但是将代码重定位到只有25 ...

  8. android 中的Context(一)

    context的功能如此强大,它是activity的父类. public abstract class Context { ... public abstract Object getSystemSe ...

  9. PHP电影小爬虫(2)

    学习了别人的爬虫后自己改的一个,算是又回顾了一下php的使用 我们来利用simple_html_dom的采集数据实例,这是一个PHP的库,上手很容易.simple_html_dom 可以很好的帮助我们 ...

  10. [转]mysql触发器的作用及语法

    转自:http://blog.csdn.net/cloudday/article/details/6905590 触发器是一种特殊的存储过程,它在插入,删除或修改特定表中的数据时触发执行,它比数据库本 ...