【LuoguP3444】[POI2006]ORK-Ploughing

题目链接

题目描述

Byteasar想耕种他那块矩形的田，他每次能耕种矩形的一边（上下左右都行），在他每次耕完后，剩下的田也一定是矩形，每块小区域边长为11 ，耕地的长宽分别为mm 和nn ，不幸的是Byteasar只有一匹老弱的马，从马开始耕地开始，只有当它耕完了一边才会停下休息。但有些地会非常难耕以至于马会非常的累，因此Byteasar需要特别小心。当耕完了一边之后，马可以停下来休息恢复体力。每块地耕种的难度不一，但是Byteasar都非常清楚。我们将地分成m×n 块单位矩形——我们用坐标(i,j) 来定义它们。每块地都有一个整数t(i,j)
​ ,来定义(i,j)的耕种难度。所以每次马耕一边地时的难度就是所有它耕种的地的难度总和，对于这匹虚弱的马而言，这个值不能超过他的体力值。Byteasar想知道在马不死掉的情况下最少需要耕多少次才能把地耕完。

题解

看了好久题似乎只能想到暴力和$O(n^4)的DP$?
这种情况下很有可能正解是贪心的吧。
考虑最好情况下要耕多少次，显然是$min(m,n)$，但也显然由于体力限制是难以达到的。
但是这也给了我们一个启发，因为上面的方案肯定是只竖着耕或只横着耕，感性分析也不难发现，符合条件的最优情况一定是尽量只沿一个方向耕(由特殊到一般)
但我们不知道是竖着耕为主还是横着耕为主，于是可以分别贪心。

这里以竖着耕为主为例:
既然是竖着耕为主，那么就左耕耕，右耕耕。但是左右都耕不了了耕哪里呢?
如果$2^n$枚举显然会TLE。于是又有一个经典的方法，可以先枚举一下在上面耕的最大次数(类似于迭代加深),这样最优解一定会被我们遍历到并且复杂度仅为$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);
}