URAL 1996 Cipher Message 3 (FFT + KMP)
转载请注明出处,谢谢http://blog.csdn.net/ACM_cxlove?viewmode=contents by---cxlove
题意 :给出两个串A , B,每个串是若干个byte,A串的每个byte的最后一个bit是可以修改的。问最少修改多少,使得B串是A的一个子串。
2013年NEERC的题。。。。。。。感觉[buaa]sd0061教我做这题。
NEERC是毛子题,但是这套题感觉除了题面很难读之外,并不是很难。。。
目前为止,J题貌似全队都没看题,E题还在WA中,其它题都没啥问题。
做法:前7位都是不能修改的,所以要完全匹配,那么先按前7位KMP一下,记录匹配的位置。之后就是统计以这些位置开始,需要修改多少位,然后 取最小值。
所以提取出最后一个bit,成了两个01串,统计hamming distance。
记得若干年前就和队友讨论过这个问题,不过当时不是01串,觉得不可在n * m之内解决,然后 就放弃了。结果这个01串是可以 在nlgn解决 的。吓傻。。
两个01串a , b。将b反转之后,求一次卷积,便可以得到a串中以i为起始位置,与b进行匹配有多少个位置同为1。
那么接下来把a,b的01反转一下,再求一次卷积,就可以得到原串中为多少个位置同为0。就得到了hamming距离。
卷积求的是c[i + j] = a[i] * b[j]。由于我们将b进行了反转,所以a串中以i为起始位置的。
c[i + m - 1] = a[i + 0] * b[m - 0 - 1] + a[i + 1] * b[m - 1 - 1] + …… a[i + j] * b[m - j - 1] + …… a[i + m - 1] * b[m - (m - 1) - 1]。
如果a , b中同为1,乘积为1,否则为0,这样就统计出了有多少位同为1了。吓傻。。。太神了。
预处理出这个之后,再枚举之前KMP所得到的匹配位置 就可以了。
#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <vector>
#include <algorithm>
#include <queue>
#include <set>
using namespace std;
typedef long long LL;
typedef pair<int , int> pii;
const int N = 1000005;
//FFT copy from kuangbin
const double pi = acos (-1.0);
// Complex z = a + b * i
struct Complex {
double a, b;
Complex(double _a=0.0,double _b=0.0):a(_a),b(_b){}
Complex operator + (const Complex &c) const {
return Complex(a + c.a , b + c.b);
}
Complex operator - (const Complex &c) const {
return Complex(a - c.a , b - c.b);
}
Complex operator * (const Complex &c) const {
return Complex(a * c.a - b * c.b , a * c.b + b * c.a);
}
};
//len = 2 ^ k
inline void change (Complex y[] , int len) {
for (int i = 1 , j = len / 2 ; i < len -1 ; i ++) {
if (i < j) swap(y[i] , y[j]);
int k = len / 2;
while (j >= k) {
j -= k;
k /= 2;
}
if(j < k) j += k;
}
}
// FFT
// len = 2 ^ k
// on = 1 DFT on = -1 IDFT
inline void FFT (Complex y[], int len , int on) {
change (y , len);
for (int h = 2 ; h <= len ; h <<= 1) {
Complex wn(cos (-on * 2 * pi / h), sin (-on * 2 * pi / h));
for (int j = 0 ; j < len ; j += h) {
Complex w(1 , 0);
for (int k = j ; k < j + h / 2 ; k ++) {
Complex u = y[k];
Complex t = w * y [k + h / 2];
y[k] = u + t;
y[k + h / 2] = u - t;
w = w * wn;
}
}
}
if (on == -1) {
for (int i = 0 ; i < len ; i ++) {
y[i].a /= len;
}
}
}
int n , m , a[N] , b[N];
int c[N] , d[N] , next[N] , p[N];
char str[10];
void get_next (int *a , int l) {
int i = 0 , j = -1;
next[i] = j;
while (i < l) {
if (j == -1 || a[i] == a[j]) {
i ++; j ++;
next[i] = j;
}
else j = next[j];
}
}
int pos[N] , cnt;
void match (int *a , int la , int *b , int lb) {
int i = 0 , j = 0;
while (i < la) {
if (j == -1 || a[i] == b[j]) {
i ++;j ++;
if (j == lb) {
pos[cnt ++] = i - lb;
j = next[j];
}
}
else j = next[j];
}
}
Complex x1[N] , x2[N];
void gao (int *a , int *b) {
int len = max (n , m);
int l = 1;
while (l < len * 2) l <<= 1;
for (int i = 0 ; i < n ; i ++)
x1[i] = Complex (a[i] , 0);
for (int i = n ; i < l ; i ++)
x1[i] = Complex (0 , 0);
FFT (x1 , l , 1);
for (int i = 0 ; i < m ; i ++)
x2[i] = Complex (b[i] , 0);
for (int i = m ; i < l ; i ++)
x2[i] = Complex (0 , 0);
FFT (x2 , l , 1);
for (int i = 0 ; i < l ; i ++)
x1[i] = x1[i] * x2[i];
FFT (x1 , l , -1);
for (int i = 0 ; i <= n - m ; i ++) {
p[i] += (int)(x1[i + m - 1].a + 0.5);
}
}
int main () {
while (scanf ("%d %d" , &n , &m) != EOF) {
for (int i = 0 ; i < n ; i ++) {
scanf ("%s" , str);
a[i] = 0;
for (int j = 0 ; j < 7 ; j ++)
a[i] = a[i] * 2 + str[j] - '0';
c[i] = str[7] - '0';
}
for (int i = 0 ; i < m ; i ++) {
scanf ("%s" , str);
b[i] = 0;
for (int j = 0 ; j < 7 ; j ++)
b[i] = b[i] * 2 + str[j] - '0';
d[m - i - 1] = str[7] - '0';
}
get_next (b , m);
cnt = 0;
match (a , n , b , m);
if (cnt == 0) puts ("No");
else {
memset (p , 0 , sizeof(p));
puts ("Yes");
gao (c , d);
for (int i = 0 ; i < n ; i ++)
c[i] = c[i] ^ 1;
for (int i = 0 ; i < m ; i ++)
d[i] = d[i] ^ 1;
gao (c , d);
int ans = m + 1 , idx = -1;
for (int i = 0 ; i < cnt ; i ++) {
if (m - p[pos[i]] < ans)
ans = m - p[pos[i]] , idx = pos[i];
}
printf ("%d %d\n" , ans , idx + 1);
}
}
return 0;
}
URAL 1996 Cipher Message 3 (FFT + KMP)的更多相关文章
- URAL 1996 Cipher Message 3
题目 神题. 记得当初DYF和HZA讲过一个FFT+KMP的题目,一直觉得很神,从来没去做. 没有真正理解FFT的卷积. 首先考虑暴力. 只考虑前7位 KMP 找出所有 B 串可以匹配 A 串的位置. ...
- URAL 1996. Cipher Message 3(KMP+fft)
传送门 解题思路 因为要完全匹配,所以前七位必须保证相同,那么就可以把前7位提出来做一遍\(kmp\)匹配,最后的答案一定在这些位置里.考虑最后一位,可以把最后一位单独取出来,要计算的是最后一位相同的 ...
- Ural 1996 Cipher Message 3 (生成函数+FFT)
题面传送门 题目大意:给你两个$01$串$a$和$b$,每$8$个字符为$1$组,每组的最后一个字符可以在$01$之间转换,求$b$成为$a$的一个子串所需的最少转换次数,以及此时是从哪开始匹配的. ...
- URAL1996 Cipher Message 3(KMP + FFT)
题目 Source http://acm.timus.ru/problem.aspx?space=1&num=1996 Description Emperor Palpatine has be ...
- URAL 1654 Cipher Message 解题报告
题目链接:http://acm.timus.ru/problem.aspx?space=1&num=1654 题意:简单的理解就是,把一个序列中相邻的且是偶数个相同的字符删除,奇数个的话就只保 ...
- ural Cipher Message
Cipher Message Time Limit:1000MS Memory Limit:65536KB 64bit IO Format:%I64d & %I64u Desc ...
- Gym 100285G Cipher Message 3
题意 给\(N,M(N,M \le 250000)\)的两个由8位二进制表示的两个序列,允许改变每个数字的第8位的数值(即0→1,1→0),求改变最少次数使得长为\(M\)的序列为长为\(N\)的连续 ...
- 北京集训TEST16——图片加密(fft+kmp)
题目: Description CJB天天要跟妹子聊天,可是他对微信的加密算法表示担心:“微信这种加密算法,早就过时了,我发明的加密算法早已风靡全球,安全性天下第一!” CJB是这样加密的:设CJB想 ...
- HDU4300 Clairewd’s message(拓展kmp)
Problem Description Clairewd is a member of FBI. After several years concealing in BUPT, she interce ...
随机推荐
- 小米开源监控open-falcon
小米开源监控系统Open-Falcon安装使用笔记 07net01.com 发布于 2016-10-25 18:42:03 分类:IT技术 阅读(88) 评论 前言 近期爆出Zabbix有严重bug, ...
- virtualenv基本使用
win 安装 virtualenv pip3 install virtualenv 创建虚拟环境 virtualenv env1 进入虚拟环境 env1/Scripts/activate 退出虚拟环境 ...
- 善待Erlang 代码 -- Xref 实践
Xref 是一个交叉引用工具,通过分析定义的函数间的调用关系,用于查找函数. 模块. 应用程序和版本之间的依赖关系. 通俗而言,Xref 可以检查代码中函数的调用关系.比如在 moduleA 中的 f ...
- 腾讯互娱开源分布式开发框架Pebble
构建游戏世界的Pebble 愿景:出色的游戏服务器端底层框架 现代游戏项目中,为了让更多的玩家能在一起玩,游戏服务器所需要承载的在线玩家数量越来越多.同时为了让游戏更好玩,越来越多复杂的业务逻辑都 ...
- Java File文件操作 创建文件\目录,删除文件\目录
Java手册 java.io 类 File java.lang.Object java.io.File 所有已实现的接口: Serializable, Comparable<File> p ...
- Hive-表连接
Hive只支持等值连接,即ON子句中使用等号连接,不支持非等值连接. Hive内置的数据存储类型,TextFile, SequenceFile, ORC(列式存储) 如果连接语句中有WHERE子句,会 ...
- CFGym 101505I 题解
一.题目链接 http://codeforces.com/gym/101505 二.题意 这题其实主要就是题意,理解题意后,就是水题了.我想了下,主要原因就是这几点: 1.题意太过英文化,很多句子不能 ...
- Python Web框架——Flask
简介 Flask是一个基于Python开发并且依赖jinja2模板和Werkzeug WSGI服务的一个微型框架,对于Werkzeug本质是Socket服务端,其用于接收http请求并对请求进行预处理 ...
- Shell 函数库
1.为什么要定义函数库 经常使用的重复代码封装成函数文件 一般不直接执行,而是由其他脚本调用 2.编写一个函数库,该函数库实现以下几个函数. 1.加法函数:add 2.减法函数:reduce 3.乘法 ...
- python 2.7中安装mysql
在python中进行安装mysql模块,但是怎么都不能导入mysql模块,出错如下所示: [root@python ~]# python Python 2.7.11 (default, Apr 5 2 ...