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

来源:互联网 发布:matlab的分块矩阵运算 编辑:程序博客网 时间:2024/05/22 16:57

   分割矩阵

                   (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;}



 

原创粉丝点击