Niuke IOIウィークリーコンテスト23-人気のグループD。一連の小さなL(dp状態+列挙係数)

元のタイトルリンク

アイデア:

この質問を見たとき、数論の知識を使って最適化する必要があると感じたので、最初に暴力を書くことにしました。

dp [i]が、i番目の項目が列挙されたときの最長の適切なシーケンス長を表すと仮定すると、転送は、80を超えるLISのn ^ 2転送を模倣します。

その後、最適化されません。

ユジュは、要因を列挙することで最適化できると述べました。

///dp[i]表示以数i为结尾的最长的好的序列长度
///dp1[i]表示因数包括i的数结尾的最长好的序列长度

この方法で転送する場合は、forの2番目のレベルをsqrt(n)の複雑さに変更し、プロセスでdpをdp1で更新します。因子を列挙した後、この数のdpを使用して因子のdp1を更新します。
ps:
1。最初の要件のため、元の配列をソートします
。2。係数が1の場合、転送は実行されないことに注意してください。

事実は、dpのステータスがいかに重要であるかを証明しています(エスケープ

コード:

#pragma GCC optimize(2)
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<ll,ll>PLL;
typedef pair<int,int>PII;
typedef pair<double,double>PDD;
#define I_int ll
inline ll read()
{
    
    
    ll x=0,f=1;
    char ch=getchar();
    while(ch<'0'||ch>'9')
    {
    
    
        if(ch=='-')f=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
    {
    
    
        x=x*10+ch-'0';
        ch=getchar();
    }
    return x*f;
}
#define read read()
#define closeSync ios::sync_with_stdio(0);cin.tie(0);cout.tie(0)
#define multiCase int T;cin>>T;for(int t=1;t<=T;t++)
#define rep(i,a,b) for(int i=(a);i<=(b);i++)
#define repp(i,a,b) for(int i=(a);i<(b);i++)
#define per(i,a,b) for(int i=(a);i>=(b);i--)
#define perr(i,a,b) for(int i=(a);i>(b);i--)
ll ksm(ll a,ll b,ll p)
{
    
    
    ll res=1;
    while(b)
    {
    
    
        if(b&1)res=res*a%p;
        a=a*a%p;
        b>>=1;
    }
    return res;
}
#define PI acos(-1)
#define x first
#define y second
const int maxn=1e6+7,inf=0x3f3f3f3f;
int n,a[maxn],dp[maxn],dp1[maxn],w[maxn];
int gcd(int a,int b){
    
    
    return b==0?a:gcd(b,a%b);
}
int main()
{
    
    
    int t=read;
    while(t--){
    
    
        n=read;
        rep(i,1,n) a[i]=read;
        sort(a+1,a+1+n);
///dp[i]表示以数i为结尾的最长的好的序列长度
///dp1[i]表示因数包括i的数结尾的最长好的序列长度
        memset(dp,0,sizeof dp);
        memset(dp1,0,sizeof dp1);
        int res=1;
        for(int i=1;i<=n;i++){
    
    
            dp[a[i]]=1;
            int idx=0;
            if(a[i]==1) continue;
            for(int j=1;j*j<=a[i];j++){
    
    
                if(a[i]%j==0){
    
    
                    int x1=a[i]/j,x2=j;
                    if(x1==x2){
    
    
                        if(x1!=1) dp[a[i]]=max(dp[a[i]],dp1[x1]+1),w[++idx]=x1;
                    }
                    else{
    
    
                        if(x1!=1) dp[a[i]]=max(dp[a[i]],dp1[x1]+1),w[++idx]=x1;
                        if(x2!=1) dp[a[i]]=max(dp[a[i]],dp1[x2]+1),w[++idx]=x2;
                    }
                }
            }
            for(int j=1;j<=idx;j++) dp1[w[j]]=max(dp1[w[j]],dp[a[i]]);///枚举a[i]的因子
        }
        for(int i=1;i<=n;i++) res=max(res,dp[a[i]]);
        printf("%d\n",res);
    }
    return 0;
}

おすすめ

転載: blog.csdn.net/weixin_45675097/article/details/114433520