poj 1054 讨厌的青蛙

这个问题看起来很复杂，其实目的很简单：帮助农民找到为害最大的青蛙。也就是要找到
一条穿越稻田的青蛙路径，这个路径上被踩踏的水稻不少于其他任何青蛙路径上被踩踏的水
稻数。当然，整个稻田中也可能根本就不存在青蛙路径。问题的关键是：找到穿越稻田的全
部青蛙路径。任何一条穿越稻田的青蛙路径L，至少包括3 棵被踩踏的水稻。假设其中前两
棵被踩踏的水稻分别是(X1，Y1)、(X2，Y2)，那么：
 令dx=X2-X1、dy=Y2-Y1；X0=X1-dx、Y0=Y1- dy；X3=X2 + dx、Y3=Y2 + dy
 (X0，Y0)位于稻田之外，青蛙从该位置经一跳后进入稻田、踩踏位置(X1，Y1)上的水稻
 (X3，Y3)位于稻田之内，该位置是L 上第3 棵被青蛙踩踏的水稻
 Xi=X0 + idx、Yi=Y1 + idy(i3)，如果(Xi，Yi)位于稻田之内，则(Xi，Yi)上的水稻必被
青蛙踩踏
根据上述规则，只要知道一条青蛙路径上的前两棵被踩踏的水稻，就可以找到该路径上其
他的水稻。为了找到全部的青蛙路径，只要从被踩踏的水稻中，任取两棵水稻(X1，Y1)、(X2，
Y2)，判断(X1，Y1)、(X2，Y2)是否能够作为一条青蛙路径上最先被踩踏的两颗水稻。
解决方案
这个问题的描述中，最基本的元素是被踩踏的水稻。在程序中要选择一个合适的数据结构，
来表达这个基本元素。这个数据结构是否合适的标准是：在程序中要表达这个元素时，能否
用一个单词或者短语，即用一个变量来表示。
struct PLANT {//描述一棵被踩踏的水稻
int x; //水稻的行号
int y; //水稻的列号
}
这个问题的主要计算是：从被踩踏的水稻中选择两棵(X1，Y1)、(X2，Y2)。判断它们是否
能够作为一条青蛙路径上最先被踩踏的两颗水稻。(X1，Y1)、(X2，Y2)唯一确定了蛙跳的方
向和步长，从(X2，Y2)开始，沿着这个方向和步长在稻田内走。每走一步，判断所到达位置
上(X，Y)的水稻是否被踩踏，直到走出稻田为止。如果在某一步上，(X，Y)没有被踩踏，
则表明(X1，Y1)、(X2，Y2)是一条青蛙路径上最先被踩踏的两颗水稻的假设不成立。这个判
断的算法在问题求解过程中要反复使用，它的效率成为决定整个计算效率的关键。
 用一个PLANT 型的数组plants[5001]表示全部被踩踏的水稻
 将plants 中的元素按照行/列序号的升序(或者降序)排列
 采用二分法查找plants 中是否有值为(X，Y)的元素：将(X，Y)与plants 中间的元素比较，
(1)相等，表明找到了元素；(2)比plants 中间元素的小，继续在plants 的前半部寻找；(3)
比plants 中间元素的大，继续在plants 的后半部寻找。
采用上述方法判断每走一步所到达位置上(X，Y)的水稻是否被踩踏，最多只要比较Log2N，
其中N 是稻田中被踩踏水稻的总量。

#include <iostream> #include <cstdio> #include <cstdlib> using namespace std; int r, c, n; struct PLANT { int x, y; }; PLANT plants[5001]; PLANT plant; int myCompare( const void *ele1, const void *ele2 ); int searchPath(PLANT secPlant, int dX, int dY); int main( ) { int i, j, dX, dY, pX, pY, steps, max = 2; scanf("%d%d", &r, &c); scanf("%d", &n); for ( i = 0; i < n; i++ ) scanf("%d %d", &plants[i].x, &plants[i].y); qsort(plants, n, sizeof(PLANT), myCompare); for ( i = 0; i < n - 2; i++ ) for ( j = i + 1; j < n - 1; j++ ) { dX = plants[j].x - plants[i].x; dY = plants[j].y - plants[i].y; pX = plants[i].x - dX; pY = plants[i].y - dY; if ( pX <= r && pX >= 1 && pY <= c && pY >= 1 ) continue; if ( plants[i].x + max * dX > r ) break; pY = plants[i].y + max * dY; if ( pY > c || pY < 1) continue; steps = searchPath( plants[j], dX, dY ); if ( steps > max ) max = steps; } if ( max == 2 ) max = 0; printf("%d/n", max); return 0; } int myCompare( const void *ele1, const void *ele2 ) { PLANT *p1, *p2; p1 = (PLANT*) ele1; p2 = (PLANT*) ele2; if ( p1->x == p2->x ) return( p1->y - p2->y ); return ( p1->x - p2->x ); } int searchPath( PLANT secPlant, int dX, int dY ){ PLANT plant; int steps; plant.x = secPlant.x + dX; plant.y = secPlant.y + dY; steps = 2; while ( plant.x <= r && plant.x >= 1 && plant.y <= c && plant.y >= 1 ) { if ( !bsearch(&plant, plants, n, sizeof(PLANT), myCompare) ) { steps = 0; break; } plant.x += dX; plant.y += dY; steps++; } return( steps ); }