https://www.hackerrank.com/challenges/the-grid-search/forum

今天碰見這題，看見難度是Moderate，覺得應(yīng)該能半小時內(nèi)搞定。

讀完題目發(fā)現(xiàn)是純粹的一道子矩陣匹配問題，想想自己以前沒做過，肯定能學(xué)到新算法，于是就開搞了。

于是上網(wǎng)搜到了Rabin-Karp算法，一種基于hashing的模式匹配算法。盡管連一維的我也沒寫過，但看了思想以后覺得推廣到二維應(yīng)該也不會很難。

于是有了以下代碼，原理就是計算子矩陣的hash key。以hash key的比較代替了子矩陣的比較，這樣可以首先排除掉hash key不相等的子矩陣。

對于hash key相等的，再用樸素方法判斷子矩陣是否相等。

為什么最后還是要判斷子矩陣是否相等呢？因為hash key可能存在碰撞，即使概率不大，為了保證正確性也需要進(jìn)行檢查。

學(xué)習(xí)Rabin-Karp算法的資料在此：

http://blog.sina.com.cn/s/blog_6a09b5a70100nhnr.html

?

思路雖簡單，代碼寫起來卻各種bug，最終我花了不下兩個鐘頭才搞定。Hackerrank果然是給hacker玩的，我這水平在上面真是舉步維艱。

不過，這么搞下來，倒是有了實實在在的收獲，如果學(xué)習(xí)算法能一直保持這種節(jié)奏就好了。

下面是AC的代碼，時空復(fù)雜度均為O(N ^ 2)：

1 # 2D Rabin-Karp Algorithm 2 import re 3 4 MOD = 10 ** 9 + 7 5 6 def get2DMatrix(n, m): 7 a = [[0 for j in xrange(m)] for i in xrange(n)] 8 return a 9 10 def calcHash(a, nn, mm): 11 n = len(a) 12 m = len(a[0]) 13 14 b = 1 15 for i in xrange(mm): 16 b = b * 10 % MOD 17 b2 = 1 18 for i in xrange(nn): 19 b2 = b2 * b % MOD 20 21 h = get2DMatrix(n, m) 22 for i in xrange(n): 23 val = 0 24 for j in xrange(m): 25 val = (val * 10 + a[i][j]) % MOD 26 if j >= mm: 27 val = (val + a[i][j - mm] * (MOD - b)) % MOD 28 h[i][j] = val 29 30 h2 = get2DMatrix(n, m) 31 h2[0] = h[0][:] 32 for i in xrange(1, n): 33 for j in xrange(m): 34 h2[i][j] = (h2[i - 1][j] * b + h[i][j]) % MOD 35 if i >= nn: 36 h2[i][j] = (h2[i][j] + h[i - nn][j] * (MOD - b2)) % MOD 37 return h, h2 38 39 def equal(a, p, ai, aj): 40 np = len(p) 41 mp = len(p[0]) 42 for i in xrange(np): 43 for j in xrange(mp): 44 if a[ai + i][aj + j] != p[i][j]: 45 return False 46 return True 47 48 def solve(): 49 na, ma = map(int, re.split('\s+', raw_input().strip())) 50 a = [] 51 for i in xrange(na): 52 a.append(map(int, list(raw_input().strip()))) 53 np, mp = map(int, re.split('\s+', raw_input().strip())) 54 p = [] 55 for i in xrange(np): 56 p.append(map(int, list(raw_input().strip()))) 57 ha, h2a = calcHash(a, np, mp) 58 hp, h2p = calcHash(p, np, mp) 59 60 for i in xrange(np - 1, na): 61 for j in xrange(mp - 1, ma): 62 if h2a[i][j] != h2p[np - 1][mp - 1]: 63 continue 64 if equal(a, p, i - np + 1, j - mp + 1): 65 print('YES') 66 return 67 print('NO') 68 69 if __name__ == '__main__': 70 t = int(raw_input()) 71 for ti in xrange(t): 72 solve() 73

?

轉(zhuǎn)載于:https://www.cnblogs.com/zhuli19901106/p/4687979.html

總結(jié)

以上是生活随笔為你收集整理的Hackerrank - The Grid Search的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。