luogu4770 [NOI2018]你的名字 后缀自动机 + 线段树合并

1326433-20181127215714033-72785805.png


其实很水的一道题吧....

题意是:每次给定一个串\(T\)以及\(l, r\),询问有多少个字符串\(s\)满足,\(s\)\(T\)的子串,但不是\(S[l .. r]\)的子串


统计\(T\)本质不同的串,建个后缀自动机

然后自然的可以想到,对于每个\(T\)的子串,它对应了一个\(right\)集合

那么,它应该会被这个\(right\)集合所限制

考虑对于每个\(i\),求出最小的\(l\)使得\(T[l .. i]\)存在于\(S[l..r]\)

这个可以套个线段树转移

然后就没了.....


如果不需要统计\(T\)本质不同的串,又怎么做呢?

统计的时候乘上\(right\)集合大小就行


#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;#define ri register int
#define ll long long
#define rep(io, st, ed) for(ri io = st; io <= ed; io ++)
#define drep(io, ed, st) for(ri io = ed; io >= st; io --)#define gc getchar
inline int read() {int p = 0, w = 1; char c =  gc();while(c > '9' || c < '0') { if(c == '-') w = -1; c = gc(); }while(c >= '0' && c <= '9') p = p * 10 + c - '0', c = gc();return p * w;
}const int sid = 1005000;
const int eid = 30000000 + 5;struct SAM {int id, fa[sid], mx[sid];int go[sid][26], mc[sid];inline int newnode() {++ id;fa[id] = mx[id] = mc[id] = 0;memset(go[id], 0, sizeof(go[id]));return id;}inline void init() {id = 0;newnode();}inline int extend(int lst, int c, int pos) {int np = newnode(), p = lst;mx[np] = mx[p] + 1; mc[np] = pos;for( ; p && !go[p][c]; p = fa[p]) go[p][c] = np;if(!p) fa[np] = 1;else {int q = go[p][c];if(mx[p] + 1 == mx[q]) fa[np] = q;else {int nq = newnode(); mx[nq] = mx[p] + 1;fa[nq] = fa[q]; fa[np] = fa[q] = nq;memcpy(go[nq], go[q], sizeof(go[q]));for( ; p && go[p][c] == q; p = fa[p]) go[p][c] = nq;}}return np;}} S, T;int q, n, m, seg;
char s[sid], t[sid];
int nc[sid], ip[sid], w[sid], val[sid];
int rt[sid], ls[eid], rs[eid];inline int merge(int x, int y) {if(!x || !y) return x + y;int o = ++ seg;ls[o] = merge(ls[x], ls[y]);rs[o] = merge(rs[x], rs[y]);return o;
}inline void ins(int &o, int l, int r, int p) {o = ++ seg;if(l == r) return;int mid = (l + r) >> 1;if(p <= mid) ins(ls[o], l, mid, p);else ins(rs[o], mid + 1, r, p);
}inline bool qry(int o, int l, int r, int ml, int mr) {if(ml > r || mr < l || ml > mr || !o) return 0;if(ml <= l && mr >= r) return 1;int mid = (l + r) >> 1;if(qry(ls[o], l, mid, ml, mr)) return 1;else return qry(rs[o], mid + 1, r, ml, mr);
}inline void init() {S.init();int lst = 1;rep(i, 1, n) lst = S.extend(lst, s[i] - 'a', i);int id = S.id;rep(i, 1, id) nc[S.mx[i]] ++;rep(i, 1, n) nc[i] += nc[i - 1];rep(i, 1, id) ip[nc[S.mx[i]] --] = i;rep(i, 1, id) if(S.mc[i]) ins(rt[i], 1, n, S.mc[i]);drep(i, id, 1) {int o = ip[i], f = S.fa[o];rt[f] = merge(rt[f], rt[o]);}
}void Match(int l, int r) {int o = 1, nl = 0;rep(i, 1, m) {int c = t[i] - 'a';while(1) {int nxt = S.go[o][c], f = S.fa[o];if(nxt && qry(rt[nxt], 1, n, l + nl, r)) {nl ++; o = nxt;break;}if(!nl) break; nl --;if(nl == S.mx[f]) o = f;}w[i] = nl;}
}int main() {scanf("%s", s + 1);n = strlen(s + 1);init(); q = read();rep(i, 1, q) {scanf("%s", t + 1);m = strlen(t + 1);T.init();int lst = 1;rep(j, 1, m) lst = T.extend(lst, t[j] - 'a', j);int l = read(), r = read();Match(l, r);int id = T.id;rep(i, 1, id) nc[i] = val[i] = 0;rep(i, 1, id) nc[T.mx[i]] ++;rep(i, 1, id) nc[i] += nc[i - 1];rep(i, 1, id) ip[nc[T.mx[i]] --] = i;drep(i, id, 1) {int o = ip[i], f = T.fa[o];if(T.mc[o]) val[o] = w[T.mc[o]];val[f] = max(val[f], val[o]);}ll ans = 0;rep(i, 1, id) ans += max(T.mx[i] - max(T.mx[T.fa[i]], val[i]), 0);printf("%lld\n", ans);}return 0;
}

