Description

题库链接

给定两个只含小写字母字符串 \(s_1,s_2\) ,求出在两个字符串中各取出一个子串使得这两个子串相同的方案数。两个方案不同当且仅当这两个子串中有一个位置不同。

\(1\leq |s_1|,|s_2|\leq 200000\)

Solution

把两个串拼在一起,然后扫两遍。

第一遍统计所有的 \(B\) ,在它前面的 \(A\) 的贡献;第二遍统计所有的 \(A\) ,在它前面的 \(B\) 的贡献;

用单调栈维护一下即可。

Code

  1. #include <bits/stdc++.h>
  2. #define ll long long
  3. using namespace std;
  4. const int N = (200000+5)<<1;
  5. char ch[N];
  6. int s1, n, m, x[N<<1], y[N<<1], c[N], sa[N], rk[N], height[N];
  7. ll ans, sum[N], a[N];
  8. int s[N], top;
  9. void get() {
  10. for (int i = 1; i <= n; i++) c[x[i] = ch[i]]++;
  11. for (int i = 2; i <= m; i++) c[i] += c[i-1];
  12. for (int i = n; i >= 1; i--) sa[c[x[i]]--] = i;
  13. for (int k = 1; k <= n; k <<= 1) {
  14. int num = 0;
  15. for (int i = n-k+1; i <= n; i++) y[++num] = i;
  16. for (int i = 1; i <= n; i++) if (sa[i] > k) y[++num] = sa[i]-k;
  17. for (int i = 1; i <= m; i++) c[i] = 0;
  18. for (int i = 1; i <= n; i++) c[x[i]]++;
  19. for (int i = 2; i <= m; i++) c[i] += c[i-1];
  20. for (int i = n; i >= 1; i--) sa[c[x[y[i]]]--] = y[i];
  21. swap(x, y); x[sa[1]] = num = 1;
  22. for (int i = 2; i <= n; i++)
  23. x[sa[i]] = (y[sa[i]] == y[sa[i-1]] && y[sa[i]+k] == y[sa[i-1]+k]) ? num : ++num;
  24. if ((m = num) == n) break;
  25. }
  26. for (int i = 1; i <= n; i++) rk[sa[i]] = i;
  27. for (int i = 1, k = 0; i <= n; i++) {
  28. if (rk[i] == 0) continue;
  29. if (k) k--; int j = sa[rk[i]-1];
  30. while (i+k <= n && j+k <= n && ch[i+k] == ch[j+k]) ++k;
  31. height[rk[i]] = k;
  32. }
  33. }
  34. void cal(int ss) {
  35. top = 0;
  36. for (int i = 1; i <= n; i++) sum[i] = sum[i-1]+(ss == (sa[i] > s1) && sa[i] != s1+1);
  37. for (int i = 1; i <= n; i++) {
  38. while (top && height[s[top]] >= height[i]) --top;
  39. s[++top] = i; a[top] = a[top-1]+(sum[i-1]-sum[s[top-1]-1])*height[i];
  40. if (ss == (sa[i] <= s1) && sa[i] != s1+1) ans += a[top];
  41. }
  42. }
  43. void work() {
  44. scanf("%s", ch+1); s1 = strlen(ch+1);
  45. ch[s1+1] = '$'; scanf("%s", ch+s1+2);
  46. n = strlen(ch+1); m = 255; get();
  47. cal(0), cal(1); printf("%lld\n", ans);
  48. }
  49. int main() {work(); return 0; }

