Given a text txt[0..n-1] and a pattern pat[0..m-1], write a function search(char pat[], char txt[]) that prints all occurrences of pat[] in txt[]. You may assume that n > m.

Examples:

Input:  txt[] = "THIS IS A TEST TEXT"
        pat[] = "TEST"
Output: Pattern found at index 10

Input:  txt[] =  "AABAACAADAABAABA"
        pat[] =  "AABA"
Output: Pattern found at index 0
        Pattern found at index 9
        Pattern found at index 12
C++ programming for Searching for Patterns Set 2 KMP Algorithm

The Naive String Matching algorithm slides the pattern one by one. After each slide, it one by one checks characters at the current shift and if all characters match then prints the match.

[ad type=”banner”] Like the Naive Algorithm, Rabin-Karp algorithm also slides the pattern one by one. But unlike the Naive algorithm, Rabin Karp algorithm matches the hash value of the pattern with the hash value of current substring of text, and if the hash values match then only it starts matching individual characters. So Rabin Karp algorithm needs to calculate hash values for following strings.

1) Pattern itself.
2) All the substrings of text of length m.

Since we need to efficiently calculate hash values for all the substrings of size m of text, we must have a hash function which has following property.
Hash at the next shift must be efficiently computable from the current hash value and next character in text or we can say hash(txt[s+1 .. s+m]) must be efficiently computable from hash(txt[s .. s+m-1]) and txt[s+m] i.e., hash(txt[s+1 .. s+m])= rehash(txt[s+m], hash(txt[s .. s+m-1]) and rehash must be O(1) operation.

The hash function suggested by Rabin and Karp calculates an integer value. The integer value for a string is numeric value of a string. For example, if all possible characters are from 1 to 10, the numeric value of “122” will be 122. The number of possible characters is higher than 10 (256 in general) and pattern length can be large. So the numeric values cannot be practically stored as an integer. Therefore, the numeric value is calculated using modular arithmetic to make sure that the hash values can be stored in an integer variable (can fit in memory words). To do rehashing, we need to take off the most significant digit and add the new least significant digit for in hash value. Rehashing is done using the following formula.

[ad type=”banner”]

hash( txt[s+1 .. s+m] ) = d ( hash( txt[s .. s+m-1]) – txt[s]*h ) + txt[s + m] ) mod q

hash( txt[s .. s+m-1] ) : Hash value at shift s.
hash( txt[s+1 .. s+m] ) : Hash value at next shift (or shift s+1)
d: Number of characters in the alphabet
q: A prime number
h: d^(m-1)

