Cod sursa(job #911388)

Utilizator DEYDEY2Tudorica Andrei DEYDEY2 Data 11 martie 2013 16:27:36
Problema Potrivirea sirurilor Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.35 kb
#include <fstream>
#include <string.h>
#define MAXN 2000001
#define P 73
#define MOD1 100007
#define MOD2 100021
using namespace std;
ifstream f("strmatch.in");
ofstream g("strmatch.out");
char A[MAXN], B[MAXN],match[MAXN];
int NA,NB,hashA1, hashA2, P1, P2;
int main()
{
    f>>A>>B;
    NA=strlen(A);
    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){g<<0;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++;
    }
    g<<Nr<<'\n';
    Nr = 0;
    for (int i = 0; i < NB && Nr < 1000; i++)
        if (match[i])
        {
            Nr++;
            g<<i<<' ';
        }
    f.close();
    g.close();
    return 0;
}