题目描述
Byteasar想耕种他那块矩形的田,他每次能耕种矩形的一边(上下左右都行),在他每次耕完后,剩下的田也一定是矩形,每块小区域边长为11 ,耕地的长宽分别为mm 和nn ,不幸的是Byteasar只有一匹老弱的马,从马开始耕地开始,只有当它耕完了一边才会停下休息。但有些地会非常难耕以至于马会非常的累,因此Byteasar需要特别小心。当耕完了一边之后,马可以停下来休息恢复体力。每块地耕种的难度不一,但是Byteasar都非常清楚。我们将地分成m×n 块单位矩形——我们用坐标(i,j) 来定义它们。每块地都有一个整数t(i,j)
,来定义(i,j)的耕种难度。所以每次马耕一边地时的难度就是所有它耕种的地的难度总和,对于这匹虚弱的马而言,这个值不能超过他的体力值。Byteasar想知道在马不死掉的情况下最少需要耕多少次才能把地耕完。
题解
看了好久题似乎只能想到暴力和
O(n4)的DP
O
(
n
4
)
的
D
P
?
这种情况下很有可能正解是贪心的吧。
考虑最好情况下要耕多少次,显然是
min(m,n)
m
i
n
(
m
,
n
)
,但也显然由于体力限制是难以达到的。
但是这也给了我们一个启发,因为上面的方案肯定是只竖着耕或只横着耕,感性分析也不难发现,符合条件的最优情况一定是尽量只沿一个方向耕(由特殊到一般)
但我们不知道是竖着耕为主还是横着耕为主,于是可以分别贪心。
这里以竖着耕为主为例:
既然是竖着耕为主,那么就左耕耕,右耕耕。但是左右都耕不了了耕哪里呢?
如果
2n
2
n
枚举显然会TLE。于是又有一个经典的方法,可以先枚举一下在上面耕的最大次数(类似于迭代加深),这样最优解一定会被我们遍历到并且复杂度仅为
O((n+m)2)
O
(
(
n
+
m
)
2
)
代码如下:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<cmath>
#include<queue>
#define get_sum(a,b,c,d) (sum[c][d]-sum[a-1][d]-sum[c][b-1]+sum[a-1][b-1])
using namespace std;
inline int read()
{
int x=0;char ch=getchar();int t=1;
for(;ch>'9'||ch<'0';ch=getchar()) if(ch=='-') t=-1;
for(;ch>='0'&&ch<='9';ch=getchar()) x=(x<<1)+(x<<3)+(ch-48);
return x*t;
}
typedef long long ll;
const int N=2010;
ll sum[N][N];
ll K;
int n,m;
const int INF=1e9;
inline int Cut_l(int i,int j,int u,int v,int up,int ti)
{
if(i>u||j>v) return ti;
if(get_sum(i,j,u,j)<=K) return Cut_l(i,j+1,u,v,up,ti+1);
if(get_sum(i,v,u,v)<=K) return Cut_l(i,j,u,v-1,up,ti+1);
if(get_sum(i,j,i,v)<=K&&up) return Cut_l(i+1,j,u,v,up-1,ti+1);
if(get_sum(u,j,u,v)<=K) return Cut_l(i,j,u-1,v,up,ti+1);
return INF;
}
inline int Cut_r(int i,int j,int u,int v,int left,int ti)
{
if(i>u||j>v) return ti;
if(get_sum(i,j,i,v)<=K) return Cut_r(i+1,j,u,v,left,ti+1);
if(get_sum(u,j,u,v)<=K) return Cut_r(i,j,u-1,v,left,ti+1);
if(get_sum(i,j,u,j)<=K&&left) return Cut_r(i,j+1,u,v,left-1,ti+1);
if(get_sum(i,v,u,v)<=K) return Cut_r(i,j,u,v-1,left,ti+1);
return INF;
}
inline int qie_shu()
{
register int ans=INF;
for(register int i=0;i<=n;i++){
ans=min(ans,Cut_l(1,1,n,m,i,0));
}
return ans;
}
inline int qie_heng()
{
register int ans=INF;
for(register int i=0;i<=n;i++){
ans=min(ans,Cut_r(1,1,n,m,i,0));
}
return ans;
}
int main()
{
K=read();m=read();n=read();
for(register int i=1;i<=n;i++)
for(register int j=1;j<=m;j++)
sum[i][j]=1ll*read()+sum[i-1][j]-sum[i-1][j-1]+sum[i][j-1];
register int ans=qie_shu();
ans=min(ans,qie_heng());
printf("%d\n",ans);
}