【jzoj5093】【GDSOI2017第四轮模拟day3】【字符串匹配】【哈希】
2017-04-26 21:22
330 查看
题目大意
对于一个字符集大小为C的字符串P,我们可以将任意两种字符在P中的位置进行互换,例如P=abcba,我们交换a,b就变为bacab,交换a,d就变为dbcbd,交换可以进行任意次。若交换后P变为了字符串Q,则我们称Q与P是匹配的。
现在给定两个字符集大小为C的字符串S,T,请你求出S中有多少个连续子串与T是匹配的。
解题思路
每个位的值为与前一个相同颜色的位置差,对这个哈希即可判断是否同构,特殊处理一下关联关系即可支持加入和删除。
code
#include<cstdio> #include<algorithm> #define LL long long #define ULL unsigned long long #define fo(i,j,k) for(int i=j;i<=k;i++) using namespace std; int const mn=1e6+9,inf=1e9+7; int t,c,n,m,ans[mn],a[mn],b[mn],pre[mn],next[mn]; ULL p=1e9+7,bit[mn]; int main(){ //freopen("string.in","r",stdin); //freopen("string.out","w",stdout); freopen("d.in","r",stdin); freopen("d.out","w",stdout); scanf("%d%d",&t,&c); bit[0]=1;fo(i,1,1000000)bit[i]=bit[i-1]*p; fo(cas,1,t){ scanf("%d%d",&n,&m); fo(i,1,n)scanf("%d",&a[i]); fo(i,1,m)scanf("%d",&b[i]); fo(i,1,c)pre[i]=0; ULL hashb=0; fo(i,1,m){ int tmp=(pre[b[i]])?i-pre[b[i]]:0; hashb=hashb*p+tmp; pre[b[i]]=i; } fo(i,1,c)pre[i]=0; fo(i,1,n)next[i]=0; ULL hasha=0; fo(i,1,m){ int tmp=(pre[a[i]])?i-pre[a[i]]:0; hasha=hasha*p+tmp; next[pre[a[i]]]=i; pre[a[i]]=i; } ans[0]=0; if(hasha==hashb)ans[++ans[0]]=1; fo(i,m+1,n){ int ii=i-m,tmp=(pre[a[i]]>ii)?i-pre[a[i]]:0,tmp2=(next[ii])?next[ii]-ii:0; hasha=(hasha-tmp2*bit[i-1-next[ii]])*p+tmp; if(hasha==hashb)ans[++ans[0]]=ii+1; next[pre[a[i]]]=i; pre[a[i]]=i; } printf("%d\n",ans[0]); fo(i,1,ans[0])printf("%d ",ans[i]); printf("\n"); } return 0; }
相关文章推荐
- 【jzoj5094】【GDSOI2017第四轮模拟day3】【鸽子】【计算几何】
- NOIP2016A组模拟7.13】字符串匹配
- 一种字符串哈希方法 与o(n+m)字符串匹配
- [ 备战NOIP2016 ] 字符串匹配 - 哈希
- 【JZOJ5068】【GDSOI2017第二轮模拟】树
- 【JZOJ4624】字符串匹配
- 【jzoj5088】【GDOI2017第四轮模拟day2】【最小边权和 】【图论】
- 哈希 字符串匹配
- 【jzoj5086】【GDOI2017第四轮模拟day1】【数列】【搜索】
- UVa 489 Hangman Judge (模拟&字符串匹配)
- 【JZOJ 4624】字符串匹配
- 【NOIP2016A组模拟7.13】字符串匹配 题解+代码
- 【JZOJ5069】【GDSOI2017第二轮模拟】蛋糕
- 【JZOJ5056】【GDSOI2017模拟4.13】黑白广场
- 算法模板——哈希单模板字符串匹配
- 【jzoj5085】【GDOI2017第四轮模拟day1】【游戏】 【sg函数】【线性基】
- JZOJ 5068. 【GDSOI2017第二轮模拟】树
- [JZOJ5082].【GDSOI2017第三轮模拟】Informatics Training
- 【jzoj4624】【字符串匹配】【kmp】
- 简单字符串匹配问题-——用哈希解决