UVA4671 K-neighbor substrings FFT+字符串hash

题解:

将字符串A、B中的a和b分别以1和-1表示,对字符串B进行反转。
将A和B看成多项式,求卷积,这样的话从结果区间的[lenB1,lenA)[lenB−1,lenA)中的每一个点的值valval(lenBval)/2(lenB−val)/2代表当前位置的字串与B串的距离,然后对字串进行字符串hash,去重就是答案。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <cstring>
#include <set>
using namespace std;
double pi = acos(-1.0);
struct complex{double re,im;complex(double r = 0.0,double i = 0.0):re(r),im(i){};complex operator+(complex com){return complex(re+com.re,im+com.im);}complex operator-(complex com){return complex(re-com.re,im-com.im);}complex operator*(complex com){return complex(re*com.re-im*com.im,re*com.im+im*com.re);}
};
complex wn,wntmp;
void rader(complex arr[],int n){int num = n-1;for(int i = 0;i < n;++i){int tn = n>>1;while(num && num >= tn) num ^= tn,tn >>= 1;num |= tn;if(num > i) swap(arr[i],arr[num]);}
}
void FFT(complex cs[],int n,int f){rader(cs,n);for(int s = 1;s < n;s <<= 1){wn = complex(cos(f*2*pi/(s*2)),sin(f*2*pi/(s*2)));for(int offset = 0;offset < n;offset += s<<1){wntmp = complex(1.0,0.0);for(int i = 0;i < s;++i){complex u = cs[offset+i],v = cs[offset+i+s]*wntmp;cs[offset+i] = u + v;cs[offset+i+s] = u - v;wntmp = wntmp * wn;}}}if(f == -1)for(int i = 0;i < n;++i)cs[i].re /= n;
}
int K;
const int maxn = 600007;
char A[maxn],B[maxn];
complex csA[maxn],csB[maxn];
unsigned long long fac = 9973;
unsigned long long pow(int x){unsigned long long ans = 1,base = fac;while(x){if(x & 1)ans *= base;base *= base;x >>= 1;}return ans;
}
int main(){int cas = 0;while(cin>>K && K != -1){memset(csA,0,sizeof(csA)),memset(csB,0,sizeof(csB));cin>>A>>B;int lenA = strlen(A),lenB = strlen(B);for(int i = 0;i < lenB/2;++i) swap(B[i],B[lenB-i-1]);int len = 1;while(len < lenA || len < lenB) len <<= 1;len <<= 1;for(int i = 0;i < lenA;++i) csA[i].re = A[i] == 'a'?1:-1;FFT(csA,len,1);for(int i = 0;i < lenB;++i) csB[i].re = B[i] == 'a'?1:-1;FFT(csB,len,1);for(int i = 0;i < len;++i) csA[i] = csA[i]*csB[i];FFT(csA,len,-1);set<unsigned long long> st;unsigned long long hash = 0,base = pow(lenB-1);for(int i = 0;i < lenB;++i) hash = hash*fac + (A[i] == 'a');long long ans = 0; for(int i = lenB-1;i < lenA;++i){int dis= (lenB - int(csA[i].re+100000.5) + 100000)/2;if(dis <= K) st.insert(hash),ans++;hash = (hash - base * (A[i-lenB+1] == 'a'))*fac+(A[i+1] == 'a');}printf("Case %d: %d\n",++cas,st.size());}return 0;
}

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

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

相关文章

微软Build 2018展示Visual Studio功能:跨系统云编程

5 月 8 日凌晨消息&#xff0c;微软 Build 2018 开发者大会在今天正式来开帷幕。很明显整场发布会被分为了上下两个部分&#xff0c;上半场如果说是聊一些技术实现和愿景&#xff0c;那下半场就应该说是实战说明了。跨平台云端编程了解下微软在现场展示了 Visual Studio&#x…

【DP】花店橱窗布置 (ssl 1626/luogu 1854)

