USACO算法序列十七——range

    题目：http://www.nocow.cn/index.php/Translate:USACO/range

    寻找正方形内部所有的正方形块，包括任意小块。

    如果遍历每个节点，把每个节点当做正方形块的起始点，然后寻找正方形，会超时的。

    但是这道题可以通过动态规划的思想来考虑。map[i][j]示正方形右下角i,j节点的最大正方形的边长，则可以得出以下公式

    map[i][j]=min（map[i-1][j],map[i][j-1],map[i-1][j-1]）+ 1;

    因此就可以统计出所有的正方形块了。代码如下：

#include <iostream> #include <fstream> #define SIZE 250 using namespace std; int length; int map[SIZE][SIZE]; int result[SIZE+1] = {0}; int calRange(int i,int j) { int mini = map[i-1][j]; if (map[i-1][j] > map[i][j-1]) { mini = map[i][j-1]; } if (mini > map[i-1][j-1]) { mini = map[i-1][j-1]; } return mini + 1; } void calResult(int s) { for (int i=2; i <=s; i ++) { result[i] ++; } } int main() { ifstream fin("range.in"); ofstream fout("range.out"); fin >> length; //输入数据 char c; for (int i=0; i < length; i ++) { for (int j=0; j < length; j ++) { fin >> c; map[i][j] = c - '0'; } } //动态规划 //棋盘表示下节点为矩阵的方格数 for(int i=1; i<length; i ++) { for (int j=1; j<length; j ++) { if(map[i][j] == 1) { map[i][j] = calRange(i,j); if (map[i][j] > 1) { calResult(map[i][j]); } }//end if }//end for }//end for for (int i=2; i <=length; i ++) { if (result[i] != 0) { fout << i << " " << result[i] << endl; } } return 0; }

    运行效果如下：

Executing... Test 1: TEST OK [0.011 secs, 3256 KB] Test 2: TEST OK [0.000 secs, 3256 KB] Test 3: TEST OK [0.000 secs, 3256 KB] Test 4: TEST OK [0.000 secs, 3256 KB] Test 5: TEST OK [0.011 secs, 3256 KB] Test 6: TEST OK [0.000 secs, 3256 KB] Test 7: TEST OK [0.032 secs, 3256 KB] All tests OK.