题目



输入格式

第一行两个整数N,M表示待检查的作文数量,和小强的标准作文库

的行数

接下来M行的01串,表示标准作文库

接下来N行的01串,表示N篇作文

输出格式

N行,每行一个整数,表示这篇作文的Lo 值。

输入样例

1 2

10110

000001110

1011001100

输出样例

4

提示

输入文件不超过1100000字节

注意:题目有改动,可识别的长度不小于90%即可,而不是大于90%

题解

想来练练SAM,却跪在了单调队列DP上。。。QAQ

根据后缀数组进行多串匹配时,用一个未出现的字符将各串连接起来形成一个串

我们可以照搬到后缀自动机来,对m个模板串建立后缀自动机

此时我们可以在后缀自动机上走一遍求出匹配串每个位置为结尾所能匹配的最大长度【见代码init】

记为len[i],即表示匹配的子串\([i - len[i] + 1,i]\)可以与模板串匹配

有了这些,我们如何计算\(L_0\)?

仔细思考发现L具有单调性,即在长度至少为L时选择的最小切割方案下,对长度小于L的限制仍然满足,长度小于L的可以照搬L的切割方法而满足条件

这个时候我们就可以二分了

对于长度L,我们可以用DP检验,令f[i]表示前i个字符在L限制下能匹配的最大长度

首先f[i]至少等于f[i - 1],因为i可以单独被化为出来而转移到f[i - 1]【也就是说f[i]单调不下降,这个待会会用到】

我们有状态转移方程:\(f[i] = max{f[j] + i - j},j\in [i - len[i],i - L]\)

解释:由于L的限制,要切割出以i结尾产生贡献的区间,必须在i - L及之前,而i最多匹配到i - len[i] + 1,所以\(j\in [i - len[i],i - L]\)

最后的匹配长度就是i割出这段的长度 + f[j]

这样做是\(O(n^2logn)\)的,不能满足题意

考虑优化

我们将方程变形\(f[i] - i = f[j] - j\),左边是只与i有关的量,与决策无关

右边是只与j有关的量,由f[j]的定义,f[j] - j表示前j个没匹配的个数【负数】,由于没匹配的个数不会变少,所以f[j] - j单调递减

也就是说决策具有单调性,且只与j有关,可以采用单调队列优化

【一定要好好写队列QAQ,WA了几发】

呼啦啦~~~搞完了【CTSC的题目哇。】

#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long int
#define REP(i,n) for (int i = 1; i <= (n); i++)
#define Redge(u) for (int k = h[u],to; k; k = ed[k].nxt)
#define BUG(s,n) for (int i = 1; i <= (n); i++) cout<<s[i]<<' '; puts("");
using namespace std;
const int maxn = 3000005,maxm = 100005,INF = 1000000000;
int pre[maxn],ch[maxn][3],step[maxn],cnt,last;
char s[maxn];
void ins(int x){
int p = last,np = ++cnt;
last = np; step[np] = step[p] + 1;
while (p && !ch[p][x]) ch[p][x] = np,p = pre[p];
if (!p) pre[np] = 1;
else {
int q = ch[p][x];
if (step[q] == step[p] + 1) pre[np] = q;
else {
int nq = ++cnt; step[nq] = step[p] + 1;
memcpy(ch[nq],ch[q],sizeof(ch[q]));
pre[nq] = pre[q]; pre[q] = pre[np] = nq;
while (ch[p][x] == q) ch[p][x] = nq,p = pre[p];
}
}
}
int q[maxn],head,tail,f[maxn],N,len[maxn];
bool check(int L){
head = tail = f[0] = 0;
for (int i = 1; i <= N; i++){
f[i] = f[i - 1];
int l = i - len[i],r = i - L;
if (r >= 0){
while (head < tail && f[q[tail - 1]] - q[tail - 1] < f[r] - r) tail--;
q[tail++] = r;
}
while (head < tail && q[head] < l) head++;
if (head < tail) f[i] = max(f[i],f[q[head]] - q[head] + i);
}
return 10 * f[N] >= 9 * N;
}
void init(){
int u = 1,id,ans = 0;
for (int i = 1; i <= N; i++){
id = s[i] - '0';
if (ch[u][id]) len[i] = ++ans,u = ch[u][id];
else {
while (u != 1 && !ch[u][id]) u = pre[u];
if (u == 1) len[i] = ans = 0;
else len[i] = ans = step[u] + 1,u = ch[u][id];
}
}
}
int main(){
int n,m; last = cnt = 1;
scanf("%d%d",&n,&m);
while (m--){
scanf("%s",s); N = strlen(s);
for (int i = 0; i < N; i++) ins(s[i] - '0');
ins(2);
}
while (n--){
scanf("%s",s + 1); N = strlen(s + 1);
init();
int l = 0,r = N,mid,ans = 0;
while (l <= r){
mid = l + r >> 1;
if (check(mid)) l = mid + 1,ans = mid;
else r = mid - 1;
}
printf("%d\n",ans);
}
return 0;
}

