Substring【HDU-5769】【后缀数组】
程序员文章站
2022-04-17 14:07:30
...
题目链接
题意:要让我们求包含字符X的本质不同的子串的个数。
思路:我们可以求以每个点作为子串的首字符时候的造成的贡献,那么此时首先就是需要减去前面的height,还要再去找到后面的字符X的最近的位置,可以是自己本身就是字符X,然后就可以用这两个来进行计算贡献了。
#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <limits>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <bitset>
#include <unordered_map>
#include <unordered_set>
#define lowbit(x) ( x&(-x) )
#define pi 3.141592653589793
#define e 2.718281828459045
#define INF 0x3f3f3f3f
#define HalF (l + r)>>1
#define lsn rt<<1
#define rsn rt<<1|1
#define Lson lsn, l, mid
#define Rson rsn, mid+1, r
#define QL Lson, ql, qr
#define QR Rson, ql, qr
#define myself rt, l, r
using namespace std;
typedef unsigned long long ull;
typedef unsigned int uit;
typedef long long ll;
const int maxN = 1e5 + 7;
struct SA
{
int n, m;
int s[maxN];
int y[maxN], x[maxN], c[maxN], sa[maxN], rk[maxN], height[maxN];
inline void get_SA()
{
for(int i=1; i<=m; i++) c[i] = 0; //桶的初始化
for(int i=1; i<=n; i++) ++c[x[i] = s[i]];
for(int i=2; i<=m; i++) c[i] += c[i - 1]; //利用差分前缀和的思想知道每个关键字最多是在第几名
for(int i=n; i>=1; i--) sa[c[x[i]]--] = i;
for(int k=1; k<=n; k<<=1)
{
int num = 0;
for(int i=n - k + 1; i<=n; i++) y[++num] = i;
for(int i=1; i<=n; i++) if(sa[i] > k) y[++num] = sa[i] - k; //是否可以作为第二关键字
for(int i=1; i<=m; i++) c[i] = 0;
for(int i=1; i<=n; i++) c[x[i]]++; //因为上一次循环已经求出这次的第一关键字了
for(int i=2; i<=m; i++) c[i] += c[i - 1];
for(int i=n; i>=1; i--) //在同一第一关键字下,按第二关键字来排
{
sa[c[x[y[i]]]--] = y[i];
y[i] = 0;
}
swap(x, y);
x[sa[1]] = 1; num = 1;
for(int i=2; i<=n; i++)
{
x[sa[i]] = (y[sa[i]] == y[sa[i - 1]] && y[sa[i] + k] == y[sa[i - 1] + k]) ? num : ++num;
}
if(num == n) break;
m = num;
}
}
inline void get_height()
{
int k = 0;
for(int i=1; i<=n; i++) rk[sa[i]] = i;
for(int i=1; i<=n; i++)
{
if(rk[i] == 1) continue; //第一名的height为0
if(k) k--; //height[i] >= height[i - 1] - 1
int j = sa[rk[i] - 1];
while(j + k <= n && i + k <= n && s[i + k] == s[j + k]) k++;
height[rk[i]] = k;
}
}
inline void clear()
{
n = 0; m = 200;
}
} sa;
char ch[5], s[maxN];
int nex[maxN];
ll ans;
int main()
{
int T; scanf("%d", &T);
for(int Cas=1; Cas<=T; Cas++)
{
scanf("%s%s", ch, s);
sa.clear();
int N = (int)strlen(s);
for(int i=0; i<N; i++)
{
sa.s[++sa.n] = s[i];
}
for(int i=N - 1, pos = N + 1; i>=0; i--)
{
if(s[i] == ch[0])
{
pos = i + 1;
}
nex[i + 1] = pos;
}
sa.s[++sa.n] = 'a' + 30;
sa.get_SA();
sa.get_height();
// for(int i=1; i<=N; i++) printf("%d%c", sa.sa[i], i == N ? '\n' : ' ');
// for(int i=1; i<=N; i++) printf("%d%c", sa.rk[i], i == N ? '\n' : ' ');
ans = 0;
for(int i=1, ith, beg_pos; i<=N; i++)
{
ith = sa.sa[i];
beg_pos = ith + sa.height[i];
if(nex[ith] == N + 1) continue;
if(beg_pos >= nex[ith])
{
ans += N - beg_pos + 1;
}
else ans += N - nex[ith] + 1;
}
printf("Case #%d: %lld\n", Cas, ans);
}
return 0;
}
/*
1
a
aabaa
ans:10
*/
下一篇: 机场柜台职员
推荐阅读
-
BZOJ2882: 工艺(后缀数组)
-
POJ1743 Musical Theme(后缀数组 二分)
-
【总结】后缀数组
-
洛谷P5108 仰望半月的夜空(后缀数组)
-
BZOJ4698: Sdoi2008 Sandy的卡片(后缀数组 二分)
-
洛谷P3763 [TJOI2017]DNA(后缀数组 RMQ)
-
BZOJ4516: [Sdoi2016]生成魔咒(后缀数组 set RMQ)
-
SPOJ1811 LCS - Longest Common Substring(后缀自动机)
-
文件夹下(包含子文件夹和文件),取文件夹和子文件夹下所有后缀为JPG的文件的,路径和文件名 ,把路径和文件名放在数组中
-
BZOJ2882: 工艺(后缀数组)