洛谷P3625 - [APIO2009]采油区域
2024-09-01 20:49:05
Description
给出一个\(n\times m(n,m\leq1500)\)的矩阵,从中选出\(3\)个互不相交的\(k\times k\)方阵,使得被选出的数的和最大。
Solution
奇怪做法...
三个矩形分别在三个部分中,把矩形划分成三部分只有这六种。首先搞出\(s[i][j]\)表示以\((i,j)\)为右下角的\(k\times k\)方阵的和,然后搞出\(f_1[i][j]\)表示\((1,1)-(i,j)\)中\(s\)的最大值,\(f_2[i][j]\)表示\((1,m)-(i,j)\)中\(s\)的最大值,\(f_3[i][j]\)表示\((n,m)-(i,j)\)中\(s\)的最大值,\(f_4[i][j]\)表示\((n,1)-(i,j)\)中\(s\)的最大值。枚举横竖划分在哪就可以解决四种。
平行的那两种搞出行/列最大值然后瞎搞即可。
时间复杂度\(O(nm)\)。
Code
//[APIO2009]Oil
#include <cstdio>
const int N=2000;
inline int max(int x,int y) {return x>y?x:y;}
int n,m,k,a[N][N];
int pre[N][N],s[N][N],f1[N][N],f2[N][N],f3[N][N],f4[N][N],row[N],col[N];
int main()
{
scanf("%d%d%d",&n,&m,&k);
int ans;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
scanf("%d",&a[i][j]);
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
pre[i][j]=pre[i-1][j]+pre[i][j-1]-pre[i-1][j-1]+a[i][j];
for(int i=k;i<=n;i++)
for(int j=k;j<=m;j++)
s[i][j]=pre[i][j]-pre[i-k][j]-pre[i][j-k]+pre[i-k][j-k];
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
f1[i][j]=max(s[i][j],max(f1[i-1][j],f1[i][j-1]));
for(int i=1;i<=n;i++)
for(int j=m;j>=1;j--)
f2[i][j]=max(s[i][j+k-1],max(f2[i-1][j],f2[i][j+1]));
for(int i=n;i>=1;i--)
for(int j=m;j>=1;j--)
f3[i][j]=max(s[i+k-1][j+k-1],max(f3[i+1][j],f3[i][j+1]));
for(int i=n;i>=1;i--)
for(int j=1;j<=m;j++)
f4[i][j]=max(s[i+k-1][j],max(f4[i+1][j],f4[i][j-1]));
for(int i=k;i<=n-k;i++)
for(int j=k;j<=m-k;j++)
{
ans=max(ans,f1[i][j]+f2[i][j+1]+f3[i+1][1]); //┴
ans=max(ans,f2[i][j+1]+f3[i+1][j+1]+f4[1][j]); //├
ans=max(ans,f3[i+1][j+1]+f4[i+1][j]+f1[i][m]); //┬
ans=max(ans,f4[i+1][j]+f1[i][j]+f2[n][j+1]); //┤
}
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
row[i]=max(row[i],s[i][j]),col[j]=max(col[j],s[i][j]);
for(int i=k;i<=n-k-k;i++)
for(int j=i+k,mid=row[j];j<=n-k;j++,mid=max(mid,row[j]))
ans=max(ans,f1[i][m]+mid+f3[j+1][1]);
for(int i=k;i<=n-k-k;i++)
for(int j=i+k,mid=col[j];j<=n-k;j++,mid=max(mid,col[j]))
ans=max(ans,f1[n][i]+mid+f3[1][j+1]);
printf("%d\n",ans);
return 0;
}
P.S.
写的我好难受...
最新文章
- 基于Ruby的Watir-WebDriver自动化测试框架
- XML Schema and XMLspy notes
- MyBatis学习(三)
- 【转载】css3 content 生成内容
- oracle的关闭过程(各个模式关闭)
- C# sogou地图API应用总结(二)
- docs/pcs/rest/file data apis list - 百度开发者中心
- Windows 7的 磁盘管理中,某个磁盘或分区,突然变成只读。
- ubuntu下使用openocd+jlink进行STM32开发调试
- Cassandra存储time series类型数据时的内部数据结构?
- ASP.NET没有魔法——ASP.NET Identity的加密与解密
- 推荐!手把手教你使用Git(转载)
- BOX
- angular $cookies、$cookieStore
- *args和**kwargs
- CentOS7 systemctl tomcat常用配置
- Intellij Idea15 快捷键设置大全
- HTML元素被定义为块级元素或内联元素。那么什么是块级元素,什么是内联元素呢
- C# Web开发中弹出对话框的函数[转载]
- 串口-CreateFile的使用