ZOJ 4028 LIS(贪心)

本文介绍了一种构造最长严格上升子序列的算法,通过预处理和贪心策略确保序列符合给定条件。讨论了如何根据序列元素的值确定其在序列中的位置,并给出具体实现代码。

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

Description

有一个长度为nn的序列a1,...,an,以fifi表示以aiai结尾的最长严格上升子序列长度,先给出fifiaiai的取值范围[li,ri][li,ri],给出满足条件的序列a1,...,ana1,...,an,保证有解

Input

第一行一整数TT表示用例组数,每组用例首先输入一整数n表示序列长度,之后输入nn个整数f1,...,fn,最后nn行每行输入两个整数li,ri表示aiai的取值范围

(1n105,1fin,0liri2109)(1≤n≤105,1≤fi≤n,0≤li≤ri≤2⋅109)

Output

输出满足条件的aa序列

Sample Input

4
6
1 2 3 2 4 3
0 5
2 4
3 3
1 2
3 5
1 5
5
1 2 1 3 1
100 200
200 300
200 400
400 500
100 500
7
1 2 3 1 1 4 2
0 3
0 3
0 3
0 3
0 3
0 3
0 3
2
1 1
1 2
2 3

Sample Output

1 2 3 2 5 3
200 300 200 500 200
0 1 2 0 0 3 1
2 2

Solution

fi=1,说明a1,...,ai1a1,...,ai−1均不小于aiai,那么所有fi=1fi=1的位置的取值应该是一个不增的序列,贪心的让最后一个fi=1fi=1的位置取最小值lili,之后从后往前把所有fi=1fi=1的位置都取最小值(两个限制:所给取值限制以及保证所有fi=1fi=1的位置取值序列不增的限制),这样一定对其他位置取值影响最小且合法

fi=x>1fi=x>1,说明存在1j<i1≤j<i,使得fj=x1fj=x−1aj<aiaj<ai,也即第ii个位置的取值要大于前i个位置中fj=x1fj=x−1的位置取值最小值,由于所有fj=x1fj=x−1的位置取值不增,故第ii个位置的取值只要大于它前面第一个fj=x1的位置的取值即可,同样的,所有fi=xfi=x的位置取值不增,从后往前让每个位置取最小值即可(三个限制:所给取值限制,保证所有fi=xfi=x的位置取值序列不增的限制,大于该位置前面第一个fj=x1fj=x−1的位置取值的限制)

故对于每个位置ii,预处理pre[i]表示第ii个位置前第一个满足fj=fi1的位置jj,把所有位置按f取值升序排列,同样ff取值的位置按下标降序排列,以此序一个个按之前分析确定取值即可

Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
typedef pair<int,int>P;
const int INF=0x3f3f3f3f,maxn=100005;
struct node
{
    int l,r,f,id;
    bool operator<(const node&b)const
    {
        if(f!=b.f)return f<b.f;
        return id>b.id;
    }
}a[maxn];
int T,n,ans[maxn],Next[maxn],pos[maxn];
int main()
{
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        for(int i=1;i<=n;i++)scanf("%d",&a[i].f),a[i].id=i;
        for(int i=1;i<=n;i++)scanf("%d%d",&a[i].l,&a[i].r);
        memset(Next,0,sizeof(Next));
        for(int i=1;i<=n;i++)
        {
            pos[i]=Next[a[i].f-1];
            Next[a[i].f]=i;
        }
        sort(a+1,a+n+1);
        int L,F;
        for(int i=1;i<=n;i++)
        {
            if(pos[a[i].id])a[i].l=max(a[i].l,ans[pos[a[i].id]]+1); 
            if(i==1||a[i].f!=a[i-1].f)
            {
                F=a[i].f,L=a[i].l;
                ans[a[i].id]=L;
                continue;
            }
            ans[a[i].id]=L=max(a[i].l,L);
        }
        for(int i=1;i<=n;i++)printf("%d%c",ans[i],i==n?'\n':' ');
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值