[HAOI 2016]找相同字符的更多相关文章

  1. [bzoj 4566][Haoi 2016]找相同字符

    传送门 Description 给定两个字符串,求出在两个字符串中各取出一个子串使得这两个子串相同的方案数.两个方案不同当且仅当这两个子串中有一个位置不同. Solution 考虑用SAM,采用较为轻 ...

  2. BZOJ4566&&lg3181 HAOI找相同字符(广义后缀自动机)

    BZOJ4566&&lg3181 HAOI找相同字符(广义后缀自动机) 题面 自己找去 HINT 给定两个文本串,问从两个串中各取一个非空子串,使这俩子串相同,问方案有多少种.我的思路 ...

  3. BZOJ 4566: [Haoi2016]找相同字符 [后缀自动机]

    4566: [Haoi2016]找相同字符 Time Limit: 20 Sec  Memory Limit: 256 MBSubmit: 275  Solved: 155[Submit][Statu ...

  4. 【BZOJ4566】找相同字符(后缀数组)

    [BZOJ4566]找相同字符(后缀数组) 题面 BZOJ 题解 后缀数组的做法,应该不是很难想 首先看到两个不同的串,当然是接在一起求\(SA,height\) 那么,考虑一下暴力 在两个串各枚举一 ...

  5. 【BZOJ4566】找相同字符(后缀自动机)

    [BZOJ4566]找相同字符(后缀自动机) 题面 BZOJ 题解 看到多串处理,\(SA\)就连起来 \(SAM???\) 单串建自动机 然后其他串匹配 对于一个串建完\(SAM\)后 另一个串在\ ...

  6. bzoj4566 / P3181 [HAOI2016]找相同字符

    P3181 [HAOI2016]找相同字符 后缀自动机 (正解应是广义后缀自动机) 并不会广义后缀自动机. 然鹅可以用普通的后缀自动机.   我们先引入一个问题:算出从一个串内取任意两个不重合子串完全 ...

  7. 【BZOJ4566】[HAOI2016]找相同字符

    [BZOJ4566][HAOI2016]找相同字符 题面 给定两个字符串,求出在两个字符串中各取出一个子串使得这两个子串相同的方案数.两个方案不同当且仅当这两个子串中有一个位置不同. 其中\(1\le ...

  8. [BZOJ4566][Haoi2016]找相同字符 后缀自动机+dp

    4566: [Haoi2016]找相同字符 Time Limit: 20 Sec  Memory Limit: 256 MBSubmit: 1212  Solved: 694[Submit][Stat ...

  9. 【BZOJ4566】[Haoi2016]找相同字符 后缀数组+单调栈

    [BZOJ4566][Haoi2016]找相同字符 Description 给定两个字符串,求出在两个字符串中各取出一个子串使得这两个子串相同的方案数.两个方案不同当且仅当这两 个子串中有一个位置不同 ...

随机推荐

  1. PyCharm选择性忽略PEP8代码风格警告信息

    用了几天的PyCharm,发现确实在编写Python代码上非常好用,但有一点体验不太好,就是代码编写时要按照PEP8代码风格编写,不然会有波浪线的警告信息.解决方法如下: 方法一: 将鼠标移到提示的地 ...

  2. MySQL缓存参数优化(转)

    MySQL 数据库性能优化之缓存参数优化 数据库属于 IO 密集型的应用程序,其主要职责就是数据的管理及存储工作.而我们知道,从内存中读取一个数据库的时间是微秒级别,而从一块普通硬盘上读取一个IO是在 ...

  3. CUDA[3] Samples for accessing shared/global memory

    memory model: programming model: Source:  Udacity Class CS344

  4. MySQL常用内置变量

      MySQL用很多常用的内置变量,掌握这些内置变量后对于我们快速获取当前MySQL的配置有很大帮助,下面就来列举几个常用的变量. 查看当前MySQL版本号信息.show variables like ...

  5. Jquery - 添加属性、添加class、添加Css

    一.设置属性: 方式一  jQuery 代码: $("img").attr({ src: "test.jpg", alt: "Test Image&q ...

  6. 770. Basic Calculator IV

    Given an expression such as expression = "e + 8 - a + 5" and an evaluation map such as {&q ...

  7. 696. Count Binary Substrings

    Give a string s, count the number of non-empty (contiguous) substrings that have the same number of ...

  8. 二进制加法Java实现

    二进制整数的Java实现 任意两个二进制数(不论在什么位置)相加,只可能出现4种情况.它们是: 0+0=0 1+0=0+1=1 1+1=10=0+向高一位的进位1 1+1+1=11=1+向高一位的进位 ...

  9. Scala微服务架构 二

    三. Scala的Macro(宏) Scala Macros对scala函数库编程人员来说是一项不可或缺的编程工具,可以通过它来解决一些用普通编程或者类层次编程(type level programm ...

  10. Java 虚拟机的内存结构

    Java虚拟机运行时数据区 整个程序执行过程中,JVM会用一段空间来存储程序执行期间需要用到的数据和相关信息,这段空间一般被称作为Runtime Data Area(运行时数据区),也就是我们常说的J ...