P1578 奶牛浴场
1.P1077 [NOIP2012 普及组] 摆花2.P3842 [TJOI2007] 线段3.P1164 小A点菜4.P2392 kkksc03考前临时抱佛脚5.P1102 A-B 数对6.P1002 [NOIP2002 普及组] 过河卒7.P163 银行贷款8.P1182 数列分段 Section II9.[NOIP2015 提高组] 跳石头10.P1734 最大约数和11.P8772 [蓝桥杯 2022 省 A] 求和12.P1216 [USACO1.5] [IOI1994]数字三角形 Number Triangles13.P5602 小 E 与美食14.P3817 小A的糖果15.P1130 红牌16.P1196 [NOI2002] 银河英雄传说17.P1955 [NOI2015] 程序自动分析18.P1621 集合19.P5250 【深基17.例5】木材仓库20.P2184 家谱21.[USACO16DEC] Cities and States S22.P4653 [CEOI2017] Sure Bet
23.P1578 奶牛浴场
24.P1387 最大正方形25.小鸟的设备26.P1439 【模板】最长公共子序列27.P1091 [NOIP2004 提高组] 合唱队形28.P1020 [NOIP1999 提高组] 导弹拦截29.P2340 [USACO03FALL] Cow Exhibition G30.P2285 [HNOI2004] 打鼹鼠31.P4995 跳跳!32.P3817 小A的糖果33.P5019 [NOIP2018 提高组] 铺设道路34.P1969 [NOIP2013 提高组] 积木大赛35.P1094 [NOIP2007 普及组] 纪念品分组36.P4447 [AHOI2018初中组] 分组最大子子矩形问题
首先一些概念
1.有效子矩阵: 内部不包含任何障碍点,且边界与坐标轴平行的子矩阵
2.极大子矩阵:一个有效子矩阵,如果不包含它,且比它大的有效子矩阵,则为极大有效子矩阵
3.最大有效子矩阵:所有有效矩阵中最大面积的子矩阵
极大化思想
定理1. 有一个障碍点的矩形中的最大子矩阵一定是极大子矩阵
定理2.一个极大矩阵的四边一定不能向外扩展,即每条边碰到障碍点或者边界
算法1
由于极大子矩形的4条边一定覆盖障碍点(或者边界) ——— 因此我们可以枚举每一个障碍点
1.把矩形的四个边界点也加入到障碍点的集合当中
s[++k].x = 0, s[k].y = 0; //左上边界
s[++k].x = 0,s[k].y = m; //右上边界
s[++k].x = n, s[k].y = 0; //左下边界
s[++k].x = n, s[k].y = m; //右下边界
2.对横坐标进行排序
3.枚举极大子矩阵的左边界,从左到右扫描每一个障碍点
4.不断修改可行的上下边界,从而得到所有以这个定点为左边界的极大子矩阵
一开始上下边界取的是牛场(矩形)的上下边界
更新上下边界
1.当前点的纵坐标比定点大 ,即当前点位于定点的上方,修改上边界为当前点的纵坐标
2.当前点位于定点的下方,则修改下边界为当前带你的纵坐标
3.当前点和定点位于同一位置,则结束扫描
for(int i=1;i<=k;i++)
{
x1 = s[i].x;
y1 = 0;
y2 = m;
for(int j=i+1;j<=k;j++)
{
x2 = s[j].x;
if(s[i].y > s[j].y) y1 = max(y1,s[j].y);
else y2 = min(y2,s[j].y);
ans = max(ans,(x2-x1) * (y2-y1));
}
}
5.枚举极大子矩阵的右边界,从右到左扫描每一个障碍点
由于从左到右扫描后可能会出现遗漏的情况,如极大子矩形的左边界是牛场的左边界,有边界覆盖一个障碍点
6.更新答案
ans = max(ans, (上边界 - 下边界 )* 横坐标之差)
7.还有,极大子矩阵的左右边界与矩形的左右边界重合的情况;
我们将按照纵坐标进行排序,枚举相邻的两个纵坐标之间的极大子矩形
for(int i=1;i<=k-1;i++){
ans = max(ans,n*(s[i+1].y-s[i].y));
}
C++ Code
#include <bits/stdc++.h>
using namespace std;
const int N = 3e5+10;
int n,m,k;
int ans;
struct Point {
int x,y;
}s[N];
bool cmp1(Point a, Point b){
if(a.x != b.x) return a.x < b.x;
return a.y < b.y;
}
bool cmp2(Point a, Point b){
if(a.y != b.y) return a.y < b.y;
return a.x < b.y;
}
int main(){
cin>>n>>m>>k;
for(int i=1;i<=k;i++){
int x,y;
cin>>x>>y;
s[i] = {x,y};
}
//加入矩形的四个边界点
s[++k].x = 0, s[k].y = 0;
s[++k].x = 0, s[k].y = m;
s[++k].x = n, s[k].y = 0;
s[++k].x = n, s[k].y = m;
int x1,x2,y1,y2;
sort(s+1,s+k+1,cmp1);
for(int i=1;i<=k;i++){
x1 = s[i].x;
y1 = 0, y2 = m;
for(int j=i+1;j<=k;j++){
x2 = s[j].x;
ans = max(ans,(x2-x1)*(y2-y1));
if(s[i].y > s[j].y) y1 = max(y1,s[j].y);
else y2 = min(y2,s[j].y);
}
}
for(int i=k;i>=1;i--){
x1 = s[i].x;
y1=0,y2=m;
for(int j=i-1;j>=1;j--){
x2 = s[j].x;
ans = max(ans,(x1-x2)*(y2-y1)); //右边减去左边
if(s[i].y > s[j].y) y1 = max(y1,s[j].y);
else y2 = min(y2,s[j].y);
}
}
sort(s+1,s+1+k,cmp2);
for(int i=1;i<=k-1;i++){
ans = max(ans,n*(s[i+1].y-s[i].y));
}
cout<<ans;
return 0;
}
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型,支持深度思考和联网搜索!
· 基于 Docker 搭建 FRP 内网穿透开源项目(很简单哒)
· ollama系列01:轻松3步本地部署deepseek,普通电脑可用
· 25岁的心里话
· 按钮权限的设计及实现