分割矩阵 (二分范围[L,R)）

   分割矩阵

                   (browine.c/cpp/pas)

【问题描述】

    有N*M的一个非负整数矩阵。现在要把矩阵分成A*B块。矩阵先水平地切A-1刀，把矩阵划分成A块。然后再把剩下来的每一块独立地切竖着B-1刀。每块的价值为块上的数字和。求一种方案，使得最小价值块的价值最大。

【输入格式】

    第一行四个整数N,M,A,B。

    接下来N行，每行M个非负整数。代表这个矩阵

【输出格式】

    一个数字。最小价值块的价值。

【样例输入】

5 4 4 2

1 2 21

3 1 1 1

2 0 1 3

1 1 1 1

1 1 11

【样例输出】

    3

样例解释见图片

【数据规模】

   1<=A<=N<=500

   1<=B<=M<=500

    其他数字小于4000。

 二分的同志请注意了

m=(l+r)/2 <-----这句是永远也枚不到L与R的

#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<iostream>
#include<functional>
#include<algorithm>
using namespace std;
#define MAXN (500+10)
#define MAXM (500+10)
#define MAXT (2000000+10)
int n,m,t1,t2,a[MAXN][MAXM],sum[MAXN][MAXM]={0};
bool is_ok(int l,int r,int _m)
{
  int tot=0,p=0;
  for (int i=1;i<=m;i++)
  {
    tot+=sum[r][i]-sum[l-1][i];
    if (tot>=_m) {tot=0;p++;}
  }
  if (p>=t2) return 1;
  else return 0;
}
bool is_ok_(int _m)
{
  int p=0,l=1;
  for (int i=1;i<=n;i++)
  {
    if (is_ok(l,i,_m)) {l=i+1;p++;}
  }
  if (p>=t1) return 1;
  else return 0;
}
int main()
{
  freopen("browine.in","r",stdin);
  freopen("browine.out","w",stdout);
  scanf("%d%d%d%d",&n,&m,&t1,&t2);
  for (int i=1;i<=n;i++)
    for (int j=1;j<=m;j++)
    {
      scanf("%d",&a[i][j]);
      sum[i][j]=sum[i-1][j]+a[i][j];
    }
  /*
  for (int i=1;i<=n;i++)
    for (int j=1;j<=m;j++)
    {
      printf("%d ",sum[i][j]);
    }
  */
//  cout<<(is_ok_(4));
  int l=1,r=1,ans=0;
  for (int j=1;j<=m;j++) r+=sum[n][j];
  
  for (int i=1;i<=60;i++)
  {
    int m_=(l+r)/2;
    if (is_ok_(m_)) {l=ans=m_;}
    else r=m_;        
  }
  printf("%d\n",ans);
    
//  while (1);
  return 0;
}