Cod sursa(job #512310)

Utilizator antoanelaAntoanela Siminiuc antoanela Data 13 decembrie 2010 19:47:33
Problema Potrivirea sirurilor Scor 78
Compilator cpp Status done
Runda Arhiva educationala Marime 1 kb
#include <cstdio>
#include <cstring>
#define p 79
#define mod 100007
#define mod2 100021
#define nmax 2000010

int na, nb, hasha, hash, nr, sol[1010];
char a[nmax], b[nmax];


int main()
{
	freopen("strmatch.in","r",stdin);
	freopen("strmatch.out","w",stdout);
	scanf("%s %s",&a,&b);
	na=strlen(a);
	nb=strlen(b);
	if (na>nb)
	{
		printf("0\n");
		return 0;
	}
	int hasha=0, hash=0, hash2=0, hasha2=0;
	int i, pn=1, pn2=1;
	for (i=0; i<na; i++)
	{
		hasha=(hasha*p+a[i])%mod;
		hasha2=(hasha2*p+a[i])%mod2;
		if (i>0) pn=(pn*p)%mod;
		if (i>0) pn2=(pn2*p)%mod2;
		hash=(hash*p+b[i])%mod;
		hash2=(hash2*p+b[i])%mod2;
	}
	if (hasha==hash && hasha2==hash2) sol[++nr]=1;
	for (i=na; i<nb; i++)
	{
		hash=((hash-(b[i-na]*pn)%mod +mod)*p+b[i])%mod;
		hash2=((hash2-(b[i-na]*pn2)%mod2 +mod2)*p+b[i])%mod2;
		if (hash==hasha && hasha2==hash2) 
		{
			nr++;
			if (nr<=1000) sol[nr]=i-na+1;
		}
	}
	printf("%d\n",nr);
	if (nr>1000) nr=1000;
	for (i=1; i<=nr; i++) printf("%d ",sol[i]);
}