600字范文,内容丰富有趣,生活中的好帮手!
600字范文 > Lis算法如何用C语言实现? – 网络

Lis算法如何用C语言实现? – 网络

时间:2022-06-05 19:51:39

相关推荐

Lis算法如何用C语言实现? – 网络

问题描述:给出一个序列a1,a2,a3,a4,a5,a6,a7….an,求它的一个子序列(设为s1,s2,…sn),使得这个子序列满足这样的性质,s1<s2<s3<…<sn并且这个子序列的长度最长。输出这个最长的长度。(为了简化该类问题,大家将诸如最长下降子序列及最长不上升子序列等问题都看成同一个问题,其实仔细思考就会发现,这其实只是<符号定义上的问题,并不影响问题的实质)

例如有一个序列:1735948,它的最长上升子序列就是1348长度为4.

算法1(n^2):大家依次遍历整个序列,每一次求出从第一个数到当前这个数的最长上升子序列,直至遍历到最后一个数字为止,然后再取dp数组里最大的那个即为整个序列的最长上升子序列。大家用dp[i]来存放序列1-i的最长上升子序列的长度,那么dp[i]=max(dp[j])+1,(j∈[1,i-1]);显然dp[1]=1,大家从i=2开始遍历后面的元素即可。

下面是模板:

//最长上升子序列(n^2)模板

//入口参数:1.数组名称2.数组长度(注意从1号位置开始)

template<classT>

intLIS(Ta[],intn)

{

inti,j;

intans=1;

intm=0;

int*dp=newint[n+1];

dp[1]=1;

for(i=2;i<=n;i++)

{

m=0;

for(j=1;j<i;j++)

{

if(dp[j]>m&&a[j]<a[i])

m=dp[j];

}

dp[i]=m+1;

if(dp[i]>ans)

ans=dp[i];

}

returnans;

}

算法2(nlogn):维护一个一维数组c,并且这个数组是动态扩展的,初始大小为1,c[i]表示最长上升子序列长度是i的所有子串中末尾最小的那个数,根据这个数字,大家可以比较知道,只要当前考察的这个数比c[i]大,那么当前这个数一定能通过c[i]构成一个长度为i+1的上升子序列。当然大家希望在C数组中找一个尽量靠后的数字,这样大家得到的上升子串的长度最长,查找的时候使用二分搜索,这样时间复杂度便下降了。模板如下:

//最长上升子序列nlogn模板

//入口参数:数组名+数组长度,类型不限,结构体类型可以通过重载运算符实现

//数组下标从1号开始。

/**//////////////////////////BEGIN_TEMPLATE_BY_ABILITYTAO_ACM////////////////////////////

template<classT>

intbsearch(Tc[],intn,Ta)

{

intl=1,r=n;

while(l<=r)

{

intmid=(l+r)/2;

if(a>c[mid]&&a<=c[mid+1])returnmid+1;//>&&<=换为:>=&&<

elseif(a<c[mid])r=mid-1;

elsel=mid+1;

}

}

template<classT>

intLIS(Ta[],intn)

{

inti,j,size=1;

T*c=newT[n+1];

int*dp=newint[n+1];

c[1]=a[1];dp[1]=1;

for(i=2;i<=n;++i)

{

if(a[i]<=c[1])j=1;//<=换为:<

elseif(a[i]>c[size])

j=++size;//>换为:>=

else

j=bsearch(c,size,a[i]);

c[j]=a[i];dp[i]=j;

}

returnsize;

}

本内容不代表本网观点和政治立场,如有侵犯你的权益请联系我们处理。
网友评论
网友评论仅供其表达个人看法,并不表明网站立场。