• codeforces632D. Longest Subsequence (最小公倍数)


    You are given array a with n elements and the number m. Consider some subsequence of a and the value of least common multiple (LCM) of its elements. Denote LCM as l. Find any longest subsequence of a with the value l ≤ m.

    A subsequence of a is an array we can get by erasing some elements of a. It is allowed to erase zero or all elements.

    The LCM of an empty array equals 1.

    Input

    The first line contains two integers n and m (1 ≤ n, m ≤ 106) — the size of the array a and the parameter from the problem statement.

    The second line contains n integers ai (1 ≤ ai ≤ 109) — the elements of a.

    Output

    In the first line print two integers l and kmax (1 ≤ l ≤ m, 0 ≤ kmax ≤ n) — the value of LCM and the number of elements in optimal subsequence.

    In the second line print kmax integers — the positions of the elements from the optimal subsequence in the ascending order.

    Note that you can find and print any subsequence with the maximum length.

    Examples
    input
    7 8
    6 2 9 2 7 2 3
    
    output
    6 5
    1 2 4 6 7
    
    input
    6 4
    2 2 2 3 3 3
    
    output
    2 3
    

    1 2 3

    题意:给你n个数,让你找出最长的子序列,使得子序列的最小公倍数小于等于m,并把序列顺序输出来。

    思路:这题容易想到用dp[i][j]表示当前遍历到第i个数,最小公倍数为j的最长子序列的长度是多少。状态转移为dp[i][lcm(a[i],k)]=dp[i-1][k]+1,但是这样的复杂度就爆了,所以要换种思路。我们可以枚举1~m中的每一个数作为最小公倍数,然后数出n个数中多少个数是它的因子,然后取最大就行了,这里要把a[i]的个数记下来,不然时间复杂度就爆。

    #include<iostream>
    #include<stdio.h>
    #include<stdlib.h>
    #include<string.h>
    #include<math.h>
    #include<vector>
    #include<map>
    #include<set>
    #include<queue>
    #include<stack>
    #include<string>
    #include<algorithm>
    #define inf 99999999
    #define pi acos(-1.0)
    #define maxn 1000050
    #define MOD 1000000007
    using namespace std;
    typedef long long ll;
    typedef long double ldb;
    int len[maxn],a[maxn],cnt[maxn];
    
    int main()
    {
        int n,m,i,j;
        while(scanf("%d%d",&n,&m)!=EOF)
        {
            for(i=1;i<=m;i++)len[i]=0,cnt[i]=0;
            for(i=1;i<=n;i++){
                scanf("%d",&a[i]);
                if(a[i]<=m)cnt[a[i] ]++;
            }
            for(i=1;i<=m;i++){
                for(j=i;j<=m;j+=i){
                    len[j]+=cnt[i];
                }
            }
            int maxx=0;int maxlen=0;
            for(i=1;i<=m;i++){
                if(len[i]>maxlen){
                    maxx=i;
                    maxlen=len[i];
                }
    
            }
            if(maxx==0){
                printf("1 0
    ");continue;
            }
    
            printf("%d %d
    ",maxx,maxlen);
            int flag=1;
            for(i=1;i<=n;i++){
                if(maxx%a[i]==0){
                    if(flag){
                        flag=0;printf("%d",i);
                    }
                    else printf(" %d",i);
                }
            }
            printf("
    ");
        }
        return 0;
    }
    



  • 相关阅读:
    近期前端中的 一些常见的面试题
    一道前端学习题
    前端程序员容易忽视的一些基础知识
    web前端工程师入门须知
    Web前端知识体系精简
    面试分享:一年经验初探阿里巴巴前端社招
    抽象类、抽象函数/抽象方法详解
    C#语法-虚方法详解 Virtual 虚函数
    面向对象语言:继承关系教程
    C#动态创建Xml-LinQ方式
  • 原文地址:https://www.cnblogs.com/herumw/p/9464536.html
Copyright © 2020-2023  润新知