Pagini recente » Cod sursa (job #2215507) | Cod sursa (job #2254890) | Cod sursa (job #2386643) | Cod sursa (job #290183) | Cod sursa (job #2619232)
#define _CRT_SECURE_NO_WARNINGS
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
void computeLPSArray(char* pat, int M, int* lps)
{
// length of the previous longest prefix suffix
int len = 0;
lps[0] = 0; // lps[0] is always 0
// the loop calculates lps[i] for i = 1 to M-1
int i = 1;
while (i < M) {
if (pat[i] == pat[len]) {
len++;
lps[i] = len;
i++;
}
else // (pat[i] != pat[len])
{
// This is tricky. Consider the example.
// AAACAAAA and i = 7. The idea is similar
// to search step.
if (len != 0) {
len = lps[len - 1];
// Also, note that we do not increment
// i here
}
else // if (len == 0)
{
lps[i] = 0;
i++;
}
}
}
}
int v[1000001];
int r;
void KMPSearch(char* pat, char* txt)
{
int M = strlen(pat);
int N = strlen(txt);
// create lps[] that will hold the longest prefix suffix
// values for pattern
int *lps=(int*)malloc(sizeof(int)*(M+1));
// Preprocess the pattern (calculate lps[] array)
computeLPSArray(pat, M, lps);
int i = 0; // index for txt[]
int j = 0; // index for pat[]
while (i < N) {
if (pat[j] == txt[i]) {
j++;
i++;
}
if (j == M) {
v[r++] = i - j;
j = lps[j - 1];
}
// mismatch after j matches
else if (i < N && pat[j] != txt[i]) {
// Do not match lps[0..lps[j-1]] characters,
// they will match anyway
if (j != 0)
j = lps[j - 1];
else
i = i + 1;
}
}
}
// Fills lps[] for given patttern pat[0..M-1]
int main()
{
char* a = (char*)malloc(20000000);
char* b = (char*)malloc(20000000);
FILE* fin = fopen("strmatch.in", "r");
FILE* fout = fopen("strmatch.out", "w");
fscanf(fin, "%s%s", a,b);
KMPSearch(a, b);
fprintf(fout, "%d\n", r);
for (int i = 0; i < r; i++)
fprintf(fout, "%d ", v[i]);
}