花店橱窗布置花店橱窗布置花店橱窗布置 ssl 1626 luogu 1854 题目大意&#xff1a; 有n朵花和m个花瓶&#xff0c;每一朵花插在不同的花瓶里会有不同的美丽度&#xff0c;花瓶和花必须按原来的顺序摆放&#xff0c;求最大的美丽值 Description 假设你想以最美观的方式布置…

P4819-[中山市选]杀人游戏【tarjan】

正题 题目链接:https://www.luogu.com.cn/problem/P4819 题目大意 nnn个人&#xff0c;一个杀手&#xff0c;搜查一个平民可以知道他认识的人的身份&#xff0c;搜查杀手就会死&#xff0c;求最优情况下警察的最低死亡概率。 解题思路 先用tarjantarjantarjan搜出强连通&…

hdu4609 3idiots 三角形计数 FFT

直接上代码 #include <iostream> #include <cstdio> #include <cmath> #include <cstring> #include <algorithm> using namespace std; double pi acos(-1.0); struct complex{double re,im;complex(double r 0.0,double i 0.0):re(r),im(i…

利用Skywalking-netcore监控你的应用性能

SkywalkingSkyWalking开源项目由吴晟于2015年创建&#xff0c;同年10月在GitHub上作为个人项目开源。SkyWalking项目的核心目标&#xff0c;是针对微服务、Cloud Native、容器化架构&#xff0c;提供应用性能监控&#xff08;APM&#xff09;和分布式调用链追踪能力。2017年11月…

【背包】逃亡的准备 (ssl 1236)

逃亡的准备逃亡的准备逃亡的准备 ssl 1236 类似一样题目&#xff08;除数组范围外&#xff09;:ssl 2289&#xff08;庆功会&#xff09; 题目大意&#xff1a; 有n个物品&#xff0c;每个物品可以选l[i]个,每个的代价为a[i],价值为b[i]&#xff0c;在代价不大于m的情况下&a…

P3648-[APIO2014]序列分割【斜率优化】

正题 题目链接:https://www.luogu.com.cn/problem/P3648 题目大意 nnn个数字的序列&#xff0c;分割kkk次&#xff0c;每次的权值是左右两块数字的乘积。求最大权值和分割方案。 解题思路 显然分割顺序不会影响结果&#xff0c;一个分割方式的答案是每一块与其他块的乘积之和…

codeforces E. Jamie and Tree LCA+dfs序+线段树

题解&#xff1a; 写起来还稍微有点麻烦。 dfs序线段树可以维护子树的整体修改和查询。 因此&#xff0c;这道题我们要往子树上靠。 我们首先从1号点进行dfs遍历&#xff0c;顺便求出点的dfs序和深度&#xff0c;然后我们采用倍增的思想&#xff0c;可以预处理出每个点的祖先…

async/await 的基本实现和 .NET Core 2.1 中相关性能提升

前言这篇文章的开头&#xff0c;笔者想多说两句&#xff0c;不过也是为了以后再也不多嘴这样的话。在日常工作中&#xff0c;笔者接触得最多的开发工作仍然是在 .NET Core 平台上&#xff0c;当然因为团队领导的开放性和团队风格的多样性&#xff08;这和 CTO 以及主管的个人能…

【DP】砝码称重 (ssl 1072)

砝码称重 ssl 1072 题目大意&#xff1a; 有6种砝码&#xff0c;每种的个数分别为a1,a2,a3,a4,a5,a6,请问可以平出多少种重量(0除外) Description 设有1g、2g、3g、5g、10g、20g的砝码各若干枚&#xff08;其总重<1000&#xff09;&#xff0c; 要求&#xff1a; 输入…

P5025-[SNOI2017]炸弹【tarjan,线段树优化建图】

正题 题目链接:https://www.luogu.com.cn/problem/P5025 题目大意 .nnn个炸弹&#xff0c;每个在xxx位置处&#xff0c;范围为rrr。定义fif_ifi​表示第iii个炸弹爆炸能连锁的炸弹数&#xff0c;要求输出∑i1nfi∗i\sum_{i1}^nf_i*ii1∑n​fi​∗i 解题思路 将每个炸弹可以炸…

