SPOJ MAXMATCH - Maximum Self-Matching (FFT)
题目链接:MAXMATCH - Maximum Self-Matching
Description
You're given a string s consisting of letters 'a', 'b' and 'c'.
The matching function \(m_s( i )\) is defined as the number of matching characters of s and its i-shift. In other words, \(m_s( i )\) is the number of characters that are matched when you align the 0-th character of s with the i-th character of its copy.
You are asked to compute the maximum of \(m_s( i )\) for all i ( 1 <= i <= |s| ). To make it a bit harder, you should also output all the optimal i's in increasing order.
Input
The first and only line of input contains the string s. \(2 \le |s| \le 10^5\).
Output
The first line of output contains the maximal \(m_s( i )\) over all i.
The second line of output contains all the i's for which \(m_s( i )\) reaches maximum.
Example
Input:
caccacaa
Output:
4
3
Explanation:
caccacaa
caccacaa
The bold characters indicate the ones that match when shift = 3.
Solution
题意
给定一个字符串 \(s\) (下标从 \(0\) 开始,只包含 'a', 'b', 'c'),让 \(s\) 与 \(s\) 匹配,下标从 \(1\) 移动到 \(|s|\),每次匹配时的相同的字符个数记为 \(m_s( i )\),求 \(m_s( i )\) 的最大值,以及最大值所匹配的所有位置。
比如 ababa
ababa
ababa
ababa
ababa
ababa
\(m_s( i )\) 分别为 \(0, 3, 0, 1\),最大值为 \(3\)。
思路
FFT
字符串匹配问题。
设模式串为 \(p\),目标串为 \(t\),\(f[k]\) 为模式串从目标串第 \(k\) 位开始匹配的结果。
对 \(a\),\(b\),\(c\) 分开求。
首先判断 \(a\) 的情况,将字符串转化为 01 串,比如 ababa 转为 10101。
那么
\]
\]
将模式串倒置,有
\]
令 \(j = |s| - 1 - i + k\),有
\]
用 \(FFT\) 求一下卷积即可。
对于 \(b\) 和 \(c\) 的求法相同。
Code
#include <bits/stdc++.h>
using namespace std;
const double PI = acos(-1);
const double eps = 1e-8;
typedef complex<double> Complex;
const int maxn = 2e6 + 10;
Complex p[maxn], t[maxn];
Complex a[maxn], b[maxn], c[maxn];
int ans[maxn];
string str;
int n;
int bit = 2, rev[maxn];
void get_rev(){
memset(rev, 0, sizeof(rev));
while(bit <= n + n) bit <<= 1;
for(int i = 0; i < bit; ++i) {
rev[i] = (rev[i >> 1] >> 1) | (bit >> 1) * (i & 1);
}
}
void FFT(Complex *arr, int op) {
for(int i = 0; i < bit; ++i) {
if(i < rev[i]) swap(arr[i], arr[rev[i]]);
}
for(int mid = 1; mid < bit; mid <<= 1) {
Complex wn = Complex(cos(PI / mid), op * sin(PI / mid));
for(int j = 0; j < bit; j += mid<<1) {
Complex w(1, 0);
for(int k = 0; k < mid; ++k, w = w * wn) {
Complex x = arr[j + k], y = w * arr[j + k + mid];
arr[j + k] = x + y, arr[j + k + mid] = x - y;
}
}
}
}
int main() {
ios::sync_with_stdio(false);
cin.tie(0);
cin >> str;
n = str.size();
for(int i = 0; i < n; ++i) {
p[n - i - 1] = str[i] == 'a' ? 1 : 0;
t[i] = p[n - i - 1];
}
get_rev();
FFT(p, 1); FFT(t, 1);
for(int i = 0; i < bit; ++i) {
a[i] = p[i] * t[i];
}
FFT(a, -1);
for(int i = 0; i < n; ++i) {
p[n - i - 1] = str[i] == 'b' ? 1 : 0;
t[i] = p[n - i - 1];
}
for(int i = n; i < bit; ++i) {
p[i] = 0;
t[i] = 0;
}
FFT(p, 1); FFT(t, 1);
for(int i = 0; i < bit; ++i) {
b[i] = p[i] * t[i];
}
FFT(b, -1);
for(int i = 0; i < n; ++i) {
p[n - i - 1] = str[i] == 'c' ? 1 : 0;
t[i] = p[n - i - 1];
}
for(int i = n; i < bit; ++i) {
p[i] = 0;
t[i] = 0;
}
FFT(p, 1); FFT(t, 1);
for(int i = 0; i < bit; ++i) {
c[i] = p[i] * t[i];
}
FFT(c, -1);
int maxa = 0;
for(int i = 1; i < n; ++i) {
ans[i] = (int)(a[n - 1 + i].real() / bit + 0.5) + (int)(b[n - 1 + i].real() / bit + 0.5) + (int)(c[n - 1 + i].real() / bit + 0.5);
maxa = max(maxa, ans[i]);
}
vector<int> pos;
for(int i = 1; i < n; ++i) {
if(ans[i] == maxa) {
pos.push_back(i);
}
}
cout << maxa << endl;
for(int i = 0; i < pos.size(); ++i) {
cout << pos[i] << " ";
}
cout << endl;
return 0;
}
SPOJ MAXMATCH - Maximum Self-Matching (FFT)的更多相关文章
- Maximum Bipartite Matching
算法旨在用尽可能简单的思路解决这个问题.理解算法也应该是一个越看越简单的过程,当你看到算法里的一串概念,或者一大坨代码,第一感觉是复杂,此时最好还是从样例入手.通过一个简单的样例,并编程实现,这个过程 ...
- SPOJ TSUM Triple Sums(FFT + 容斥)
题目 Source http://www.spoj.com/problems/TSUM/ Description You're given a sequence s of N distinct int ...
- SPOJ - VFMUL - Very Fast Multiplication FFT加速高精度乘法
SPOJ - VFMUL:https://vjudge.net/problem/SPOJ-VFMUL 这是一道FFT求高精度的模板题. 参考:https://www.cnblogs.com/Rabbi ...
- 2018.11.18 spoj Triple Sums(容斥原理+fft)
传送门 这次fftfftfft乱搞居然没有被卡常? 题目简述:给你nnn个数,每三个数ai,aj,ak(i<j<k)a_i,a_j,a_k(i<j<k)ai,aj,ak( ...
- 牛客网暑期ACM多校训练营(第三场)DEncrypted String Matching fft
题意:给你一个解密后的字符串,给你加密方式,加密过程可能出错,字符可能加减1,然后给你一个字符串,要求匹配个数(其实我也不太懂具体怎么加密解密,反正你把给你的前两个字符串用第三个加密一下,然后搞可以有 ...
- SPOJ VFMUL - Very Fast Multiplication (FFT)
题目链接:VFMUL - Very Fast Multiplication Description Multiply the given numbers. Input n [the number of ...
- SPOJ 3693 Maximum Sum(水题,记录区间第一大和第二大数)
#include <iostream> #include <stdio.h> #include <algorithm> #define lson rt<< ...
- 转债---Pregel: A System for Large-Scale Graph Processing(译)
转载:http://duanple.blog.163.com/blog/static/70971767201281610126277/ 作者:Grzegorz Malewicz, Matthew ...
- Pregel: A System for Large-Scale Graph Processing(译)
[说明:Pregel这篇是发表在2010年的SIGMOD上,Pregel这个名称是为了纪念欧拉,在他提出的格尼斯堡七桥问题中,那些桥所在的河就叫Pregel.最初是为了解决PageRank计算问题,由 ...
随机推荐
- java.net.BindException: Address already in use: 解决方法
java.net.BindException: Address already in use: 解决方法 1. 执行cmd 2. cmd命令模式下输入netstat -ano,然后找到占用端口的那 ...
- JSON.stringify,JSON.parse方法
var obj={name:'zhangsan',age:'18'};/** js对象--->JSON字符串* JSON.stringify(js对象) --转化为--> JSON字符串* ...
- java多线程学习笔记(四)
上一节讲到Synchronized关键字,synchronized上锁的区域:对象锁=方法锁/类锁 本节补充介绍一下synchronized锁重入: 关键字synchronized拥有锁重入的功能,也 ...
- docker内的服务无法获取用户真实IP
原文:blog.baohaipeng.top 背景:MySQL数据库和Redis运行在宿主机上(Linux),server运行在docker内,web运行在Nginx内(Nginx运行在docker内 ...
- 使用org-mode写cnblogs博客
使用org-mode写cnblogs博客 */--> pre.src {background-color: #002b36; color: #839496;} pre.src {backgrou ...
- im开发总结:netty的使用
最近公司在做一个im群聊的开发,技术使用得非常多,各种代码封装得也是十分优美,使用到了netty,zookeeper,redis,线程池·,mongdb,lua,等系列的技术 netty是对nio的一 ...
- call和apply实现的继承
<!DOCTYPE html> <html> <head lang="en"> <meta charset="UTF-8&quo ...
- TreeMap和TreeSet在排序时如何比较元素,Collections工具类中的sort()方法如何比较元素
TreeSet和TreeMap排序时比较元素要求元素对象必须实现Comparable接口 Collections的sort方法比较元素有两种方法: 元素对象实现Comparable接口 实体类Dog ...
- 案例-3D旋转木马
<!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...
- SpringMVC 与axis2 的整合(服务器端)
1,新建一个web project项目,项目的目录 如下: 2,导入需要的jar包,本例用的是axis2 1.7.3版本. 3,创建接口. package com.dsdl.hcm.webServic ...