求最大子序列的线性时间复杂度的解法的理解

本文探讨了《C++语言描述(第3版)_MarkAllenWeiss》中提到的最大子序列和问题,并详细解析了书中介绍的第四种求解方法。通过逆向思考和逐步分析,展示了该算法如何仅通过一次扫描就能找到最大子序列和,突出了其作为在线算法的独特优势。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

今天在看《C++语言描述(第3版)_Mark Allen Weiss》时,看到 第二章的四种求解最大子序列的方法,对于第四种解法,书上有言:“这是许多聪明算法的典型:运行时间是明显的,但正确性则不明显。对于。。。;然而即便如此许多仍然不信服。” 好吧,这句话激起了我的兴趣。

//最大子序列必然是序列中间的一部分:去掉序列的头部啊和尾部b,且知道,去掉的两部分必定小于0. 那么问题来了,中间这部分怎么找?
int maxSubSum4( const vector<int> & a )
{
    int maxSum = 0;
    int thisSum = 0;

    for(int j = 0; j <a.size(); j++ )
    {
        thisSum += a[j] ;

        if(thisSum > maxSum)  //maxSum只记录最大和,thisSum忠实记录每一和,只有大于maxSum才会被记录
            maxSum = thisSum;
        else if( thisSum < 0 )  //将a置为0,则不影响最大的和
            thisSum = 0;
    }
}

算法的好处:只对数据进行一次扫描,一旦 a[i] 被读入并被处理,它就不需要被记忆,更棒的是:在任意时候算法都能对它已经读入的数据给出子序列问题的正确答案,具有这种特性的算法叫做联机算法(on-line algorithm)

我对这种解法的解释:

我们不妨采取逆向理解:假设我们已经找到了一个数列的 最大 子序列 a[m]...........a[n]  (n >m) , 在此序列的前缀或是后缀,如果有不为负(也不为零)的子序列存在,则最大子序列的和将会增加,则与我们已经找到最大子序列矛盾,同时这里也向世人证明一条定理:最大子序列的前缀或是后缀子序列都不为正!




可以处理无限数据,线性复杂度 //by史瑞 //LFSR线性移位寄存器以及异或运算OK unsigned char LFSRB_M(unsigned char *fun,unsigned char *seq,unsigned long Cont){ unsigned long x; unsigned char ch=0x00,t=0x00,*array; array=(unsigned char *)malloc(Cont*sizeof(unsigned char)); for(x=0;x<Cont;x++) *(array+x)=(*(seq+x))&(*(fun+x)); for(x=0;x<Cont*8;x++){ ch=(*array)&0x80; LeftShift(array,Cont); t^=ch; } free(array); return t; } #define word(ln) ((ln-1)/8+1) #define place(ln) ((ln-1)%8) //Berlek_Massey生成任意序列的联接多项式OK unsigned char *Berlek_Massey(unsigned char *seq,unsigned long *Rank,unsigned long Cont,FILE *fmm){ unsigned long n=0,m,lm,ln=0; int d; unsigned long lfm; unsigned char *fun,*fm,*array,ch,t=0x00; unsigned char *func,*fmc; unsigned long x,y,k; fun=(unsigned char *)malloc(sizeof(unsigned char)); *(fun)=0x00; for(x=0;x<Cont*8;x++){ if(ln!=0){ array=(unsigned char *)malloc(word(ln)); InitialDSR(array,word(ln)); for(y=0;y<ln;y++){ ch=((*(seq+(n-(y+1))/8))<<((n-(y+1))%8))&0x80; if(ch) *(array+y/8)^=(ch>>(y%8)); } t=LFSRB_M(fun,array,word(ln)); d=((((*(seq+n/8))<<(n%8))&0x80;)^t)?1:0; free(array); } else d=(((*(seq+n/8))<<(n%8))&0x80;)?1:0; if(d){ if(ln!=0){ lm=ln; func=(unsigned char *)malloc(word(ln)*sizeof(unsigned char)); memcpy(func,fun,word(ln)); if(ln<(n+1-ln)){ ln=n+1-ln; } fmc=(unsigned char *)malloc(word(ln)*sizeof(unsigned char)); InitialDSR(fmc,word(ln)); memcpy(fmc,fm,word(lfm)); for(k=0;k<n-m;k++){ RightShift(fmc,word(ln)); } ch=0x80; *(fmc+(n-m-1)/8)^=(ch>>((n-m-1)%8)); fun=(unsigned char *)realloc(fun,word(ln)*sizeof(unsigned char)); for(k=word(lm);k<word(ln);k++) *(fun+k)=0x00; XorDSR(fun,fmc,word(ln)); free(fmc); if(lm<(n+1-lm)){ m=n; lfm=lm; fm=(unsigned char *)realloc(fm,word(lm)); memcpy(fm,func,word(lfm)); } free(func); } else{ ln=n+1; m=n; if(m!=0){ fm=(unsigned char *)malloc(word(m)*sizeof(unsigned char)); InitialDSR(fm,word(m)); lfm=m; } else{ fm=(unsigned char *)malloc(sizeof(unsigned char)); *fm=0x00; lfm=1; } fun=(unsigned char *)realloc(fun,word(ln)); InitialDSR(fun,word(ln)); ch=0x80; ch>>=(place(ln)); *(fun+(word(ln)-1))^=ch; } } n++; printf("\t<%d,%d>",n,ln); fprintf(fmm,"\t<n,d,ln>=<%d,%d,%d>",n,d,ln); if(n%3==0){ printf("\n"); fprintf(fmm,"\n"); } } printf("\nFn="); // for(k=0;k<word(ln);k++) // FromBytetoBit(*(fun+k)); *Rank=ln; return fun; }
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值