题意:问a串中有多少种字符串集合B中没有的连续子串。

a的长度10^5,B中的总长度为不超过10^5.

解法:后缀数组题目;后缀数组能够非常easy算出来一个串中有多少种子串。

把a和B集合连起来。求一次不同子串数量,然后减掉B相互连起来的数量。

在求时候,要减掉含有链接符的子串,方法是扫一遍,枚举最后出现的连接符。

代码:

/******************************************************
* @author:xiefubao
*******************************************************/
#pragma comment(linker, "/STACK:102400000,102400000")
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <queue>
#include <vector>
#include <algorithm>
#include <cmath>
#include <map>
#include <set>
#include <stack>
#include <string.h>
//freopen ("in.txt" , "r" , stdin);
using namespace std;#define eps 1e-8
#define zero(_) (abs(_)<=eps)
const double pi=acos(-1.0);
typedef long long LL;
const LL INF=0x3FFFFFFF;const int MAX=300010;
int n, num[MAX];
char s[MAX];
int sa[MAX], rank[MAX], height[MAX];//sa[i]表示排名第i的后缀的位置,height[i]表示后缀SA[i]和SA[i-1]的最长公共前缀
int wa[MAX], wb[MAX], wv[MAX], wd[MAX];
/* *suffix array *倍增算法  O(n*logn) *待排序数组长度为n,放在0~n-1中,在最后面补一个0 *da(str ,n+1,sa,rank,height,  ,   );//注意是n+1; *比如: *n   = 8; *num[]   = { 1, 1, 2, 1, 1, 1, 1, 2, $ };注意num最后一位为0,其它大于0 *rank[]  = { 4, 6, 8, 1, 2, 3, 5, 7, 0 };
rank[0~n-1]为有效值,rank[n]必然为0无效 值 *sa[]    = { 8, 3, 4, 5, 0, 6, 1, 7, 2 };sa[1~n]为有效值。sa[0]必然为n是无效值 *height[]= { 0, 0, 3, 2, 3, 1, 2, 0, 1 };height[2~n]为有效值 * */int t1[MAX],t2[MAX],c[MAX];//求SA数组须要的中间变量,不须要赋值 //待排序的字符串放在s数组中,从s[0]到s[n-1],长度为n,且最大值小于m,
//除s[n-1]外的全部s[i]都大于0。r[n-1]=0 //函数结束以后结果放在sa数组中
bool cmp(int *r,int a,int b,int l)
{return r[a] == r[b] && r[a+l] == r[b+l];
}
void da(int str[],int n,int m)
{n++;int i, j, p, *x = t1, *y = t2;     //第一轮基数排序,假设s的最大值非常大,可改为高速排序for(i = 0; i < m; i++)c[i] = 0;for(i = 0; i < n; i++)c[x[i] = str[i]]++;for(i = 1; i < m; i++)c[i] += c[i-1];for(i = n-1; i >= 0; i--)sa[--c[x[i]]] = i;for(j = 1; j <= n; j <<= 1){p = 0;         //直接利用sa数组排序第二keywordfor(i = n-j; i < n; i++)y[p++] = i;//后面的j个数第二keyword为空的最小for(i = 0; i < n; i++) if(sa[i] >= j) y[p++] = sa[i] - j;         //这样数组y保存的就是依照第二keyword排序的结果//基数排序第一keywordfor(i = 0; i < m; i++)c[i] = 0;for(i = 0; i < n; i++)c[x[y[i]]]++;for(i = 1; i < m; i++)c[i] += c[i-1];for(i = n-1; i >= 0; i--)sa[--c[x[y[i]]]] = y[i];        //依据sa和x数组计算新的x数组swap(x,y);p = 1;x[sa[0]] = 0;for(i = 1; i < n; i++)x[sa[i]] = cmp(y,sa[i-1],sa[i],j)?p-1:p++;if(p >= n)break;m = p;//下次基数排序的最大值}int k = 0;n--;for(i = 0; i <= n; i++)rank[sa[i]] = i;for(i = 0; i < n; i++){if(k)k--;j = sa[rank[i]-1];while(str[i+k] == str[j+k])k++;height[rank[i]] = k;}
}void da1(int *r, int n, int m)            //  倍增算法0(nlgn)。
{int i, j, p, *x = wa, *y = wb, *t;for(i = 0; i < m; i ++) wd[i] = 0;for(i = 0; i < n; i ++) wd[x[i]=r[i]] ++;for(i = 1; i < m; i ++) wd[i] += wd[i-1];for(i = n-1; i >= 0; i --) sa[-- wd[x[i]]] = i;for(j = 1, p = 1; p < n; j *= 2, m = p){for(p = 0, i = n-j; i < n; i ++) y[p ++] = i;for(i = 0; i < n; i ++) if(sa[i] >= j) y[p ++] = sa[i] - j;for(i = 0; i < n; i ++) wv[i] = x[y[i]];for(i = 0; i < m; i ++) wd[i] = 0;for(i = 0; i < n; i ++) wd[wv[i]] ++;for(i = 1; i < m; i ++) wd[i] += wd[i-1];for(i = n-1; i >= 0; i --) sa[-- wd[wv[i]]] = y[i];for(t = x, x = y, y = t, p = 1, x[sa[0]] = 0, i = 1; i < n; i ++){x[sa[i]] = cmp(y, sa[i-1], sa[i], j) ? p - 1: p ++;}}
}LL get(int* p,int m,int hh)
{da(p,m,hh+2);LL ans=0;for(int i=2; i<=m; i++){ans+=m-sa[i-1]-height[i];}ans+=m-sa[m];LL last=-1;for(int i=0; i<m; i++){if(p[i]>=39){if(last==-1)last=i;else{ans-=(i-last)*(last+1);last=i;}}}if(last!=-1)ans-=(m-last)*(last+1);return ans;
}
int main()
{int t;cin>>t;int kk=1;while(t--){scanf("%d",&n);scanf("%s",s);int len=strlen(s);int sum=0;for(int i=0; i<len; i++)num[sum++]=s[i]-'a'+1;int hh=39;num[sum++]=hh++;for(int i=0; i<n; i++){scanf("%s",s);int le=strlen(s);for(int j=0; j<le; j++){num[sum++]=s[j]-'a'+1;}if(i!=n-1)num[sum++]=hh++;}num[sum]=0;LL ans1=get(num,sum,hh);LL ans2=get(num+len+1,sum-len-1,hh);printf("Case %d: ",kk++);cout<<ans1-ans2<<endl;}return 0;
}

