农田个数【DP】

题目：

你的老家在河北农村。过年时，你回老家去拜年。你家有一片NM农田，将其看成一个NM的方格矩阵，有些方格是一片水域。你的农村伯伯听说你是学计算机的，给你出了一道题： 他问你：这片农田总共包含了多少个不存在水域的正方形农田。

两个正方形农田不同必须至少包含下面的两个条件中的一条：

边长不相等

左上角的方格不是同一方格

输入:

输入数据第一行为两个由空格分开的正整数N、M（1<=m 第2行到第N+1行每行有M个数字（0或1），描述了这一片农田。0表示这个方格为水域，否则为农田（注意：数字之间没有空格，而且每行不会出现空格）

样例输入:

3 3                                               
1 1 0
1 1 0
0 0 0     

输出:

满足条件的正方形农田个数。

样例输出:

5

思路：

(1)设F[I，J]表示以方格（I，J）为右下角，可以得到的最大无水正方形边长，那么显然如果（I，J）是水域，F[I，J]=0；否则F[I，J]=Min{F[I-1，J]，F[I，J-1]，F[I-1，J-1]}+1

(2)求出了F数组的值后，我们可以用F1[I]表示F数组中，值为I的个数。

(3)显然，假定边长为I的正方形数目为Sum[I]，那么有Sum[I]=Sum[I+1]+F1[I]

(4)最后只要算出Sum数组各个值的和为问题的解。

代码：

#include<cstdio>
#include<iostream>
using namespace std;
int n,m,sum,f[1001][1001];//设变量。
char c;
int main()
{
	scanf("%d%d",&n,&m);//输入。
	for(int i=1;i<=n;i++)
	for(int j=1;j<=m;j++)
	{
		c=getchar();//对空格的处理
		while(!(c>='0')&&(c<='9'))c=getchar();//转换。
		if (c==49) 
		f[i][j]=min(f[i-1][j-1],min(f[i-1][j],f[i][j-1]))+1;//动态转移方程
		sum=sum+f[i][j];//累加。
	}
	printf("%d",sum);//输出。
	return 0;
}