题目大意:

给定A,B两种字符串,问他们当中的长度大于k的公共子串的个数有多少个

这道题目本身理解不难,将两个字符串合并后求出它的后缀数组

然后利用后缀数组求解答案

这里一开始看题解说要用栈的思想,觉得很麻烦就不做了,后来在比赛中又遇到就后悔了,到今天看了很久才算看懂

首先建一个栈,从栈底到栈顶都保证是单调递增的

我们用一个tot记录当前栈中所有项和一个刚进入的子串匹配所能得到的总的子串的数目(当然前提是,当前进入的子串height值比栈顶还大,那么和栈中任意一个子串匹配都保持当前栈中记录的那时候入栈的height值)

但是若height不比栈顶大,说明从栈顶开始到刚好比它小的这一段tot有多加的部分,这部分就是height值多出来的那块,然后把这部分都视作height值为当前的height值,因为后面子串进入,它的height值总是取决于那段区间的最小值,所以不会产生影响,这样就可以把所有比当前height大的都弹出栈,这样就达到了O(n)的复杂度

这里用q[][]手写栈

q[i][0]表示栈中第i号元素记录时候的height值,q[i][1]表示在这个height值上覆盖了q[i][1]个子串

  1 #include <cstdio>
  2 #include <cstring>
  3 #include <iostream>
  4 using namespace std;
  5 #define INF 0x3f3f3f3f
  6 #define ll long long
  7 const int MAXN = 2*100010;
  8 int sa[MAXN] , rank[MAXN] , height[MAXN];
  9 int wa[MAXN] , wb[MAXN] , wsf[MAXN] , wv[MAXN];
 10 int a[MAXN] , k;
 11 char str1[MAXN] , str2[MAXN];
 12 int q[MAXN][2];
 13
 14 int cmp(int *r , int a , int b , int l)
 15 {
 16     return r[a]==r[b] && r[a+l]==r[b+l];
 17 }
 18
 19 void getSa(int *r , int *sa , int n , int m)
 20 {
 21     int *x = wa , *y = wb , *t;
 22     for(int i=0 ; i<m ; i++) wsf[i]=0;
 23     for(int i=0 ; i<n ; i++) wsf[x[i]=r[i]]++;
 24     for(int i=1 ; i<m ; i++) wsf[i]+=wsf[i-1];
 25     for(int i=n-1 ; i>=0 ; i--) sa[--wsf[x[i]]] = i;
 26
 27     int i,j,p=1;
 28     for(j=1 ; p<n ; j*=2 , m=p)
 29     {
 30         for(p=0 , i=n-j ; i<n ; i++) y[p++] = i;
 31         for(i=0 ; i<n ; i++) if(sa[i]>=j) y[p++] = sa[i]-j;
 32
 33         for(i=0 ; i<n ; i++) wv[i]=x[y[i]];
 34         for(i=0 ; i<m ; i++) wsf[i]=0;
 35         for(i=0 ; i<n ; i++) wsf[wv[i]]++;
 36         for(i=1 ; i<m ; i++) wsf[i]+=wsf[i-1];
 37         for(i=n-1 ; i>=0 ; i--) sa[--wsf[wv[i]]] = y[i];
 38
 39         for(t=x , x=y , y=t , x[sa[0]]=0 , p=1 , i=1; i<n ; i++)
 40             x[sa[i]] = cmp(y , sa[i-1] , sa[i] , j)?p-1:p++;
 41     }
 42     return ;
 43 }
 44
 45 void callHeight(int *r , int *sa , int n)
 46 {
 47     for(int i=0 ; i<=n ; i++) rank[sa[i]]=i;
 48     int i , j , k=0;
 49     for(i=0 ; i<n ; height[rank[i++]]=k)
 50         for(j=sa[rank[i]-1] , k?k--:0 ; r[i+k]==r[j+k] ; k++) ;
 51     return;
 52 }
 53
 54 ll solve(int len1 , int len2)
 55 {
 56     ll ans = 0;
 57     //B串中的子串不断匹配rank比其高的A子串
 58     int top = 0;
 59     ll tot =0 , cnt = 0;
 60     for(int i=1 ; i<=len1+len2+1 ; i++){
 61         if(height[i]<k){
 62             top = tot = 0;
 63             continue;
 64         }
 65         cnt = 0;
 66         if(sa[i-1]<len1){
 67             cnt ++;
 68             tot += height[i]-k+1;
 69         }
 70         while(top&&height[i]<=q[top][0]){
 71             tot -= q[top][1]*(q[top][0]-height[i]);
 72             cnt += q[top][1];
 73             top--;
 74         }
 75         q[++top][0] = height[i];
 76         q[top][1] = cnt;
 77         if(sa[i]>len1) ans+=tot;
 78     }
 79     //A串中的子串不断匹配rank比其高的B子串
 80     tot = top = 0;
 81     for(int i=1 ; i<=len1+len2+1 ; i++){
 82         if(height[i]<k){
 83             top = tot = 0;
 84             continue;
 85         }
 86         cnt = 0;
 87         if(sa[i-1]>len1){
 88             cnt ++;
 89             tot += height[i]-k+1;
 90         }
 91         while(top&&height[i]<=q[top][0]){
 92             tot -= q[top][1]*(q[top][0]-height[i]);
 93             cnt += q[top][1];
 94             top--;
 95         }
 96         q[++top][0] = height[i];
 97         q[top][1] = cnt;
 98         if(sa[i]<len1) ans+=tot;
 99     }
100     return ans;
101 }
102
103 int main()
104 {
105    // freopen("a.in" , "r" , stdin);
106
107     while(scanf("%d" , &k) , k)
108     {
109         scanf("%s%s" , str1 , str2);
110         int len1 = strlen(str1) , len2 = strlen(str2);
111         for(int i=0 ; i<len1 ; i++) a[i] = (int)str1[i];
112         a[len1] = 259;
113         for(int i=0 ; i<len2 ; i++) a[i+len1+1] = (int)str2[i];
114         a[len1+len2+1] = 0;
115
116         getSa(a , sa , len1+len2+2 , 260);
117         callHeight(a , sa , len1+len2+1);
118
119       //  for(int i=0 ; i<len1+len2+2 ; i++) cout<<"rank i: "<<i<<" "<<rank[i]<<endl;
120       //  for(int i=1 ; i<len1+len2+2 ; i++) cout<<"xixi: "<<height[i]<<endl;
121         ll ans = solve(len1 , len2);
122         printf("%I64d\n" , ans);
123     }
124     return 0;
125 }

