Pagini recente » Cod sursa (job #1398095) | Cod sursa (job #2330134) | Cod sursa (job #2556425) | Cod sursa (job #3183794) | Cod sursa (job #1328177)
// Rabin Karp (2 bases)
#include <fstream>
#include <string.h>
#include <algorithm>
#include <iostream>
using namespace std;
fstream fin("strmatch.in", ios::in);
fstream fout("strmatch.out", ios::out);
#define MAX 2000005
#define MOD1 100007
#define MOD2 100021
#define P 73
char A[MAX], B[MAX];
int NA, NB;
int hashA1, hashA2, P1, P2, sol[1001];
char match[MAX];
int main()
{
int P1=1, P2=1;
fin.getline(A, MAX); NA=strlen(A);
fin.getline(B, MAX); NB=strlen(B);
P1 = P2 = 1;
hashA1 = hashA2 = 0;
for (int i = 0; i < NA; i++)
{
hashA1 = (hashA1 * P + A[i]) % MOD1;
hashA2 = (hashA2 * P + A[i]) % MOD2;
if (i != 0)
P1 = (P1 * P) % MOD1,
P2 = (P2 * P) % MOD2;
}
if (NA > NB)
{
printf("0\n");
return 0;
}
int hash1 = 0, hash2 = 0;
for (int i = 0; i < NA; i++)
hash1 = (hash1 * P + B[i]) % MOD1,
hash2 = (hash2 * P + B[i]) % MOD2;
int Nr = 0;
if (hash1 == hashA1 && hash2 == hashA2)
match[0] = 1, Nr++;
for (int i = NA; i < NB; i++)
{
hash1 = ((hash1 - (B[i - NA] * P1) % MOD1 + MOD1) * P + B[i]) % MOD1;
hash2 = ((hash2 - (B[i - NA] * P2) % MOD2 + MOD2) * P + B[i]) % MOD2;
if (hash1 == hashA1 && hash2 == hashA2)
match[ i - NA + 1 ] = 1, Nr++;
}
fout<<Nr<<"\n";
Nr = 0;
for (int i = 0; i < NB && Nr < 1000; i++)
if (match[i])
{
Nr++;
fout<<i<<" ";
}
fout<<"\n";
fin.close();
fout.close();
return 0;
}