0%

c-找不同的子串

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
#include <cstdio>
#include <iostream>
#include <string>
#include <algorithm>
#include <cmath>
#include <cstring>

using namespace std;
#define MAXN 1011

int n,k;//n=strlen(s);

int Rank[MAXN];
int tmp[MAXN];
char s[MAXN];
int lcp[MAXN],sa[MAXN];

/*使用Rank对sa排序*/
bool cmpSa(int i, int j)
{
if(Rank[i] != Rank[j])return Rank[i] < Rank[j];
else
{ /*以下的Rank[t]。已经是以t开头长度小于等于k/2的,
sa[i]的名次。仅仅是以i开头的后缀。而长度不同*/
int ri = i+k <=n? Rank[i+k]:-1;
int rj = j+k <= n ? Rank[j+k]:-1;
return ri <rj;
}
}

/*计算SA*/
void consa()
{
/*n=strlen(s); 必要时注明*/
/*初始化sa和rank保证两点
1、Rank[i]表示下标为i的是第几大,必须表示出相对大小。能够直接用字符代表其大小
2、sa[1...n]值为1..n*/
for(int i=0;i<=n;i++){
sa[i]=i;Rank[i] = i < n?s[i]:-1;
}

/*利用长度为k的字符串对长度为2*k的字符串排序*/
for(k=1;k<=n;k*=2)/*注意此代码中k是全局变量 别乱用,循环必须从1開始,由于0*2=0*/
{
sort(sa,sa+n+1,cmpSa);
tmp[sa[0]] = 0; /*此时tmp仅仅是暂存rank*/
for(int i=1;i<=n;i++){
tmp[sa[i]] = tmp[sa[i-1]] +(cmpSa(sa[i-1],sa[i])?1:0);
/*这一句非常关键,等号右側的sa[i]在此循环里表示第i大的长度小于等于k/2的字符串,
从而求出第i大的长度小于等于k的字符串的sa[i]*/
}
for(int i=0;i<=n;i++){
Rank[i] = tmp[i];
}
}
}

void construct_lcp()
{
//n=strlen(s);
for(int i=0; i<=n; i++)Rank[sa[i]]=i;

int h=0;
lcp[0]=0;
for(int i=0;i<n;i++)
{
int j=sa[Rank[i]-1];

if(h>0)h--;
for(; j+h<n && i+h<n; h++)
{
if(s[j+h]!=s[i+h])break;
}
lcp[Rank[i]-1]=h;
}
}

int main()
{
int t,ans;
scanf("%d",&t);
while(t--)
{
ans=0;
scanf("%s",s);
n=strlen(s);
consa();
construct_lcp();
for(int i=1;i<=n;i++)
{
ans+=n-sa[i]-lcp[i-1];
}
printf("%d\n",ans);
}
return 0;

return 0;
}
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
#include<iostream>
#define maxn 1024
#include<set>
using namespace std;
long long sum=0,f=1;
set<string> sr[maxn];
int main()
{
string s1;
cin>>s1;
for(int i=0;i<s1.length();i++)
{
for(int j=1;j<=s1.length();j++)
{
string s2=s1.substr(i,j);
if(s2.length()==j)
sr[j-1].insert(s2);
}
}
for(int i=0;i<s1.length();i++)
sum+=sr[i].size();
cout<<sum;
/*
set<string>::iterator int1=sr[0].begin();
set<string>::iterator int2=sr[0].end();

for(;int1!=int2;int1++)
cout<<*int1<<" ";
return 0;
*/
}