转载于:https://www.cnblogs.com/CSU3901130321/p/4516109.html

POJ 3415 后缀数组+单调栈相关推荐

  1. POJ - 3415 Common Substrings(后缀数组+单调栈)

    题目链接:点击查看 题目大意:给出两个字符串,再给出一个k,问两个字符串中长度大于等于k的公共子串有多少个(种类可重复) 题目分析:因为涉及到了子串问题,先用后缀数组跑出height数组来,接下来如果 ...

  2. [Ahoi2013]差异[后缀数组+单调栈]

    链接 解题思路:很明显前面∑1<=i<j<=nlen(Ti)+len(Tj)\sum_{1<=i<j<=n}len(T_i)+len(T_j)∑1<=i< ...

  3. [BZOJ 3238] [AHOI 2013] 差异 【后缀数组 + 单调栈】

    题目链接:BZOJ - 3238 题目分析 显然,这道题就是求任意两个后缀之间的LCP的和,这与后缀数组的联系十分明显. 求出后缀数组后,求出字典序相邻两个后缀的LCP,即 Height 数组. 那么 ...

  4. BZOJ3879: SvT【后缀数组+单调栈】

    Description (我并不想告诉你题目名字是什么鬼) 有一个长度为n的仅包含小写字母的字符串S,下标范围为[1,n]. 现在有若干组询问,对于每一个询问,我们给出若干个后缀(以其在S中出现的起始 ...

  5. [bzoj3238]差异(后缀数组+单调栈)

    显然我们可以先把len(Ti)+len(Tj)的值先算出来,再把LCP减去.所有len(Ti)+len(Tj)的值为n*(n-1)*(n+1)/2,这个随便在纸上画一画就可以算出来的. 接下来问题就是 ...

  6. [BZOJ3238][AHOI2013]差异 [后缀数组+单调栈]

    题目地址 - GO-> 题目大意: 给定一个长度为 nn 的字符串SS,令TiTi表示它从第ii个字符开始的后缀,求以下这个式子的值: ∑1≤i<j≤nlen(Ti)+len(Tj)−2× ...

  7. 【BZOJ3879】SvT,后缀数组+单调栈维护sum

    Time:2016.08.15 Author:xiaoyimi 转载注明出处谢谢 如果有不明白的地方可以在下面评论问我 传送门 思路: 建立后缀数组求出Height 如果说对每次询问暴力求LCP累加的 ...

  8. bzoj 3238: [Ahoi2013]差异(后缀数组+单调栈)

    3238: [Ahoi2013]差异 Time Limit: 20 Sec  Memory Limit: 512 MB Submit: 3443  Solved: 1562 [Submit][Stat ...

  9. 【HAOI2016/BZOJ4566】找相同字符 后缀数组+单调栈

    原题走这里 鉴于我实在不是很懂单调栈和单调队列这一系列东西,所以我决定稍微具体讲一下单调栈. 恩,本题实质上就是求两个字符串的公共子串数,其中只要出现位置不同,就算是不同的子串. 处理多个字符串的经典 ...

最新文章

  1. Java FAQ(6)
  2. dplyr 数据操作 常用函数(2)
  3. MySQL 的存储引擎
  4. InstallShild的研究,msde2000,.netframwork,ScriptProject与ScriiptMSIProject的区别
  5. Java编程:Java的反射机制中的 getComponentType() 方法
  6. linux 从github拉取更新_关于拉取请求
  7. 飞冰:Iceworks 自定义模板支持布局定制(v2.3.0 版本)
  8. 这些藏在成都的 NB 互联网公司
  9. 搜狗推送工具-搜狗批量提交软件
  10. Python 爬虫超详细讲解(零基础入门,老年人都看的懂)
  11. oracle 同义词public,oracle中private同义词跟public同义词
  12. Linux下ppp拨号+电信3G模块
  13. Python本地文件合并(csv)
  14. 廊坊金彩教育:如何提高店铺评分
  15. 百度网盘微信小程序文件 同步到百度网盘APP或客户端
  16. 百度地图API基础操作--导航篇
  17. 易语言 vb c 那个写的程序运行快,VB好还是易语言
  18. html代码打猎小游戏,html贪吃蛇小游戏课件-附源代码
  19. Apache的性能解读
  20. 【java并发探赜索隐】判断标志、死锁、守护进程

热门文章

  1. linux windows死机,linux死机的处理方式
  2. 从 0 到 1 实现 React 系列 —— 组件和 state|props
  3. Exchange Server 2016管理系列课件17.创建和使用会议室邮箱
  4. 【转载】C#中回滚TransactionScope的使用方法和原理
  5. python程序编译成exe格式
  6. 如何夺回Windows 10系统权限?一分钟就够
  7. Spring-IOC
  8. 10倍提升应用性能的10个建议
  9. 二叉查找树的类模板实现
  10. 堆结构导致数据文件不能收缩