[pastacode lang=”cpp” manual=”%2F*%20Following%20program%20is%20a%20C%2B%2B%20implementation%20of%20Rabin%20Karp%0AAlgorithm%20given%20in%20the%20CLRS%20book%20*%2F%0A%23include%3Cstdio.h%3E%0A%23include%3Cstring.h%3E%0A%20%0A%2F%2F%20d%20is%20the%20number%20of%20characters%20in%20input%20alphabet%0A%23define%20d%20256%0A%20%0A%2F*%20pat%20-%3E%20pattern%0A%20%20%20%20txt%20-%3E%20text%0A%20%20%20%20q%20-%3E%20A%20prime%20number%0A*%2F%0Avoid%20search(char%20pat%5B%5D%2C%20char%20txt%5B%5D%2C%20int%20q)%0A%7B%0A%20%20%20%20int%20M%20%3D%20strlen(pat)%3B%0A%20%20%20%20int%20N%20%3D%20strlen(txt)%3B%0A%20%20%20%20int%20i%2C%20j%3B%0A%20%20%20%20int%20p%20%3D%200%3B%20%2F%2F%20hash%20value%20for%20pattern%0A%20%20%20%20int%20t%20%3D%200%3B%20%2F%2F%20hash%20value%20for%20txt%0A%20%20%20%20int%20h%20%3D%201%3B%0A%20%0A%20%20%20%20%2F%2F%20The%20value%20of%20h%20would%20be%20%22pow(d%2C%20M-1)%25q%22%0A%20%20%20%20for%20(i%20%3D%200%3B%20i%20%3C%20M-1%3B%20i%2B%2B)%0A%20%20%20%20%20%20%20%20h%20%3D%20(h*d)%25q%3B%0A%20%0A%20%20%20%20%2F%2F%20Calculate%20the%20hash%20value%20of%20pattern%20and%20first%0A%20%20%20%20%2F%2F%20window%20of%20text%0A%20%20%20%20for%20(i%20%3D%200%3B%20i%20%3C%20M%3B%20i%2B%2B)%0A%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20p%20%3D%20(d*p%20%2B%20pat%5Bi%5D)%25q%3B%0A%20%20%20%20%20%20%20%20t%20%3D%20(d*t%20%2B%20txt%5Bi%5D)%25q%3B%0A%20%20%20%20%7D%0A%20%0A%20%20%20%20%2F%2F%20Slide%20the%20pattern%20over%20text%20one%20by%20one%0A%20%20%20%20for%20(i%20%3D%200%3B%20i%20%3C%3D%20N%20-%20M%3B%20i%2B%2B)%0A%20%20%20%20%7B%0A%20%0A%20%20%20%20%20%20%20%20%2F%2F%20Check%20the%20hash%20values%20of%20current%20window%20of%20text%0A%20%20%20%20%20%20%20%20%2F%2F%20and%20pattern.%20If%20the%20hash%20values%20match%20then%20only%0A%20%20%20%20%20%20%20%20%2F%2F%20check%20for%20characters%20on%20by%20one%0A%20%20%20%20%20%20%20%20if%20(%20p%20%3D%3D%20t%20)%0A%20%20%20%20%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%2F*%20Check%20for%20characters%20one%20by%20one%20*%2F%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20(j%20%3D%200%3B%20j%20%3C%20M%3B%20j%2B%2B)%0A%20%20%20%20%20%20%20%20%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20(txt%5Bi%2Bj%5D%20!%3D%20pat%5Bj%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20break%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%7D%0A%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%2F%2F%20if%20p%20%3D%3D%20t%20and%20pat%5B0…M-1%5D%20%3D%20txt%5Bi%2C%20i%2B1%2C%20…i%2BM-1%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20(j%20%3D%3D%20M)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20printf(%22Pattern%20found%20at%20index%20%25d%20%5Cn%22%2C%20i)%3B%0A%20%20%20%20%20%20%20%20%7D%0A%20%0A%20%20%20%20%20%20%20%20%2F%2F%20Calculate%20hash%20value%20for%20next%20window%20of%20text%3A%20Remove%0A%20%20%20%20%20%20%20%20%2F%2F%20leading%20digit%2C%20add%20trailing%20digit%0A%20%20%20%20%20%20%20%20if%20(%20i%20%3C%20N-M%20)%0A%20%20%20%20%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20t%20%3D%20(d*(t%20-%20txt%5Bi%5D*h)%20%2B%20txt%5Bi%2BM%5D)%25q%3B%0A%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%2F%2F%20We%20might%20get%20negative%20value%20of%20t%2C%20converting%20it%0A%20%20%20%20%20%20%20%20%20%20%20%20%2F%2F%20to%20positive%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20(t%20%3C%200)%0A%20%20%20%20%20%20%20%20%20%20%20%20t%20%3D%20(t%20%2B%20q)%3B%0A%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20%7D%0A%7D” message=”C++” highlight=”” provider=”manual”/]

Output:
Pattern found at index 0
Pattern found at index 10
The average and best case running time of the Rabin-Karp algorithm is O(n+m), but its worst-case time is O(nm). Worst case of Rabin-Karp algorithm occurs when all characters of pattern and text are same as the hash values of all the substrings of txt[] match with hash value of pat[]. For example pat[] = “AAA” and txt[] = “AAAAAAA”.

Please write comments if you find anything incorrect, or you want to share more information about the topic discussed above.

[ad type=”banner”]