题意:求出字符串中包含了某个字符的字符序列不一样的数量。

思路:其实主要的是找出每个被包含字符的数量,假设除了目标字符之外的所有字符都不一样,那么应该就很好求了,但是显然不可能,所以我们可以枚举每一个起点,个数应该是从他的下一个字符是目标字符起的所有数量,但是通过观察我们可以发现这样计算我们又会多计算了一部分,例如a , abbabbabb 在计算第四个和第七个时,我们会多计算了a, ab, abb 或者计算第二位和第五位时多计算了bba,bbab,bbabb,我们可以这是就是相当于后缀数组里面的height数组。

#include<bits/stdc++.h>
using namespace std; typedef long long ll;
const int MAXN = ;
bool cmp(int *r, int a, int b, int l) {
return r[a] == r[b] && r[a + l] == r[b + l];
}
int t1[MAXN], t2[MAXN], c[MAXN]; void da(int str[], int sa[], int rk[], int height[], int n, int m) {
n++;
int i, j, p, *x = t1, *y = t2;
//第一轮基数排序,如果 s 的最大值很大,可改为快速排序
for(i = ; i < m; i++)
c[i] = ;
for(i = ; i < n; i++)
c[x[i] = str[i]]++;
for(i = ; i < m; i++)
c[i] += c[i - ];
for(i = n - ; i >= ; i -- )
sa[ -- c[x[i]]] = i;
for(j = ; j <= n; j <<= ) {
p = ;
//直接利用 sa 数组排序第二关键字
for(i = n - j; i < n; i++)
y[p++] = i;//后面的 j 个数第二关键字为空的最小
for(i = ; i < n; i++)
if(sa[i] >= j)
y[p++] = sa[i] - j;
//这样数组 y 保存的就是按照第二关键字排序的结果
//基数排序第一关键字
for(i = ; i < m; i++)
c[i] = ;
for(i = ; i < n; i++)
c[x[y[i]]]++;
for(i = ; i < m; i++)
c[i] += c[i - ];
for(i = n - ; i >= ; i -- )
sa[ -- c[x[y[i]]]] = y[i];
//根据 sa 和 x 数组计算新的 x 数组
swap(x, y);
p = ;
x[sa[]] = ;
for(i = ; i < n; i++)
x[sa[i]] = cmp(y, sa[i - ], sa[i], j) ? p - : p++;
if(p >= n)
break;
m = p;//下次基数排序的最大值
}
int k = ;
n -- ;
for(i = ; i <= n; i++)
rk[sa[i]] = i;
for(i = ; i < n; i++) {
if(k) k -- ;
j = sa[rk[i] - ];
while(str[i + k] == str[j + k])
k++;
height[rk[i]] = k;
}
}
int rk[MAXN], height[MAXN], sa[MAXN], in[MAXN];
char s[MAXN], t[MAXN];
vector<int>vec;
vector<int>:: iterator it; int main() {
int T;scanf("%d", &T);
for(int ncase = ; ncase <= T; ncase ++){
scanf("%s%s",t, s);
int n = strlen(s);
vec.clear();
for(int i = ; i < n; i ++){
in[i] = s[i];
if(s[i] == t[]) vec.push_back(i);
}in[n] = ;
da(in, sa, rk, height, n, );
ll ans = ;
for(int i = ; i <= n; i ++){
it = lower_bound(vec.begin(), vec.end(), sa[i]);
printf("%d %d %d\n", sa[i], *it, height[i]);
if(it == vec.end()) continue;
int tmp = n - sa[i] - max(height[i], *it - sa[i]);
if(tmp > ) ans += tmp;
printf("Case #%d: %lld\n",ncase, ans);
}printf("\n");
}
return ;
}

