Maximum Subarray

本文探讨了最大子序列和问题,并对比了直接求解法与分治算法两种方案。直接求解法虽简单易懂但效率较低,而分治算法尽管实现稍显复杂,但在处理大规模数据时更高效。
这周学了分治算法,挑了LeetCode 上的问题,问题内容如下:


Find the contiguous subarray within an array (containing at least one number) which has the largest sum.

For example, given the array [-2,1,-3,4,-1,2,1,-5,4],
the contiguous subarray [4,-1,2,1] has the largest sum = 6.


最大子序列是分治算法当中的经典问题了。但这个问题其实首先看到我是不会直接用分治算法去计算的,先考虑的还是直接把子项和依次相加作出比较,复杂度就是O(N^2),算法难度很小。可以直接通过for循环的嵌套完成,假设是在数组A中,A的大小为A_size。

int MaxSum = 0, sum; 
for(int i = 0; i < A_size; i++) {
    sum = 0;
    for(int j = i; j < A_size; j++) {
        sum += A[j];
        if(sum > MaxSum) {
            MaxSum = sum;
        }
    }
} 

也可以把解法优化为,首先判断首项是否小于0,因为在首项小于0的时候一定不会是最大子序列,所以第二个循环可以跳过小于0的项,会简单一些。


使用分治算法的原理也很容易理解,分治算法的基本思想就是化大为小,把一个大问题化为很多个相同解法的小问题,在这个求最大子序列的问题中,我们可以把整个序列分为左右两部分,左边部分从最右开始往左相加,右边部分从最左开始往右相加。可以猜想到,最大的数列会出现在左边、右边或者是中间的部分(即左边末尾加上右边开头)。而我们分成的每个部分里面也可以重复进行这个分割,直到只剩最后一个数字。

也就是说我们需要在这当中记录三个数字,一个是左边最大的子序列和,一个是右边最大的子序列和,一个是从左边末尾起相加最大的和加上右边开头起相加的最大的和,来进行比较,挑出最大的作为答案。

int Max(int maxleft, int maxright, int maxmiddle) {
    if(maxleft < maxright) {
    maxleft = maxright;
    }
    if(maxleft < maxmiddle) {
    maxleft = maxmiddle;
    }
    return maxleft;
}

int MaxArraySum(int A[], int left, int right) {
    //边界处,当递归进行到只有一个数字时,判断数字是否大于0,大于0则留下,小于0则不要
    if(left == right) {
        return A[left];
    } else {
        return 0;
    }

    //递归
    int center, maxLeft, maxRight;
    center = (left + right) / 2;
    maxLeft = MaxArraySum(A, left, center);
    maxRight = MaxArraySum(A, center + 1, right);
    //这里是对中间部分大小的计算
    int maxLeftBorderSum = 0, maxRightBorderSum = 0;
    int LeftBorderSum = 0, RightBorderSum = 0;
    int i;
    for (i = center; i >= left; i--) {
        LeftBorderSum += A[i];
        if(maxLeftBorderSum < LeftBorderSum) {
            maxLeftBorderSum = LeftBorderSum;
        }
    }
    for (i = center + 1; i <= right; i++) {
        RightBorderSum += A[i];
        if(maxRightBorderSum < RightBorderSum) {
            maxRightBorderSum = RightBorderSum;
        }
    }
    int MaxSum;
    MaxSum = Max(maxLeft, maxRight, maxLeftBorderSum + maxRightBorderSum);
    return MaxSum;
}


不过我还是觉得这个问题用分治算法有点麻烦,因为如果要输出最大子序列的话,考虑第一种方法会比较好输出,只需要记录下max出现时的起点A[i]和终点A[j]则能很快输出这个序列。


一、 内容概要 本资源提供了一个完整的“金属板材压弯成型”非线性仿真案例,基于ABAQUS/Explicit或Standard求解器完成。案例精确模拟了模具(凸模、凹模)与金属板材之间的接触、压合过程,直至板材发生塑性弯曲成型。 模型特点:包含完整的模具-工件装配体,定义了刚体约束、通用接触(或面面接触)及摩擦系数。 材料定义:金属板材采用弹塑性材料模型,定义了完整的屈服强度、塑性应变等真实应力-应变数据。 关键结果:提供了成型过程中的板材应力(Mises应力)、塑性应变(PE)、厚度变化​ 云图,以及模具受力(接触力)曲线,完整再现了压弯工艺的力学状态。 二、 适用人群 CAE工程师/工艺工程师:从事钣金冲压、模具设计、金属成型工艺分析与优化的专业人员。 高校师生:学习ABAQUS非线性分析、金属塑性成形理论,或从事相关课题研究的硕士/博士生。 结构设计工程师:需要评估钣金件可制造性(DFM)或预测成型回弹的设计人员。 三、 使用场景及目标 学习目标: 掌握在ABAQUS中设置金属塑性成形仿真的全流程,包括材料定义、复杂接触设置、边界条件与载荷步。 学习如何调试和分析大变形、非线性接触问题的收敛性技巧。 理解如何通过仿真预测成型缺陷(如减薄、破裂、回弹),并与理论或实验进行对比验证。 应用价值:本案例的建模方法与分析思路可直接应用于汽车覆盖件、电器外壳、结构件等钣金产品的冲压工艺开发与模具设计优化,减少试模成本。 四、 其他说明 资源包内包含参数化的INP文件、CAE模型文件、材料数据参考及一份简要的操作要点说明文档。INP文件便于用户直接修改关键参数(如压边力、摩擦系数、行程)进行自主研究。 建议使用ABAQUS 2022或更高版本打开。显式动力学分析(如用Explicit)对计算资源有一定要求。 本案例为教学与工程参考目的提供,用户可基于此框架进行拓展,应用于V型弯曲
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值