转载于:https://www.cnblogs.com/reverymoon/p/10029355.html

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/364796.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

centos-7.2 node.js免编译安装

cd /usr/local/wget https://npm.taobao.org/mirrors/node/v8.9.3/node-v8.9.3-linux-x64.tar.gz tar -zxvf node-v8.9.3-linux-x64.tar.gz //已编译可以直接运行./bin下面的命令rm -rf node-v8.9.3-linux-x64.tar.gz //解压完删包 // 建立全局快捷方式 ln -s 源命令文件 快…

团队计划会议

跟航哥想了挺多要做什么&#xff0c;要完成什么&#xff0c;以什么为主要功能 提出了几个想法&#xff0c;并做了投票 最后决定一起做一个跑腿软件 最初的任务量&#xff1a; 跟航哥商量两个人一人负责两个模块 航哥负责管理员和下单 我负责接单跟其他琐碎的小功能 呐&#xff…

在基于Spring MVC的应用程序中配置favicon.ico

Favicon是与您的网站相关的图标&#xff08;favicon.ico&#xff09;。 并非每个网站都在使用favicon。 但是大多数浏览器并不关心它&#xff0c;反正他们都要求它。 当图标图标不在适当位置时&#xff0c;服务器将返回不必要的404 Not Found错误。 在典型的Spring MVC应用程序…

vue.js java php_准吗?Java程序员喜欢AngularJS,PHP程序员喜欢Vue.js!

编程语言与框架或者库之间有联系是很正常的事情&#xff0c;如果我们告诉你&#xff0c;使用某一种编程语言或技术的开发人员可能更喜欢某个框架&#xff0c;你会作何反应呢&#xff1f;Stack Overflow根据网站内最常访问的标签将开发人员分为多个组&#xff0c;并检查了每组每…

BAJT高级Java面试题

答对这些面试题&#xff0c;PASS 掉 80 % 的竞争者 hashcode相等两个类一定相等吗?equals呢?相反呢? 介绍一下集合框架? hashmap hastable 底层实现什么区别?hashtable和concurrenthashtable呢? hashmap和treemap什么区别?低层数据结构是什么? 线程池用过吗都有什么…

结构化日志:出错时你最想要的好朋友

目录 介绍什么是日志&#xff1f;Grab中日志的状况为什么改变&#xff1f;结构化日志支持不同格式的多写开发中类似生产环境的日志因果顺序但为什么要结构化记日志&#xff1f;原文&#xff1a;Structured Logging: The Best Friend You’ll Want When Things Go Wrong 介绍 在…

在vue项目中添加特殊字体

这里的特殊字体&#xff0c;指的是一般用户电脑未安装到本地的字体&#xff0c;要引入这样的字体&#xff0c;首先需要把字体文件下载下来。 就像上图这样的&#xff0c;ttf格式的&#xff0c;然后在项目里添加它。 然后我们在font.css里用font-face规则引入这个字体文件并命名…

使用Spring WS创建合同优先的Web服务

1引言 本文介绍了如何使用来实现和测试SOAP Web服务 Spring Web Services项目 。 本示例使用JAXB2进行&#xff08;取消&#xff09;编组。 为了开发服务&#xff0c;我将使用合同优先的方法&#xff0c;该方法首先定义服务合同&#xff0c;然后基于该合同实施服务。 本文分为…

java中的常用日期类_Java中的常用日期类说明

日期类常用的有三个&#xff0c;Date类&#xff0c;Calendar(日历)类和日期格式转换类(DateFormat)Date类中的大部分的方法都已经过时&#xff0c;一般只会用到构造方法取得系统当前的时间。public class DateDemo {public static void main(String[] args) {Date date new Da…

