bzoj4641 基因改造 KMP / hash
依稀记得,$NOIP$之前的我是如此的弱小....
完全不会$KMP$的写法,只会暴力$hash$....
大体思路为把一个串的哈希值拆成$26$个字母的位权
即$hash(S) = \sum\limits_{a} a * \sum w^i * [s[i] == a]$
通过记录每个字母第一次出现的位置,用$26$的时间来确定$f$是什么
然后通过确定的$f$计算出$f$是正确的时候的$hash$值,和原串的$hash$值比较
复杂度$O(26n)$
自然取模....
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
namespace remoon {
#define ri register int
#define ll long long
#define ull unsigned 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 --)
}
using namespace std;
using namespace remoon; #define sid 500050
char s[], t[sid]; int n, m, tim, f[];
int num[], vis[], nxt[], tot;
ull val[], wei[sid];
ull seed = ; inline void Init() {
wei[] = ;
rep(i, , m) wei[i] = wei[i - ] * seed;
rep(i, , m) {
int le = t[i];
if(!vis[le]) nxt[++ tot] = i, vis[le] = ;
val[le] += wei[m - i];
}
tim ++;
} inline void Solve() {
ull now = , tval = ;
rep(i, , m) now += s[i] * wei[m - i];
rep(i, m, n) { tim ++;
int flag = ; rep(j, , tot) {
int v = s[i - m + nxt[j]];
if(vis[v] == tim) { flag = ; break; }
if(vis[v] != tim) vis[v] = tim;
f[j] = v;
} if(!flag) {
tval = ;
rep(j, , tot)
tval += f[j] * val[t[nxt[j]]];
if(tval == now) write(i - m + );
} now -= s[i - m + ] * wei[m - ];
now *= seed; now += s[i + ];
}
} int main() {
scanf("%s", s + ); n = strlen(s + );
scanf("%s", t + ); m = strlen(t + );
Init(); Solve();
return ;
}
现在我明白了$KMP$是非常伟大的算法....
对于此题而言,考虑每个字符的上一个字符离当前字符的距离,这可以成为一个新串
然后比对新串即可
特别的,如果上一个字符出现的位置超过了匹配长度,那么我们也要视作合法
但是,我们发现这种匹配满足有前效性,没有后效性,因此可以用$KMP$
复杂度$O(n)$,十分的优秀
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
namespace remoon {
#define ri register int
#define rep(io, st, ed) for(ri io = st; io <= ed; io ++)
#define drep(io, ed, st) for(ri io = ed; io >= st; io --)
}
using namespace std;
using namespace remoon; const int sid = ; int n, m;
char s[sid], t[sid];
int lst[], S[sid], T[sid], nxt[sid]; inline bool match(int x, int y) {
if(y > m) return ;
if(x == T[y]) return ;
if(!T[y] && x >= y) return ;
return ;
} int main() {
scanf("%s", s + ); scanf("%s", t + );
n = strlen(s + ); m = strlen(t + ); rep(i, , n) {
S[i] = lst[s[i]] ? i - lst[s[i]] : ;
lst[s[i]] = i;
} memset(lst, , sizeof(lst));
rep(i, , m) {
T[i] = lst[t[i]] ? i - lst[t[i]] : ;
lst[t[i]] = i;
} for(ri i = , j = ; i <= m; i ++) {
while(j && !match(T[i], j + )) j = nxt[j];
if(match(T[i], j + )) j ++;
nxt[i] = j;
} for(ri i = , j = ; i <= n; i ++) {
while(j && !match(S[i], j + )) j = nxt[j];
if(match(S[i], j + )) j ++;
if(j == m) printf("%d\n", i - m + );
}
return ;
}
bzoj4641 基因改造 KMP / hash的更多相关文章
- BZOJ4641 基因改造[KMP]
这道题以前好像模拟的时候做过,当时不会做,于是用hash水过去了.. 正解是KMP,还是用当前字符与上一次相同字符位置的距离表示数组,于是数值相等时就代表相似.第一次出现用INF代替. 然后要匹配有多 ...
- 【BZOJ4641】基因改造 KMP
[BZOJ4641]基因改造 Description "人类智慧的冰峰,只有萌萌哒的我寂寞地守望." --TB TB正走在改造人类智慧基因的路上.TB发现人类智慧基因一点也不萌萌哒 ...
- Codeforces 1090J $kmp+hash+$二分
题意 给出两个字符串\(s\)和\(t\),设\(S\)为\(s\)的任意一个非空前缀,\(T\)为\(t\)的任意一个非空前缀,问\(S+T\)有多少种不同的可能. Solution 看了一圈,感觉 ...
- 【bzoj4641】基因改造 特殊匹配条件的KMP
题目描述 如果两个长度相等的字符串,如果存在一种字符的一一映射,使得第一个字符串的所有字符经过映射后与第二个字符串相同,那么就称它们“匹配”.现在给出两个串,求第一个字符串所有长度等于第二个字符串的长 ...
- 【BZOJ3940】【BZOJ3942】[Usaco2015 Feb]Censoring AC自动机/KMP/hash+栈
[BZOJ3942][Usaco2015 Feb]Censoring Description Farmer John has purchased a subscription to Good Hoov ...
- 【POJ2185】【KMP + HASH】Milking Grid
Description Every morning when they are milked, the Farmer John's cows form a rectangular grid that ...
- HDU 5782 Cycle(KMP+Hash)
[题目链接] http://acm.hdu.edu.cn/showproblem.php?pid=5782 [题目大意] 给出两个字符串,判断他们每一个前缀是否循环同构,循环同构的意思就是,字符串首位 ...
- Cycle (KMP + hash)
题意:给你2个串,让你判断2个字符串的前缀是否满足首尾连接形成的环是不是一样的. 思路:我们需要提前知道的是满足条件的前缀一定满足 strA = str1 + str2, strB = str2 + ...
- 7.26机房报零赛——无尽的矩阵【kmp+hash】
恩,其实大家都没有报零,反正我是蒟蒻 为了纪念我第一次打过哈希,特此写一篇题解 题目描述 从前有一个的小矩阵,矩阵的每个元素是一个字母(区分大小写),突然有一天它发生了 变异,覆盖了整个二维空间,即不 ...
随机推荐
- centos 安装memcache服务后memcahce本机连接Permission
自己手动在虚拟机下装了下memcache,整个过程真是充满波折,本身用php5.3安装memcache扩展就麻烦很多,无法通过yum直接安装,安装方法详见http://chenwei.me/blog/ ...
- npm 下载node-zookeeper包
环境:centos7(lunix) 1.安装nvm curl -o- https://raw.githubusercontent.com/creationix/nvm/v0.33.6/install. ...
- 13 JSON-RPC: a tale of interfaces
JSON-RPC: a tale of interfaces 27 April 2010 Here we present an example where Go's interfaces made i ...
- mysql自增id归0
mysql自增id归0 ALTER TABLE table_name AUTO_INCREMENT=1;
- Linux sudo 配置
1.配置sudo权限 sudo配置信息保存在 /etc/sudoers 文件中,可以使用vi修改.但使用vi修改需要给root用户设置权限,可以使用 visudo 直接修改. visudo 复制这行将 ...
- Activity总结练习
package com.example.wang.myapplication; import android.content.Intent; import android.os.Bundle; imp ...
- Adapter.notifyDataSetChanged()源码分析以及与ListView.setAdapter的区别
一直很好奇,notifyDataSetChanged究竟是重绘了整个ListView还是只重绘了被修改的那些Item,它与重新设置适配器即调用setAdapter的区别在哪里?所以特地追踪了一下源码, ...
- python使用cookie登陆网页
python2: Python 爬虫入门六之 Cookie 的使用 python3: Python3 网络爬虫 (六):Python3 使用 Cookie - 模拟登陆获取妹子联系方式 python ...
- loadrunner学习笔记之参数设置
一.关于参数的定义 在你录制程序运行的过程中,脚本生成器自动生成由函数组成的用户脚本.函数中参数的值就是在录制过程中输入的实际值. 例如,你录制了一个Web应用程序的脚本.脚本生成器生成了一个声明,该 ...
- spring boot get和post请求,以及requestbody为json串时候的处理
GET.POST方式提时, 根据request header Content-Type的值来判断: application/x-www-form-urlencoded, 可选(即非必须,因为这种情况的 ...