XJOI NOIP2016提高组冲剌题1 T2:道路规划(LIS)

XJOI NOIP2016提高组冲剌题1 T2:道路规划

题目分析:

考虑这样一个问题,对于有n个城市的已知集合,如何判断其合法性?
那么必须保证最左边的对应最右边的,那尝试把下方的序列反转一下,
则保障左边的对应左边——即求最长公共子序列,其转移方程如下,

显然时间复杂度为O(n^2),但是n<=10^5,显然是会超时的,考虑如何优化呢?
既然对于上面某一确定A序列,要使其合法,那么B序列必须保证反转后与原序列顺序一致.则设rk[i]表示i编号的城市在上方序列对应的排序号,则只要保证在将下方序列反转后所得到的序列保证是单调上升的即可,则将其转化成了求最长上升子序列,对于lis问题可以用nlogn方法求解.

代码:

#include<cstdio>
#include<iostream>
#include<algorithm>

using namespace std;

const int maxn=100000+10;

int read()
{
    char ch=getchar();int ret=0,flag=1;
    while(ch<'0'||ch>'9') {if(ch=='-') flag=-1; ch=getchar();}
    while(ch>='0'&&ch<='9') ret=ret*10+ch-'0',ch=getchar();
    return ret*flag;
}

int nc[maxn],sc[maxn];//nc[i]/sc[i]表示i号位置的北方/南方城市编号,sc是反向存储的
int nrk[maxn],srk[maxn],dp[maxn];//nrk[i]表示i编号城市的北方位置号,srk是对应位置的城市在北方的位置号(反向存储)

int main()
{
    int n=read();
    for(int i=1;i<=n;i++) nc[i]=read(),nrk[nc[i]]=i;
    for(int i=1;i<=n;i++) sc[n-i+1]=read(),srk[n-i+1]=nrk[sc[n-i+1]];
    int ans=1;dp[1]=srk[1];
    for(int i=2;i<=n;i++) {//nlogn求最长上升子序列
        int k=upper_bound(dp+1,dp+ans+1,srk[i])-dp;
        dp[k]=srk[i];ans=max(ans,k);
    }
    printf("%d",ans);
    return 0;
}