采油区域
Siruseri 政府决定将石油资源丰富的 Navalur 省的土地拍卖给私人承包商以
建立油井。被拍卖的整块土地为一个矩形区域,被划分为M×N 个小块。
Siruseri地质调查局有关于Navalur土地石油储量的估测数据。这些数据表示
为M×N 个正整数,即对每一小块土地石油储量的估计值。
为了避免出现垄断,政府规定每一个承包商只能承包一个由K×K 块相连的
土地构成的正方形区域。
AoE石油联合公司由三个承包商组成,他们想选择三块互不相交的 K×K 的
区域使得总的收益最大。
例如,假设石油储量的估计值如下:
1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1
1 8 8 8 8 8 1 1 1
1 8 8 8 8 8 1 1 1
1 8 8 8 8 8 1 1 1
1 1 1 1 8 8 8 1 1
1 1 1 1 1 1 8 8 8
1 1 1 1 1 1 9 9 9
1 1 1 1 1 1 9 9 9
如果 K = 2, AoE 公司可以承包的区域的石油储量总和为 100, 如果 K = 3,
AoE公司可以承包的区域的石油储量总和为 208。
AoE公司雇佣你来写一个程序,帮助计算出他们可以承包的区域的石油储量
之和的最大值。
输入格式
输入第一行包含三个整数 M, N, K,其中 M和N 是矩形区域的行数和列数,
K 是每一个承包商承包的正方形的大小(边长的块数) 。接下来 M行,每行有 N
个正整数表示这一行每一小块土地的石油储量的估计值。
输出格式
输出只包含一个正整数,表示 AoE 公司可以承包的区域的石油储量之和的
最大值。
采油区域
Siruseri 政府决定将石油资源丰富的 Navalur 省的土地拍卖给私人承包商以
建立油井。被拍卖的整块土地为一个矩形区域,被划分为M×N 个小块。
Siruseri地质调查局有关于Navalur土地石油储量的估测数据。这些数据表示
为M×N 个正整数,即对每一小块土地石油储量的估计值。
为了避免出现垄断,政府规定每一个承包商只能承包一个由K×K 块相连的
土地构成的正方形区域。
AoE石油联合公司由三个承包商组成,他们想选择三块互不相交的 K×K 的
区域使得总的收益最大。
例如,假设石油储量的估计值如下:
1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1
1 8 8 8 8 8 1 1 1
1 8 8 8 8 8 1 1 1
1 8 8 8 8 8 1 1 1
1 1 1 1 8 8 8 1 1
1 1 1 1 1 1 8 8 8
1 1 1 1 1 1 9 9 9
1 1 1 1 1 1 9 9 9
如果 K = 2, AoE 公司可以承包的区域的石油储量总和为 100, 如果 K = 3,
AoE公司可以承包的区域的石油储量总和为 208。
AoE公司雇佣你来写一个程序,帮助计算出他们可以承包的区域的石油储量
之和的最大值。
输入格式
输入第一行包含三个整数 M, N, K,其中 M和N 是矩形区域的行数和列数,
K 是每一个承包商承包的正方形的大小(边长的块数) 。接下来 M行,每行有 N
个正整数表示这一行每一小块土地的石油储量的估计值。
输出格式
输出只包含一个正整数,表示 AoE 公司可以承包的区域的石油储量之和的
最大值。
这种题到现在应该都是对于我们应该是很熟练的题目了。但是不知道为什么做不起。。。
代码比较长,但是确实是NOIP难度,一次AC了。
因为有三个块,因此可以分成六种情况
维护4种前缀和,然后递推算出rect1~4,表示前i行前j列最大的矩形,前i行后j列最大的矩形,后i行前j列最大的矩形,后i行后j列最大的矩形。
同时递推算出前i行右下角横坐标为j的最大矩形,前j列右下角纵坐标为i的最大矩形。
枚举两条割线就可以得出答案了。
一开始我想用悬线法用吃西瓜的那种方法来做,打完了调了很久了才发现有反例,结果错完了。总结出考试一定不能急躁,一定要考虑清楚再开动。
反例很简单就是这个,两个矩形的割线不一样
#include <cstdio>
long n;long m;long K;
long map[1510][1510];
long sum1[1510][1510];//left and above
long sum2[1510][1510];//right and above
long sum3[1510][1510];//left and below
long sum4[1510][1510];//right and below
long rect1[1510][1510];//left and above
long rect2[1510][1510];//right and above
long rect3[1510][1510];//left and below
long rect4[1510][1510];//right and below
long rectA[1510][1510];
long rectB[1510][1510];
#define MAX(a,b) (a>b?a:b)
int main()
{
freopen("oil.in","r",stdin);
freopen("oil.out","w",stdout);
scanf("%ld%ld%ld",&n,&m,&K);
for (long i=1;i<n+1;i++)
{
for (long j=1;j<m+1;j++)
{
scanf("%ld",map[i]+j);
sum1[i][j] = sum1[i][j-1]+map[i][j];
sum3[i][j] = sum1[i][j];
}
}
for (long i=1;i<n+1;i++)
{
for (long j=m;j>0;j--)
{
sum2[i][j] = sum2[i][j+1]+map[i][j];
sum4[i][j] = sum2[i][j];
}
}
for (long i=1;i<n+1;i++)
{
for (long j=1;j<m+1;j++)
{
sum1[i][j] += sum1[i-1][j];
sum2[i][j] += sum2[i-1][j];
}
}
for (long i=n;i>0;i--)
{
for (long j=m;j>0;j--)
{
sum3[i][j] += sum3[i+1][j];
sum4[i][j] += sum4[i+1][j];
}
}
for (long i=1;i<n+1;i++)
{
for (long j=1;j<m+1;j++)
{
if (i-K>=0 && j-K>=0)
{
rect1[i][j] = sum1[i][j]-sum1[i][j-K]-sum1[i-K][j]+sum1[i-K][j-K];//left and above
rectA[i][j] = rect1[i][j];
}
if (i-K>=0 && j+K-1<=m)
rect2[i][j] = sum2[i][j]-sum2[i][j+K]-sum2[i-K][j]+sum2[i-K][j+K];//right and above
if (i+K-1<=n && j-K>=0)
{
rect3[i][j] = sum3[i][j]-sum3[i][j-K]-sum3[i+K][j]+sum3[i+K][j-K];//left and blow
rectB[i][j] = rect3[i][j];
}
if (i+K-1<=n && j+K-1<=m)
rect4[i][j] = sum4[i][j]-sum4[i][j+K]-sum4[i+K][j]+sum4[i+K][j+K];//right and blow
}
}
for (long i=1;i<n+1;i++)
{
for (long j=1;j<m+1;j++)
{
rect1[i][j] = MAX(rect1[i][j],rect1[i-1][j]);
rect1[i][j] = MAX(rect1[i][j],rect1[i][j-1]);
rectA[i][j] = MAX(rectA[i][j],rectA[i][j-1]);
rectB[i][j] = MAX(rectB[i][j],rectB[i-1][j]);
}
}
for (long i=1;i<n+1;i++)
{
for (long j=m;j>0;j--)
{
rect2[i][j] = MAX(rect2[i][j],rect2[i-1][j]);
rect2[i][j] = MAX(rect2[i][j],rect2[i][j+1]);
}
}
for (long i=n;i>0;i--)
{
for (long j=1;j<m+1;j++)
{
rect3[i][j] = MAX(rect3[i][j],rect3[i+1][j]);
rect3[i][j] = MAX(rect3[i][j],rect3[i][j-1]);
}
}
for (long i=n;i>0;i--)
{
for (long j=m;j>0;j--)
{
rect4[i][j] = MAX(rect4[i][j],rect4[i+1][j]);
rect4[i][j] = MAX(rect4[i][j],rect4[i][j+1]);
}
}
long ans = 0;
for (long i=1+K;i<n-K+2;i++)
{
for (long j=1+K;j<m-K+2;j++)
{
ans = MAX(ans,rect1[i-1][j-1]+rect2[i-1][j]+rect4[i][1]);
ans = MAX(ans,rect3[i][j-1]+rect4[i][j]+rect2[i-1][1]);
ans = MAX(ans,rect1[i-1][j-1]+rect3[i][j-1]+rect4[1][j]);
ans = MAX(ans,rect2[i-1][j]+rect4[i][j]+rect3[1][j-1]);
}
}
for (long i=1+K;i<n-2*K+1;i++)
{
ans = MAX(ans,rect2[i-1][1]+rect4[i+K][1]+rectA[i+K-1][m]);
}
for (long j=1+K;j<m-2*K+1;j++)
{
ans = MAX(ans,rect3[1][j-1]+rect4[1][j+K]+rectB[n][j+K-1]);
}
printf("%ld",ans);
return 0;
}