Home On The Range_usaco3.3_dp

题目描述 Description

---

农民约翰在一片边长是N (2 <= N <= 250)英里的正方形牧场上放牧他的奶牛。(因为一些原因，他的奶牛只在正方形的牧场上吃草。)遗憾的是,他的奶牛已经毁坏一些土地。( 一些1平方英里的正方形)

农民约翰需要统计那些可以放牧奶牛的正方形牧场(至少是2x2的,在这些较大的正方形中没有一个点是被破坏的，也就是说，所有的点都是“1”)。

你的工作要在被供应的数据组里面统计所有不同的正方形放牧区域(>=2x2)的个数。当然，放牧区域可能是重叠。

输入描述 Input Description

---

第 1 行:N,牧区的边长。

第 2 到　n+1行:N个没有空格分开的字符。

0 表示 “那一个区段被毁坏了”;1 表示 ” 准备好被吃”。

输出描述 Output Description

---

输出那些存在的正方形的边长和个数，一种一行。

题解 Analysis

---

中山做过类似的题哦，不过那是道暴力(笑
设$f[k][i][j]$表示以$i,j$为左上角长为 k <script type="math/tex" id="MathJax-Element-3">k</script>的正方形个数，转移显然
第一维可以滚，内存不会爆
然而这其实是二维dp我的方法比较容易懂吧

代码 Code

---

三围维

/*
ID:wjp13241
PROG:range
LANG:C++
*/
#include <stdio.h>
using namespace std;
int f[2][251][251],map[251][251],sum[251][251];
char s[251];
int main()
{
    freopen("range.in","r",stdin);
    freopen("range.out","w",stdout);
    int n;
    scanf("%d",&n);
    for (int i=1;i<=n;i++)
    {
        scanf("%s",s);
        for (int j=1;j<=n;j++)
        {
            map[i][j]=f[1][i][j]=s[j-1]-'0';
            sum[i][j]=sum[i-1][j]+sum[i][j-1]-sum[i-1][j-1]+map[i][j];
        }
    }
    int x=1;
    for (int k=2;k<=n;k++)
    {
        x^=1;
        int ans=0;
        for (int i=1;i<=n-k+1;i++)
            for (int j=1;j<=n-k+1;j++)
                if (f[x^1][i][j]&&sum[i+k-1][j+k-1]-sum[i-1][j+k-1]-sum[i+k-1][j-1]+sum[i-1][j-1]==k*k)
                {
                    f[x][i][j]++;
                    ans++;
                }
        if (ans)
            printf("%d %d\n",k,ans);
    }
    return 0;
}

二维

/*
ID:wjp13241
PROG:range
LANG:C++
*/
#include <stdio.h>
using namespace std;
int f[1001][1001],map[1001][1001];
char s[1001];
int min(int x,int y)
{
    return x<y?x:y;
}
int main()
{
    int n,m;
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++)
    {
        scanf("%s",s);
        for (int j=1;j<=m;j++)
            f[i][j]=map[i][j]=s[j-1]-'0';
    }
    int ans=0;
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
            if (map[i][j])
            {
                f[i][j]=min(min(f[i-1][j],f[i][j-1]),f[i-1][j-1])+1;
                ans+=f[i][j];
            }
    printf("%d\n",ans);
    return 0;
}