/* Following program is a C++ implementation of Rabin Karp
Algorithm given in the CLRS book */
#include <bits/stdc++.h>
using namespace std;
// d is the number of characters in the input alphabet
#define d 256
/* pat -> pattern
txt -> text
q -> A prime number
*/
void search(char pat[], char txt[], int q)
{
int M = strlen(pat);
int N = strlen(txt);
int i, j;
int p = 0; // hash value for pattern
int t = 0; // hash value for txt
int h = 1;
// The value of h would be "pow(d, M-1)%q"
for (i = 0; i < M - 1; i++)
h = (h * d) % q;
// Calculate the hash value of pattern and first
// window of text
for (i = 0; i < M; i++)
{
p = (d * p + pat[i]) % q;
t = (d * t + txt[i]) % q;
}
// Slide the pattern over text one by one
for (i = 0; i <= N - M; i++)
{
// Check the hash values of current window of text
// and pattern. If the hash values match then only
// check for characters on by one
if ( p == t )
{
/* Check for characters one by one */
for (j = 0; j < M; j++)
{
if (txt[i+j] != pat[j])
break;
}
// if p == t and pat[0...M-1] = txt[i, i+1, ...i+M-1]
if (j == M)
cout<<"Pattern found at index "<< i<<endl;
}
// Calculate hash value for next window of text: Remove
// leading digit, add trailing digit
if ( i < N-M )
{
t = (d*(t - txt[i]*h) + txt[i+M])%q;
// We might get negative value of t, converting it
// to positive
if (t < 0)
t = (t + q);
}
}
}
/* Driver code */
int main()
{
char txt[] = "GEEKS FOR GEEKS";
char pat[] = "GEEK";
int q = 101; // A prime number
search(pat, txt, q);
return 0;
}
Ci8qIEZvbGxvd2luZyBwcm9ncmFtIGlzIGEgQysrIGltcGxlbWVudGF0aW9uIG9mIFJhYmluIEthcnAgIApBbGdvcml0aG0gZ2l2ZW4gaW4gdGhlIENMUlMgYm9vayAqLwojaW5jbHVkZSA8Yml0cy9zdGRjKysuaD4gCnVzaW5nIG5hbWVzcGFjZSBzdGQ7IAogIAovLyBkIGlzIHRoZSBudW1iZXIgb2YgY2hhcmFjdGVycyBpbiB0aGUgaW5wdXQgYWxwaGFiZXQgIAojZGVmaW5lIGQgMjU2ICAKICAKLyogcGF0IC0+IHBhdHRlcm4gIAogICAgdHh0IC0+IHRleHQgIAogICAgcSAtPiBBIHByaW1lIG51bWJlciAgCiovCnZvaWQgc2VhcmNoKGNoYXIgcGF0W10sIGNoYXIgdHh0W10sIGludCBxKSAgCnsgIAogICAgaW50IE0gPSBzdHJsZW4ocGF0KTsgIAogICAgaW50IE4gPSBzdHJsZW4odHh0KTsgIAogICAgaW50IGksIGo7ICAKICAgIGludCBwID0gMDsgLy8gaGFzaCB2YWx1ZSBmb3IgcGF0dGVybiAgCiAgICBpbnQgdCA9IDA7IC8vIGhhc2ggdmFsdWUgZm9yIHR4dCAgCiAgICBpbnQgaCA9IDE7ICAKICAKICAgIC8vIFRoZSB2YWx1ZSBvZiBoIHdvdWxkIGJlICJwb3coZCwgTS0xKSVxIiAgCiAgICBmb3IgKGkgPSAwOyBpIDwgTSAtIDE7IGkrKykgIAogICAgICAgIGggPSAoaCAqIGQpICUgcTsgIAogIAogICAgLy8gQ2FsY3VsYXRlIHRoZSBoYXNoIHZhbHVlIG9mIHBhdHRlcm4gYW5kIGZpcnN0ICAKICAgIC8vIHdpbmRvdyBvZiB0ZXh0ICAKICAgIGZvciAoaSA9IDA7IGkgPCBNOyBpKyspICAKICAgIHsgIAogICAgICAgIHAgPSAoZCAqIHAgKyBwYXRbaV0pICUgcTsgIAogICAgICAgIHQgPSAoZCAqIHQgKyB0eHRbaV0pICUgcTsgIAogICAgfSAgCiAgCiAgICAvLyBTbGlkZSB0aGUgcGF0dGVybiBvdmVyIHRleHQgb25lIGJ5IG9uZSAgCiAgICBmb3IgKGkgPSAwOyBpIDw9IE4gLSBNOyBpKyspICAKICAgIHsgIAogIAogICAgICAgIC8vIENoZWNrIHRoZSBoYXNoIHZhbHVlcyBvZiBjdXJyZW50IHdpbmRvdyBvZiB0ZXh0ICAKICAgICAgICAvLyBhbmQgcGF0dGVybi4gSWYgdGhlIGhhc2ggdmFsdWVzIG1hdGNoIHRoZW4gb25seSAgCiAgICAgICAgLy8gY2hlY2sgZm9yIGNoYXJhY3RlcnMgb24gYnkgb25lICAKICAgICAgICBpZiAoIHAgPT0gdCApICAKICAgICAgICB7ICAKICAgICAgICAgICAgLyogQ2hlY2sgZm9yIGNoYXJhY3RlcnMgb25lIGJ5IG9uZSAqLwogICAgICAgICAgICBmb3IgKGogPSAwOyBqIDwgTTsgaisrKSAgCiAgICAgICAgICAgIHsgIAogICAgICAgICAgICAgICAgaWYgKHR4dFtpK2pdICE9IHBhdFtqXSkgIAogICAgICAgICAgICAgICAgICAgIGJyZWFrOyAgCiAgICAgICAgICAgIH0gIAogIAogICAgICAgICAgICAvLyBpZiBwID09IHQgYW5kIHBhdFswLi4uTS0xXSA9IHR4dFtpLCBpKzEsIC4uLmkrTS0xXSAgCiAgICAgICAgICAgIGlmIChqID09IE0pICAKICAgICAgICAgICAgICAgIGNvdXQ8PCJQYXR0ZXJuIGZvdW5kIGF0IGluZGV4ICI8PCBpPDxlbmRsOyAgCiAgICAgICAgfSAgCiAgCiAgICAgICAgLy8gQ2FsY3VsYXRlIGhhc2ggdmFsdWUgZm9yIG5leHQgd2luZG93IG9mIHRleHQ6IFJlbW92ZSAgCiAgICAgICAgLy8gbGVhZGluZyBkaWdpdCwgYWRkIHRyYWlsaW5nIGRpZ2l0ICAKICAgICAgICBpZiAoIGkgPCBOLU0gKSAgCiAgICAgICAgeyAgCiAgICAgICAgICAgIHQgPSAoZCoodCAtIHR4dFtpXSpoKSArIHR4dFtpK01dKSVxOyAgCiAgCiAgICAgICAgICAgIC8vIFdlIG1pZ2h0IGdldCBuZWdhdGl2ZSB2YWx1ZSBvZiB0LCBjb252ZXJ0aW5nIGl0ICAKICAgICAgICAgICAgLy8gdG8gcG9zaXRpdmUgIAogICAgICAgICAgICBpZiAodCA8IDApICAKICAgICAgICAgICAgdCA9ICh0ICsgcSk7ICAKICAgICAgICB9ICAKICAgIH0gIAp9ICAKICAKLyogRHJpdmVyIGNvZGUgKi8KaW50IG1haW4oKSAgCnsgIAogICAgY2hhciB0eHRbXSA9ICJHRUVLUyBGT1IgR0VFS1MiOyAgCiAgICBjaGFyIHBhdFtdID0gIkdFRUsiOyAgCiAgICBpbnQgcSA9IDEwMTsgLy8gQSBwcmltZSBudW1iZXIgIAogICAgc2VhcmNoKHBhdCwgdHh0LCBxKTsgIAogICAgcmV0dXJuIDA7ICAKfSAgCiAg