• 最长不下降子序列


    这里写图片描述

    分析:
    在考场上的思路是对的:
    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;
    }
  • 相关阅读:
    vue中的Data为什么必须是一个函数
    单页面应用的优缺点
    数组去重
    mvvm框架
    前端计算精确度问题处理JS
    shell 修改json配置。
    ubuntu 两个文件夹合并
    fdisk、df与du的区别
    新买移动磁盘,使用前需要什么操作?
    Springboot+MybatisPlust+ControllerAdvice ;Mybatis_Plus多数据源,controller统一异常返回
  • 原文地址:https://www.cnblogs.com/wutongtong3117/p/7673450.html
Copyright © 2020-2023  润新知