转载于:https://www.cnblogs.com/zsychanpin/p/6711016.html

hdu4416 Good Article Good sentence (后缀数组)相关推荐

  1. [hdu4416 Good Article Good sentence]后缀自动机SAM

    题意:给出串A和串集合B={B1,B2,...,Bn},求串A的所有不同子串中不是B中任一串的子串的数目. 思路:把A和B中所有字符串依次拼接在一起,然后构造后缀自动机,计算每个状态的R集合元素的最大 ...

  2. hdu4416 Good Article Good sentence(多个串的本质不同的串个数,后缀自动机)

    题意: 给定串S,和若干个串T(i), 问有多少个S的本质不同的子串,在所有T(i)中没有出现过 数据范围:|S|<=1e5,,sum(|T|)<=1e5 解法: 将所有T串插入后缀自动机 ...

  3. HDU 4416 Good Article Good sentence(12年杭州 后缀数组)

    转载请注明出处,谢谢http://blog.csdn.net/acm_cxlove/article/details/7854526       by---cxlove 题目:给出一个A串,给出若干个B ...

  4. 寻找一个字符串的重复子串 后缀数组

    什么是后缀数组 令字符串 S=S[1]S[2]...S[n]S=S[1]S[2]...S[n]{\displaystyle S=S[1]S[2]...S[n]} , S[i,j]S[i,j]{\dis ...

  5. hiho一下120周 后缀数组一·重复旋律

    后缀数组一·重复旋律 时间限制:5000ms 单点时限:1000ms 内存限制:256MB 描述 小Hi平时的一大兴趣爱好就是演奏钢琴.我们知道一个音乐旋律被表示为长度为 N 的数构成的数列. 小Hi ...

  6. 字符串相关处理kmp,前缀数,后缀树,后缀数组,最长回文串,最长重复字串,最长非重复字串

    1. 最长回文串 一般用后缀数组或者后缀树可以解决, 用此方法:http://blog.csdn.net/v_july_v/article/details/6897097 预处理后缀树,使得查询LCA ...

  7. 牛客多校 - B-Suffix Array(后缀数组)

    题目链接:点击查看 题目大意:给出一个只含有 ' a ' 和 ' b ' 的字符串,再给出 B 数组的构造方法如下:对于每个位置 i 来说 如果存在一个位置 j ,使得 j < i 且 s[ j ...

  8. cf244D. Match amp; Catch 字符串hash (模板)或 后缀数组。。。

    D. Match & Catch 能够用各种方法做.字符串hash.后缀数组,dp.拓展kmp,字典树.. . 字符串hash(模板) http://blog.csdn.net/gdujian ...

  9. 后缀数组(bzoj 1031: [JSOI2007]字符加密Cipher)

    后缀数组主要功能: 长度为n的字符串总共有n个后缀,求这n个后缀的字典序 实现方法:倍增+基数排序,过程就是下面那张表 求log(n)次rank数组,每次的rank数组都可以通过上次的rank数组得出 ...

最新文章

  1. Java之构造器和构造方法的使用和意义
  2. beautiful loss function
  3. zoj 3547 The Boss on Mars
  4. mysql 8.0.11-Windows (x86, 64-bit)下载地址与安装教程
  5. 用aliplayer如何实现视频的连续播放?
  6. mysql导出可执行文件_实例讲解MYSQL注入中导出可执行文件至启动项原理(三)...
  7. 基本数据类型与格式化输出
  8. Linux环境下NodeJS和MongoDB的安装配置
  9. Eclipse中添加Android系统jar包
  10. openwrt刷回原厂固件_小米路由器4刷breed, pandavan,openwrt
  11. Linux---带你区分根目录 和 家目录
  12. 解决安装Chrome翻译插件总是报CRX_HEADER_INVALID
  13. 抖音API分析,视频批量下载-初级
  14. Ruoyi框架学习--Vue前端配置文件详解
  15. Telnet 服务 开启 步骤
  16. DBD::mysql的安装
  17. pytorch tensor求向量的模长
  18. Tian Ji -- The Horse Racing 贪心算法
  19. HIve数仓新零售项目DWD层的构建
  20. 智能控制导论 # 模糊控制 - 理论基础:模糊数学 2 模糊关系和模糊推理

热门文章

  1. java输出孪生素数对数_为什么孪生素数、波利尼亚克、哥德巴赫这些素数加减关系的猜想不能被证明?问题出在哪?...
  2. python应该安装哪一个_python安装到哪个盘
  3. js bool true false 真假比较
  4. mini2440 之——LED简单控制(裸机) c与汇编
  5. windows下安装nodejs及框架express
  6. 约瑟夫环之循环链表实现
  7. HNCU1101:马的移动---BFS
  8. 防止头文件的重复包含
  9. 嵌入式系统启动流程和启动代码的作用 !
  10. C++的三大特性:封装,继承,多态