【JZOJ3918】蛋糕【二分】

题目：

题目链接：https://jzoj.net/senior/#main/show/3918
给出一个$n\times m$的矩阵，将这个矩阵横切三下，竖切三下，使得每一块的数字之和的最小值最大。

---

思路：

枚举横切分别切在哪里。然后可以二分最小值，转变为判定性问题。
每次二分就$O(n)$暴力向右扫，看看列能否分成$4$部分，每部分每一块都$\ge mid$。
时间复杂度$O(n^{3}m\log ({\sum }_{i=1}^n{\sum }_{j=1}^{m}a[i][j]))$。
跑不满常数小OJ快。能过

---

代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

const int N=80;
int n,m,ans,map[N][N],sum[N][N],cut[5];

int count(int x,int xx,int y,int yy)
{
	return sum[xx][yy]-sum[x-1][yy]-sum[xx][y-1]+sum[x-1][y-1];
}

bool check(int mid)
{
	for (register int j=1,i=1;j<=4;j++,i++)
	{
		int last=i;
		for (;;i++)
		{
			if (i>m) return 0;
			bool flag=1;
			for (register int k=1;k<=4;k++)
				if (count(cut[k-1]+1,cut[k],last,i)<mid)
				{
					flag=0;
					break;
				}
			if (flag) break;
		}
	}
	return 1;
}

int solve()
{
	int l=0,r=sum[n][m]/16+1,mid;
	while (l<=r)
	{
		mid=(l+r)>>1;
		if (check(mid)) l=mid+1;
			else r=mid-1;
	}
	return l-1;
}

int main()
{
	scanf("%d%d",&n,&m);
	for (int i=1;i<=n;i++)
		for (int j=1;j<=m;j++)
		{
			scanf("%1d",&map[i][j]);
			sum[i][j]=sum[i-1][j]+sum[i][j-1]-sum[i-1][j-1]+map[i][j];
		}
	cut[4]=n;
	for (cut[1]=1;cut[1]<n;cut[1]++)
		for (cut[2]=cut[1]+1;cut[2]<n;cut[2]++)
			for (cut[3]=cut[2]+1;cut[3]<n;cut[3]++)  //枚举每一行在哪里切割
				ans=max(ans,solve());
	printf("%d",ans);
	return 0;
}