【JZOJ4915】最长不下降子序列

Description

Data Constraint

Solution

我们看 D≤150 ，一定有循环节。

然后找到循环节，显然每个循环节至多取一个数。

然后剩余部分加上中间的一部分循环节做一次 lis ，加上剩余循环节的个数即为答案。

Code

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define fo(i,j,k) for(int i=j;i<=k;i++)
#define fd(i,j,k) for(int i=j;i>=k;i--)
#define N 151
#define M 100001
#define ll long long
using namespace std;
int a[M];
int pos[N];
int d[M];
int main()
{
    freopen("lis.in","r",stdin);
    freopen("lis.out","w",stdout);
    ll n;
    scanf("%lld",&n);
    int t0,A,B,C,D;
    scanf("%d %d %d %d %d",&t0,&A,&B,&C,&D);
    int t=t0,now=n;
    a[++a[0]]=t;
    fo(i,2,n)
    {
        t=(A*t*t+B*t+C)%D;
        if(pos[t]) {now=i;break;}
        pos[t]=i;
        a[++a[0]]=t;
    }
    int p=now-pos[t],qq=0,w=a[0];
    fo(i,1,p)
    {
        if(a[0]+p>n) break;
        fo(j,a[0]+1,a[0]+p) a[j]=a[j-p];
        a[0]+=p;
        qq+=p;
    }
    while(qq<=w)
    {
        if(a[0]+p>n) break;
        fo(j,a[0]+1,a[0]+p) a[j]=a[j-p];
        a[0]+=p;
        qq+=p;
    }
    ll ex=0;
    if(n>a[0]) ex=(n-a[0])/p;
    int q=(n-a[0])%p;
    int mx;
    fo(i,pos[t],pos[t]+q-1) a[++a[0]]=a[i],mx=max(mx,a[i]);
    int cd=0;
    fo(i,1,a[0])
    {
        int l=0,r=cd;
        while(l+1<r)
        {
            int mid=(l+r)/2;
            if(d[mid]<=a[i]) l=mid;
            else r=mid;
        }
        int zz;
        if(d[r]<=a[i]) zz=r;
        else zz=l;
        zz++;
        d[zz]=a[i];
        if(zz>cd) cd++;
    }
    printf("%lld",cd+ex);
}