分田地

题目描述

牛牛和 15 个朋友来玩打土豪分田地的游戏，牛牛决定让你来分田地，地主的田地可以看成是一个矩形，每个位置有一个价值。分割田地的方法是横竖各切三刀，分成 16 份，作为领导干部，牛牛总是会选择其中总价值最小的一份田地， 作为牛牛最好的朋友，你希望牛牛取得的田地的价值和尽可能大，你知道这个值最大可以是多少吗？

输入描述:

每个输入包含 1 个测试用例。每个测试用例的第一行包含两个整数 n 和 m（1 <= n, m <= 75），表示田地的大小，接下来的 n 行，每行包含 m 个 0-9 之间的数字，表示每块位置的价值。

输出描述:

输出一行表示牛牛所能取得的最大的价值。

示例1

输入

4 4
3332
3233
3332
2323

输出

2

---

/*
思路：
横竖各三刀将一个矩阵分为16部分，每部分的value即为这一部分包含的所有数字之和。
我们要做的就是想一种切法，使得这16部分中value最小的那个尽可能的大。
首先很显然，每一个部分的value在0-sum之间，sum是指整个矩阵所有数字之和，这样最终的结果一定是[0, sum]中的某一个整数。
这里稍微逆向思考一下，既然不容易直接求结果，可不可以我猜一个值k
然后判断能不能通过某种切法使最小的那一块value>=k,也就是说，使16块的value都能大于等于k,
如果可以的话，我们就可以对[0, sum]这个区间进行二分查找。
二分的复杂度是log(sum).
接下来是重点：对于一个值，怎么判断能不能横竖切三刀，使16块的value都大于等于k呢
二分答案，判断可行性时暴力枚举三列的情况，然后横着贪心地扫一遍,
每当四列都满足的时候就砍一刀，满足四次即可
*/

#include<iostream>
using namespace std;
int n,m;
int a[75][75];
int sum[76][76] = {0};
inline int calc(int x1,int y1,int x2,int y2)
{
    return sum[x2][y2] - sum[x2][y1] - sum[x1][y2] + sum[x1][y1];
}
bool getArea(int k)
{
    for(int x1=1;x1<=n-3;x1++)
    {
        for(int x2=x1+1;x2<=n-2;x2++)
        {
            for(int x3=x2+1;x3<=n-1;x3++)
            {
                int count = 0;
                int left_y = 0;
                for(int right_y=1;right_y<=m;right_y++)
                {
                    if(calc(0,left_y,x1,right_y)>=k
                      && calc(x1,left_y,x2,right_y)>=k
                      && calc(x2,left_y,x3,right_y)>=k
                      && calc(x3,left_y,n,right_y)>=k)
                    {
                        count++;
                        left_y = right_y;
                    }
                }
                if(count >= 4)
                    return true;
            }
        }
    }
    return false;
}
int main()
{
    while(cin >> n >> m)
    {
        string s;
        for(int i=0;i<n;i++)
        {
            cin >> s;
            for(int j=0;j<m;j++)
            {
                a[i][j] = s[j] -'0';
            }
        }
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=m;j++)
            {
                sum[i][j] = sum[i][j-1] + sum[i-1][j] - sum[i-1][j-1] + a[i-1][j-1];
            }
        }
        int l = 0,r = sum[n][m],mid,res = 0;
        while(l <= r)
        {
            mid = (l + r) / 2;
            if(getArea(mid))
            {
                l = mid + 1;
                res = mid;
            }
            else
                r = mid - 1;
        }
        cout << res << endl;
    }
    return 0;
}