URAL 1146

URAL 1146题目大意:给出一个N*N的矩阵,每个元素的范围为[-127,+127]。求其所有子矩阵中各元素之和的最大值。TimeLimit:1000MS     MemoryLimit:65536KB     64bitIOFormat:%I64d&%I64u数据规模:N理论基础:无。题目分析:求和问题,用前缀和(每一行)+枚举(对列的起始与终点和行的起始与终点枚举)即

大家好,又见面了,我是你们的朋友全栈君。

题目大意:给出一个N*N的矩阵,每个元素的范围为[-127,+127]。求其所有子矩阵中各元素之和的最大值。

Time Limit:1000MS     Memory Limit:65536KB     64bit IO Format:%I64d & %I64u

数据规模:N<=100。

理论基础:无。

题目分析:求和问题,用前缀和(每一行)+枚举(对列的起始与终点和行的起始与终点枚举)即可,复杂度O(n^3)。因为n<=100,所以不会超时。

需要注意的是:当列宽一定时。如果前面出现最大和ans=sum,那么只要sum加上下一行没有变成<=0,那么可以肯定从出现ans的地方直到当前位置不可能存在和比ans大的子矩阵。也可以证明,后面如有子矩阵的和比ans大,肯定是在sum被减为0以后,否则:sum+这个子矩阵的和比这个子矩阵的和大。综上,每当sum被减为0以后,便可以以此处为新起点重新搜索子矩阵。

代码如下:

#include<iostream>
#include<cstring>
#include<string>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<queue>
#include<ctime>
#include<vector>
using namespace std;
typedef double db;
#define DBG 1
#define maa (1<<31)
#define mii ((1<<31)-1)
#define ast(b) if(DBG && !(b)) { printf("%d!!|\n", __LINE__); while(1) getchar(); }  //调试
#define dout DBG && cout << __LINE__ << ">>| "
#define pr(x) #x"=" << (x) << " | "
#define mk(x) DBG && cout << __LINE__ << "**| "#x << endl
#define pra(arr, a, b)  if(DBG) {\
    dout<<#arr"[] |" <<endl; \
    for(int i=a,i_b=b;i<=i_b;i++) cout<<"["<<i<<"]="<<arr[i]<<" |"<<((i-(a)+1)%8?" ":"\n"); \
    if((b-a+1)%8) puts("");\
}
template<class T> inline bool updateMin(T& a, T b) { return a>b? a=b, true: false; }
template<class T> inline bool updateMax(T& a, T b) { return a<b? a=b, true: false; }
typedef long long LL;
typedef long unsigned int LU;
typedef long long unsigned int LLU;
#define N 100
int a[N+1][N+1],sum[N+1][N+1];
int n,ans,ma,anst;
int main()
{
    memset(a,0,sizeof a);
    memset(sum,0,sizeof sum);
    while(~scanf("%d",&n))
    {
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=n;j++)
            {
                scanf("%d",a[i]+j);
                sum[i][j]=sum[i][j-1]+a[i][j];
            }
        }
        ans=maa;
        for(int i=1;i<=n;i++)
        {
            for(int j=i;j<=n;j++)
            {
                ma=maa,anst=maa;
                for(int k=1;k<=n;k++)
                {
                    if(ma>0)
                    {
                        ma=ma+sum[k][j]-sum[k][i-1];
                    }
                    else ma=sum[k][j]-sum[k][i-1];
                    if(ma>anst)anst=ma;
                }
                if(anst>ans)ans=anst;
            }
        }
        printf("%d\n",ans);
    }
	return 0;
}

其中前缀和在输入矩阵时即可完成。

by:Jsun_moon http://blog.csdn.net/jsun_moon

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/132957.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • 贝云cms内容管理系统(thinkphp5.0开源cms管理系统)

    贝云cms内容管理系统(thinkphp5.0开源cms管理系统)

    2021年11月6日
    75
  • ITIL V3与ITIL V2的价值差异[通俗易懂]

    ITIL V3与ITIL V2的价值差异[通俗易懂]之前有篇日志写了ITILv3的介绍。这里说一下v3和v2才差异 ITILV3自从2007年推出后,已经将近两年了,这两年时间足够令相关的研究者和爱好者能知晓这个名词。对于新的名词、概念或者技术推出后,往往会有三个时期存在较多的争论,最热闹的时候一般是刚推出或临近推出的时期,而后慢慢冷淡一段时间后,不少有机会深入这些新事物的人慢慢地又会有第二个阶段的争论,最后一个阶段常常是这个新事物真正

    2022年10月6日
    0
  • mysql事务回滚机制概述「建议收藏」

    mysql事务回滚机制概述

    2022年2月17日
    31
  • 无刷直流电机simulink模型(ansoft电机仿真教程)

    这段时间刚开始接触Matlab中的Simulink仿真,我就结合自己的专业,利用Simulink进行了无刷直流电机的仿真,因为Simulink工具箱里面有很多可用的模块,所以建模过程变得非常简单。在Matlab界面中new->model之后,找到SimulinkLibraryBrowser,这里面有系统自带的很多模块,接下来我们就需要找到自己所需要的模块了,了解无刷直流电机的工作原理以后…

    2022年4月11日
    57
  • ie浏览器兼容性视图灰色怎么办_浏览器没有兼容性视图设置

    ie浏览器兼容性视图灰色怎么办_浏览器没有兼容性视图设置2018-07-1010:57:17win8系统浏览一些高级网站会出现打不开的问题,其实我们需要设置兼容性视图才能正常浏览,那么怎么设置IE8浏览器兼容性视图?其实设置IE8浏览器兼容性视图的方法很简单,接下去教程和…2015-06-0516:15:42Win8系统自带的IE11浏览器有兼容模式功能,通过开启兼容模式,可以让用户在将浏览器关闭后再重启,然后在测试下输入网站即可正常打开一些兼…

    2025年9月2日
    5
  • Tracert(traceroute)&Ping 工作原理分析

    Tracert(traceroute)&Ping 工作原理分析一、tracert工作过程分析Tracert命令用IP生存时间(TTL)字段和ICMP错误消息来确定从一个主机到网络上其他主机的路由。首先,tracert送出一个TTL是1的IP数据包到目的地,当路径上的第一个路由器收到这个数据包时,它将TTL减1。此时,TTL变为0,所以该路由器会将此数据包丢掉,并送回一个「ICMPtimeexceeded」消息(

    2022年9月24日
    1

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号