落谷P1020 导弹拦截（线性DP+树状数组+Dilworth）

最新推荐文章于 2025-04-07 22:23:53 发布

原创最新推荐文章于 2025-04-07 22:23:53 发布 · 592 阅读

1 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#线段树

动态规划同时被 2 个专栏收录

3 篇文章

订阅专栏

数据结构

3 篇文章

订阅专栏

本文探讨如何使用线性动态规划（DP）和树状数组优化解决导弹拦截问题，以降低时间复杂度。第一部分介绍如何利用树状数组在O(n log n)时间内找到最长不递增序列。第二部分则引入Dilworth定理，将问题转化为求最长上升序列的元素个数，从而解决导弹拦截的优化策略。

第一问，求最长不递增序列，普通DP复杂度为n*n，考虑用树状数组降低复杂度为n*log(n)

普通DP：

    for(int i=1;i<=n;i++) f[i]=1;//以i结尾的最长不上升序列长度
    for(int i=1;i<=n;i++){
        for(int j=i+1;j<=n;j++){
            if(a[j]<=a[i])
            {
            	f[j]=max(f[j],f[i]+1); 
            	ans1=max(ans1,f[j]);//更新ans1 
			}
                
        }
        
    }

树状数组：

将导弹高度作为数组下标，用树状数组不断更新区间最大长度，由于要用的是下标从1到下标为当前导弹高度的区间的最长不上升序列长度，因此可以不用线段树，用树状数组，查询时，从开头查到该点，更新时，从该点向后更新。

# include<iostream>
# include<stdio.h>
# include<algorithm>
# define maxn 100000+2 
# include<string.h>
using namespace std;
int bit[maxn];
int lowbit(int x)
{
	return x&(-x);
}

int query(int x)
{//向前求最大 
	int s=0;
	while(x)
	{
		s=max(s,bit[x]);//更新区间最大长度 
		x-=lowbit(x);
	}
	return s;
}


void add(int i,int x)
{//向后修改、添加 
	while(i<=maxn)
	{
		bit[i]=max(bit[i],x);
		i+=lowbit(i);
	}
}


int main()
{
	int ans,tp;
	int i=1;
	int a[maxn];
	while(scanf("%d",&a[i])!=EOF) i++;
	ans=0;
	for(int j=(i-1);j>=1;j--)
	{
		tp=query(a[j])+1;//最长不上升序列 
		add(a[j],tp);
		ans=max(ans,tp);
	}
	printf("%d\n",ans);
	ans=0;
	memset(bit,0,sizeof(bit));
	for(int j=1;j<=(i-1);j++)
	{
		tp=query(a[j]-1)+1;//最长上升序列 
		add(a[j],tp);
		ans=max(ans,tp);
	}
	printf("%d",ans);
}

第二问，用dilworth定理

在有穷偏序集中，任何反链最大元素数目等于任何将集合到链的划分中链的最小数目。一个关于无限偏序集的理论指出，在此种情况下，一个偏序集具有有限的宽度w，当且仅当它可以划分为最少w条链。

人话就是：一个序列中，不上升序列的集合的最少数目=最长上升序列的元素个数。

例如：题目样例：389 207 155 300 299 170 158 65

不上升序列集合最少划分为俩：{389 300 299 170 158 65}，{207 155}

最长上升序列为 155 300，元素个数为俩

因此第二问转化为求最长上升序列元素个数即可