使用Swashbuckle构建RESTful风格文档

本次和大家分享的是Swagger to WebApi的nuget包Swashbuckle&#xff1b;因为项目需要统一api文档的风格&#xff0c;并要支持多种开发语言&#xff08;C#&#xff0c;java&#xff0c;python&#xff09;&#xff0c;所以首先想到的是swagger来构建api文档&#xff0c;本章讲解…

【dfs】【bfs】【链表】 求连通分量 (ssl 1759)

求连通分量 ssl 1759 题目大意 由n个点组成的无向图&#xff0c;求连通在一起的点数最大是多少 原题 求一个图的连通分量 Input n 顶点数(<100) 边 Output 连通分量 Sample Input 8 6 3 1 2 2 5 5 4 4 1 8 7 0 0 Sample Output 4 方法一&#xff08;dfs …

P2472-[SCOI2007]蜥蜴【网络流】

正题 题目链接:https://www.luogu.com.cn/problem/P2472 题目大意 n∗mn*mn∗m个格子&#xff0c;每个格子的石柱高度不同&#xff0c;蜥蜴可以跳到距离不超过ddd的石柱处&#xff0c;并且先前所站的石柱高度减一&#xff0c;为0则不能站&#xff0c;然后求有多少只蜥蜴不可以…

发布 Rafy .NET Standard 版本 Nuget 包

去年年中&#xff0c;Rafy 框架的源码就已经支持了 Net Standard 2.0 版本。其开源代码也已经上传到 Github 中&#xff1a;https://github.com/zgynhqf/rafy/tree/NetStandard2.0 。但是这都只是在源码层面支持 NS2.0&#xff0c;并没有发布其正式的 Nuget 包。要使用这个版本…

一些来自STL的好东西

STL: 队列: 表达式作用#include&#xff1c;queue&#xff1e;定义queue&#xff1c;int&#xff1e;x定义一个int类型的队列&#xff0c;名为xx.push(y)从队列x的对尾插入yx.pop()队列x的队头出列gx.front()g等于x的队头hx.back()h等于x的队尾x.size()队列x的长度x.empty()判…

codeforces D.MADMAX 动态规划、记忆化搜索

题意 给出一个DAG&#xff0c;每条边上有权重(权重是小写字母的ASCII码)&#xff0c;现在两位同学A和B分别位于某两点上(可以相同)&#xff0c;其中A和B轮流走&#xff0c;但是每人所走的边权不能变小&#xff0c;走到不能走为止就输。 A先走&#xff0c;询问最后谁会赢。 题解…

GDOI2020游记

DAY0DAY0DAY0 早上上课一堆作业没写完(啊这 然后中午去竞赛室之前还没找到老师请假&#xff0c;还专门又跑回去拿了准考证(还白嫖了马老师红包 然后下午就去了&#xff0c;发现没带充电线 之后发现准考证没带&#xff0c;吓死了。还好酒店有可以打印的地方&#xff0c;还好搞定…

你关心才值得分享 | K8S网络安全之访问控制技术实践

(请允许我插播下广告&#xff0c;便于其它伙伴了解趣码 Cloud Coder)还是那句话&#xff0c;你关心才值得分享~最近的一起分享就在5.10本周四晚&#xff0c;精彩千万不要错过&#xff01;Hi&#xff0c;你是不是也曾觉得K8S&#xff08; Kubernetes &#xff09;网络安全话题范…

【dfs】【链表】连通图 (ssl 1758)

连通图 ssl 1758 题目大意 有一个由n个点组成的无向图&#xff0c;检测他是否联通 原题 判断一个图是否为一个边通图 Input n 顶点 (n<100) 边 Output 1 表示连通 0 表示不边通 Sample Input 5 1 2 2 3 5 4 0 0 Sample Output 0 解题方法 用dfs链表从1开始…