HDU 5748 Bellovin（LIS）

Description
给出一个长度为n的序列a[i] (1<=i<=n)，定义f[i]为以a[i]结尾的最长上升子序列长度（严格上升），令F(a[1],…,a[n])=(f[1],…,f[n])，求字典序最小的序列b[i]，使得F(a[1],…,a[n])=F(b[1],…,b[n])，其中b[i]>=1(1<=i<=n)
Input
第一行一个整数T表示用例组数，每组用例首先输入一个整数n表示序列长度，之后n个正整数a[i]
Output
对于每组用例，输出满足条件的b[i]中字典序最小的那个
Sample Input
3
1
10
5
5 4 3 2 1
3
1 3 5
Sample Output
1
1 1 1 1 1
1 2 3
Solution
稍微修改一下LIS的模版即可得到f[i]，现在考虑字典序最小的b序列，由于b[i]为正整数，根据f[i]的定义可知对任意1<=i<=n有b[i]>=f[i]，而f序列本身就是一个满足条件的序列，故字典序最小的b序列即为f序列
Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
#define INF 0x3f3f3f3f
#define maxn 111111
int T,dp[maxn],n,a[maxn],f[maxn],b[maxn];
int LIS()
{
    for(int i=1;i<n;i++)dp[i]=INF;
    dp[0]=a[0];
    f[0]=1;
    int len=1;
    for(int i=1;i<n;i++)
    {
        if(a[i]>dp[len-1])
            f[i]=len+1,dp[len++]=a[i];
        else
        {
            int temp=lower_bound(dp,dp+n,a[i])-dp;
            dp[temp]=a[i]; 
            f[i]=temp+1;
        }
    }
    return len;
}
int main()
{
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        for(int i=0;i<n;i++)scanf("%d",&a[i]);
        LIS();
        for(int i=0;i<n;i++)
            printf("%d%c",f[i],i==n-1?'\n':' ');
    }
    return 0;
}