BZOJ2806 [Ctsc2012]Cheat 【后缀自动机 + 二分 + 单调队列优化DP】的更多相关文章

  1. BZOJ 2806 [Ctsc2012]Cheat (后缀自动机+二分+单调队列+dp)

    题目大意: 给你一堆模式串和文本串 对于每个文本串,我们可以把它不可重叠地拆分成很多子串,如果拆分出的串作为子串出现在了任何一个模式串中,我们称它是“眼熟的”,我们必须保证“眼熟的”子串总长度不小于文 ...

  2. 【BZOJ2806】[Ctsc2012]Cheat 广义后缀自动机+二分+单调队列优化DP

    [BZOJ2806][Ctsc2012]Cheat Description Input 第一行两个整数N,M表示待检查的作文数量,和小强的标准作文库的行数接下来M行的01串,表示标准作文库接下来N行的 ...

  3. BZOJ2806: [Ctsc2012]Cheat(广义后缀自动机,单调队列优化Dp)

    Description Input 第一行两个整数N,M表示待检查的作文数量,和小强的标准作文库的行数接下来M行的01串,表示标准作文库接下来N行的01串,表示N篇作文 Output N行,每行一个整 ...

  4. [BZOJ1044][HAOI2008]木棍分割 二分 + 单调队列优化dp + 滚动数组优化dp

    Description 有n根木棍, 第i根木棍的长度为Li,n根木棍依次连结了一起, 总共有n-1个连接处. 现在允许你最多砍断m个连接处, 砍完后n根木棍被分成了很多段,要求满足总长度最大的一段长 ...

  5. 2018.09.26洛谷P3957 跳房子(二分+单调队列优化dp)

    传送门 表示去年考普及组的时候失了智,现在看来并不是很难啊. 直接二分答案然后单调队列优化dp检验就行了. 注意入队和出队的条件. 代码: #include<bits/stdc++.h> ...

  6. [bzoj2806][Ctsc2012]Cheat(后缀自动机(SAM)+二分答案+单调队列优化dp)

    偷懒直接把bzoj的网页内容ctrlcv过来了 2806: [Ctsc2012]Cheat Time Limit: 20 Sec  Memory Limit: 256 MBSubmit: 1943   ...

  7. 单调队列优化DP,多重背包

    单调队列优化DP:http://www.cnblogs.com/ka200812/archive/2012/07/11/2585950.html 单调队列优化多重背包:http://blog.csdn ...

  8. bzoj1855: [Scoi2010]股票交易--单调队列优化DP

    单调队列优化DP的模板题 不难列出DP方程: 对于买入的情况 由于dp[i][j]=max{dp[i-w-1][k]+k*Ap[i]-j*Ap[i]} AP[i]*j是固定的,在队列中维护dp[i-w ...

  9. hdu3401:单调队列优化dp

    第一个单调队列优化dp 写了半天,最后初始化搞错了还一直wa.. 题目大意: 炒股,总共 t 天,每天可以买入na[i]股,卖出nb[i]股,价钱分别为pa[i]和pb[i],最大同时拥有p股 且一次 ...

随机推荐

  1. python_66_生成器2

    import time def consumer(name): print('%s准备吃包子 '%name) while True: baozi=yield print('包子[%s]来了,被[%s] ...

  2. swl模拟

    #import <Foundation/Foundation.h> #define NSLog(FORMAT, ...) printf("%s\n", [[NSStri ...

  3. Java 技术栈

    JAVA是一个面向对象的编程语言,由SUN公司的程序员所开发.它不仅吸收了C++的各种优点,而且还撇弃了C++中难以理解的概念,如多继承.指针等:因此JAVA语言具有功能强大且简单易用两个特征, JA ...

  4. 极路由安装SS,SSR,搬运,侵权删除

    打开路由器的开发者模式 使用SSH,登录到路由器后台 复制以下代码,粘贴即可 cd /tmp $ echo "127.0.0.1 hiwifiss.ml" > /etc/ho ...

  5. php实现的三个常用加密解密功能函数示例

    目录 算法一: 算法二: 算法三(改进第一个加密之后的算法) 本文实例讲述了php实现的三个常用加密解密功能函数.分享给大家供大家参考,具体如下: 算法一: //加密函数 function lock_ ...

  6. 用Python学分析 - 单因素方差分析

    单因素方差分析(One-Way Analysis of Variance) 判断控制变量是否对观测变量产生了显著影响 分析步骤 1. 建立检验假设 - H0:不同因子水平间的均值无差异 - H1:不同 ...

  7. Python 列表元素分组,比如 [1,2,3,...20]变成 [[1,2,3],[4,5,6]....](列表生成式解决)

    # 生成一个1到20的列表 a=[x for x in range(1,21)] # 把a列表切片,并赋值给b列表,x为0到20且步长为3的列表,在这里具体为[0,3,6,9,12,15,18] # ...

  8. 裸奔着造房子——对政府禁止采购Win8系统的一些看法

    前段时间有消息称政府招标的项目将禁止使用Win8系统,原因是Win8系统的安全架构将有利于暴露敏感信息给微软,而微软的老子是美利坚,老子想要知道什么,儿子当然不敢不从.因此Win8也被打入冷宫,微软多 ...

  9. Android开发——用户在屏幕上的手势识别

    个定点决定.四个属性分别为left(1),top(2),right(3),bottom(4). 数字为图上标出的距离.显然这四个属性是相对于父容器来定的,均可以通过get()方法获取. 因此很容易得出 ...

  10. 4 Template层-验证码

    1.验证码 在用户注册.登录页面,为了防止暴力请求,可以加入验证码功能,如果验证码错误,则不需要继续处理,可以减轻一些服务器的压力 使用验证码也是一种有效的防止crsf的方法 验证码效果如下图: 官网 ...