最长不下降子序列

这里写图片描述

分析:
在考场上的思路是对的:
D的范围很小,所以总有循环的一天,
这样我们处理出循环节,
需要注意的是,循环节不一定从1开始
比如说:1 4 2 5 3 2 5 3 循环节从第3位开始
我在考场上把序列缩到:不循环部分+一个循环节
长度设为n
在序列上做最长不下降子序列,
至于之后的循环节,他们对答案都可以贡献1

然而挂的很惨

看了一下大佬的ac程序
ta的n是这样的:
这里写图片描述
(我天,这么大,然而人家A了,你能说什么呢。。。)

好吧,n最大上了万位,
求最长不降只能用nlongn
然而我这个zz连最长不降nlogn都jj

看一下大佬优美的写法

int le=1;
f[1]=a[1];
for (i=2;i<=n;i++)
{
    if (a[i]>=f[le]) f[++le]=a[i];   //直接接上
    else  f[upper_bound(f+1,f+1+le,a[i])-f]=a[i];
    //upper_bound  >的第一个数 
}
return le;

%%%

最后要注意在%剩的余数中,也有可能对ans贡献1

这里写代码片
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#define ll long long

using namespace std;

ll nn;
int A,B,C,D,t;
int p[1000010],beg,xh,n;
ll a[1000010],f[1000010];

void doit()  //nlongn
{
    int i,j;
    ll ans=0;
    int le=1;
    f[1]=a[1];
    for (i=2;i<=n;i++)
    {
        if (a[i]>=f[le]) f[++le]=a[i];   //直接接上
        else  f[upper_bound(f+1,f+1+le,a[i])-f]=a[i];
        //upper_bound  >的第一个数 
    }
    if (n==nn)
    {
        printf("%d",le);    
        return;
    } 
    ans=(ll)le+(nn-n)/xh;
    int r=(nn-n)%xh;  //不构成循环的一小部分 
    if (r)
        for (i=0;i<r;i++)
            if (a[beg+r]>=f[le])  //f[le]是长度为le的最小结尾 
            {
                ans++;
                break;
            } 
    printf("%lld",ans);
}

void cl()
{
    memset(p,0,sizeof(p));
    int i,j;
    int x=t;
    p[t]=1;   //出现的位置 
    a[1]=t;
    for (i=2;i<=nn;i++)
    {
        x=A*x*x+B*x+C;
        x%=D;
        if (p[x]) 
        {
            n=i-1;
            beg=p[x];  //循环开始的地方 
            xh=i-p[x];
            break;
        }
        p[x]=i;
        a[i]=x;
    }
    n=xh*(xh<beg-1 ? beg-1:xh)+beg-1;  //循环部分:不循环部分  乘上较长的部分 
    if (n>nn) n=nn;
    for (i=beg+xh;i<=n;i++) a[i]=a[i-xh];
    for (int i=1;i<=beg+xh;i++) printf("%d ",a[i]);
    puts("");
    doit();
}

int main()
{
    //freopen("lis.in","r",stdin);  
    //freopen("lis.out","w",stdout);
    scanf("%lld",&nn);
    scanf("%d%d%d%d%d",&t,&A,&B,&C,&D);
    if (nn==1) printf("1
");
    else cl();
    return 0;
}
原文地址:https://www.cnblogs.com/wutongtong3117/p/7673450.html