博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
bzoj 4566 [Haoi2016]找相同字符——广义后缀自动机
阅读量:6235 次
发布时间:2019-06-22

本文共 1177 字,大约阅读时间需要 3 分钟。

题目:

每个后缀结尾处 ct[ ] = 1 ,按拓扑序 dp 一下就能求出 right 集合的大小。自动机上每个点的贡献就是 ( l [cr]-l [fa] ) * ct[0][cr] * ct[1][cr] , ct[0] 和 ct[1] 表示在两个字符串里分别的出现次数。

#include
#include
#include
#define ll long longusing namespace std;const int N=2e5+5,M=8e5+5,K=30;int tot=1,go[M][K],fa[M],l[M],ct[2][M],ans[M];char s[2][N]; int q[M],tx[N],vis[M];int cz(int p,int w){ int q=go[p][w],nq=++tot;l[nq]=l[p]+1; fa[nq]=fa[q];fa[q]=nq; memcpy(go[nq],go[q],sizeof go[q]); for(;p&&go[p][w]==q;p=fa[p])go[p][w]=nq; return nq;}int ins(int p,int w){ if(go[p][w]) { int q=go[p][w]; if(l[q]==l[p]+1)return q; return cz(p,w); } int np=++tot;l[np]=l[p]+1; for(;p&&!go[p][w];p=fa[p])go[p][w]=np; if(!p)fa[np]=1; else { int q=go[p][w]; if(l[q]==l[p]+1)fa[np]=q; else fa[np]=cz(p,w); } return np;}void Rsort(int n){ for(int i=1;i<=tot;i++)tx[l[i]]++; for(int i=1;i<=n;i++)tx[i]+=tx[i-1]; for(int i=1;i<=tot;i++)q[tx[l[i]]--]=i;}int main(){ scanf("%s",s[0]);scanf("%s",s[1]); int n[2];n[0]=strlen(s[0]);n[1]=strlen(s[1]); for(int t=0;t<=1;t++) for(int pr=1,i=0;i

 

转载于:https://www.cnblogs.com/Narh/p/10291156.html

你可能感兴趣的文章
02 回归算法 - 线性回归求解 θ(最大似然估计求解)
查看>>
【算法学习】算法的复杂度
查看>>
14、python常用模块
查看>>
swift4.0 CAKeyframeAnimation,抖动动画
查看>>
比特币解锁脚本中的ScriptSignature都包含了什么东西
查看>>
从零开始学java(猜数字游戏)
查看>>
xcode10 编译报错Multiple commands produce
查看>>
使用ab进行页面的压力测试
查看>>
开源大数据周刊-2018年08月03日 第95期
查看>>
React, a gentle introduction.
查看>>
值得一试的Android开发工具
查看>>
如何解决Android开发学习过程中缺乏后端接口的问题「Android,资源向」
查看>>
关闭wps自动升级
查看>>
【设计模式】面向对象六大原则
查看>>
Web 通信 之 长连接、长轮询(long polling)
查看>>
Git教程摘录
查看>>
JQuery基本知识框架思维导图(上)
查看>>
Java 数据类型
查看>>
iView 发布微信小程序 UI 组件库 iView Weapp
查看>>
运维Caron 的一条心理的os
查看>>