转载 Net多线程编程—System.Threading.Tasks.Parallel

.Net多线程编程—System.Threading.Tasks.Parallel System.Threading.Tasks.Parallel类提供了Parallel.Invoke&#xff0c;Parallel.For&#xff0c;Parallel.ForEach这三个静态方法。 1 Parallel.Invoke 尽可能并行执行所提供的每个操作&#xff0c;除非用户取消了操作。 方法…

三方面搞定http协议之“状态码”

当我们向服务器请求数据的时候&#xff0c;服务器会给我们一个反馈&#xff0c;告诉我们对待我们的请求&#xff0c;服务器处理得怎么样了&#xff0c;而这个反馈&#xff0c;是通过数字来传达的&#xff0c;这个数字就叫状态码。 状态码分为以下几种&#xff1a; 1xx&#xf…

哪个更好的选择:克隆或复制构造函数?

这就是我开始撰写本文的方式。 我已经读过很多次这样的声明&#xff1a; “当对象引用可变的最终字段时&#xff0c;克隆变得很困难。” 每次我在Google上搜索它时&#xff0c;都要了解它的确切含义&#xff0c;并且在此过程中也忘了它。 因此以为我会在此撰写博客&#xff0c;…

Fiddler教程--简介

1、开发环境host配置自己修改系统的host来回挺麻烦的 2、前后的接口调试 3、线上bugfix 4、性能分析和优化 5.等等... 工作原理 一个代理服务器地址改为 127.0.0.1:8888流模式边走边返回缓冲模式http请求完成所有的数据之后&#xff0c;才返回 界面功能介绍 1.工具栏 从下图红色…

java map是有序的吗_Java:如何初始化和填充最终的静态有序Map?

我在Java中有一个词干算法,它需要一个静态的最终HashMap< String,String>预先填写了大约30 000条记录.我需要地图按照插入的顺序保存记录(我得到一个提示,我可以使用LinkedHashMap&#xff1f;).我以为我可以在Java类文件中手动插入值,因为这是在RAM中加载它们的最快方法…

跟面向对象卯上了,看看ES6的“类”

上回我们说到ES5的面向对象&#xff0c;以及被大家公认的最佳的寄生组合式继承。时代在进步&#xff0c;在ES6中对于面向对象这个大boss理所应当地进行了一次大改&#xff0c;从原先那种比较长的写法转变为“小清新”写法。我们一起来看一下。 在ES6中是有类这个概念&#xff0…

js 变量作用域

例子 <script>var a "heh"function findLove(){console.log(a);function findforyou(){var a "you";console.log(a);}function findother(){console.log(a)}findforyou();findother();}findLove(); </script> 输出 heh you heh 例子 <scri…

Jin Ge Jin Qu hao UVA - 12563 01背包

题目&#xff1a;题目链接 思路&#xff1a;由于t最大值其实只有180 * 50 678&#xff0c;可以直接当成01背包来做&#xff0c;需要考虑的量有两个&#xff0c;时间和歌曲数&#xff0c;其中歌曲优先级大于时间&#xff0c;于是我们将歌曲数作为背包收益&#xff0c;用时间作为…

Java 8中的5个功能将改变您的编码方式

Java 8在JVM和语言级别都包含了一些非常令人兴奋的功能。 虽然最初为该发行版设想的某些功能已扩大范围或已推出到第9版&#xff0c;但实际上有数十个新功能。 许多新添加的内容在编译器&#xff0c;JVM或帮助系统级别都进行了后台改进。 这样&#xff0c;虽然我们可能会从中受…

Java相关资料分享(视频+电子书籍)

关注微信公众号【Java典籍】&#xff0c;获取百度网盘提取码 ▼微信扫一扫下图↓↓↓二维码关注 转载于:https://www.cnblogs.com/bingyimeiling/p/10279049.html

vue项目 一行js代码搞定点击图片放大缩小

一行js代码搞定xue项目需要点击图片放大缩小&#xff0c;其实主要用的是用到了vue:class的动态切换&#xff0c;内容比较简单。一开始我把维护的需求想得太复杂了&#xff0c;和测试小姐姐聊了一下才反应过来。 两个月不到跟了四个项目&#xff0c;现在是维护改bug阶段&#x…