First of all, if any of the edges of the National Park is not axis-parallel, then the National Park cannot be a union of axis-parallel rectangles. Furthermore, every rectangular plot which does not lie entirely in the National Park can be safely discarded. To compute whether a rectangle plot lies entirely inside the Nationl Park, it suffices to check if none of its sides intersects any edge of the National Park, and check if some interior point of the plot (say the center) lies inside the Park. What remains is a set of plots which lie completely inside the National Park. Since the plots are mutually disjoint, the plots cover the National Park if and only if the area of the National Park and the sum of the areas of the plots are equal. The time complexity of this algorithm is O(MN).
To compute the area of a polygon (x1,y1), (x2,y2), ..., (xN, yN) we use the known equation area = ∑ (xi-xi+1)*(yi+1+yi)/2, where we sum for i=1,2,...,N and the (N+1)-st vertex equals to the first one.
The easiest way to determine whether a point A lies inside the polygon (x1,y1), (x2,y2), ..., (xN, yN) is to find the sum of oriented angles ∑ (xi,yi)A(xi+1, yi+1). If A lies inside the polygon the sum is ± 2Π. (The sign depends on the orientation of the polygon.) Otherwise, the sum is zero.
All other computations are trivial.