Substring (后缀数组 + 计数)的更多相关文章

  1. POJ3693 Maximum repetition substring [后缀数组 ST表]

    Maximum repetition substring Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 9458   Acc ...

  2. HDU 5769 Substring 后缀数组

    Substring Problem Description ?? is practicing his program skill, and now he is given a string, he h ...

  3. hdu_1403_Longest Common Substring(后缀数组的应用)

    题目链接:hdu_1403_Longest Common Substring 题意: 给你两个字符串,然你找最长的公共子串 题解: 后缀数组的经典应用,要找两个字符串的公共子串,那么就相当于找两个串的 ...

  4. POJ3693 Maximum repetition substring 后缀数组

    POJ - 3693 Maximum repetition substring 题意 输入一个串,求重复次数最多的连续重复字串,如果有次数相同的,则输出字典序最小的 Sample input ccab ...

  5. 2016多校联合训练4 F - Substring 后缀数组

    Description ?? is practicing his program skill, and now he is given a string, he has to calculate th ...

  6. hdu 5769 Substring 后缀数组 + KMP

    http://acm.hdu.edu.cn/showproblem.php?pid=5769 题意:在S串中找出X串出现的不同子串的数目? 其中1 <= |S| < $10^5$ 官方题解 ...

  7. HDU 5679 Substring 后缀数组判重

    题意:求母串中有多少不同的包含x字符的子串 分析:(首先奉上FZU官方题解) 上面那个题就是SPOJ694 ,其实这两个题一样,原理每次从小到大扫后缀sa数组,加上新的当前后缀的若干前缀,再减去重复的 ...

  8. hdu 1403 Longest Common Substring 后缀数组 模板题

    题目链接 题意 问两个字符串的最长公共子串. 思路 加一个特殊字符然后拼接起来,求得后缀数组与\(height\)数组.扫描一遍即得答案,注意判断起始点是否分别在两个串内. Code #include ...

  9. POJ 3693 Maximum repetition substring ——后缀数组

    重复次数最多的字串,我们可以枚举循环节的长度. 然后正反两次LCP,然后发现如果长度%L有剩余的情况时,答案是在一个区间内的. 所以需要找到区间内最小的rk值. 两个后缀数组,四个ST表,$\Thet ...

随机推荐

  1. tensorflow简单记录summary方法

    虽然tf官方希望用户把 train , val 程序分开写,但实际开发中,明显写在一起比较简单舒服,但在保存数据到 summary 时, val 部分和 train 部分不太一样,会有一些问题,下面讨 ...

  2. 乾坤合一~Linux设备驱动之块设备驱动

    1. 题外话 在蜕变成蝶的一系列学习当中,我们已经掌握了大部分Linux驱动的知识,在乾坤合一的分享当中,以综合实例为主要讲解,在一个月的蜕茧成蝶的学习探索当中,觉得数据结构,指针,链表等等占据了代码 ...

  3. eclipse 查看变量或方法在什么地方被调用的快捷键和快速显示方法入参提示信息

    为了用eclipse写代码的时候,不用方向键移动光标,可以设置如下快捷键, Line Up:Alt+i   光标向上 Line Down:Alt+k 光标向下 Line Start:Alt+u  移到 ...

  4. nginx 内置变量

    http://blog.sina.com.cn/s/articlelist_1834459124_1_1.html  nginx内置变量杂谈 http://nginx.org/en/docs/http ...

  5. 常见机试题分析Java版

    1. 操作系统任务分为系统任务和用户任务两种.其中,系统任务的优先级<50,用户任务的优先级>=50且<=255.优先级大于255的为非法任务,应予以剔除.现有一任务队列task[] ...

  6. 4e00~9fa5

    <script> function init (){ var begin = 0x4e00; var end = 0x9fa5; var str = ''; for(var a=begin ...

  7. Sublime Text 插件之HTML-CSS-JS Prettify—格式化HTML CSS JS与显示函数列表

    插件名称:HTML-CSS-JS Prettify 安装步骤: 1.ctrl + shift + p 打开控制台2.输入install package,选择install package(如果已经安装 ...

  8. 阅历>感悟

    1.强扭的瓜不甜.在招聘的时候,面试官看不上你,你也不用赖着要去,你去能干好工作吗?面试官通常比你更清楚这个是事情.在比如谈恋爱,姑娘有更好的目标,不喜欢你了,决定离开你了,你再怎么挽留都是没意义的, ...

  9. hdu 2098 分拆素数和(素数)

    分拆素数和 Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submi ...

  10. oracle编码转换:AL32UTF8->ZHS16GBK

    --修改Oracle数据库字符集为utf-8: SQL>conn / as sysdba; SQL>shutdown immediate; SQL>startup mount; SQ ...