7.27日比赛 总结+题解

小T的GCD

分数：10
预期分数：80~90（由于机子比较快，预期会高几分）

题意
求$GCD(a_l,a_{l+1}...,a_r)==1$ 的 $max(r-l+1)$

求$LCM(a_l,a_{l+1}...,a_r)== \prod_{k=l}^ra_k$ 的 $max(r-l+1)$

GCD直接暴力

LCM

1.预处理出所有数的因数

2.dp(单调队列优化）

f[i]表示以i为结尾的最大长度

$f[i]=min(f[i-1]+1,i-cnt)$

cnt表示最大的前缀

#include <bits/stdc++.h>
#define maxn 1000000
using namespace std;
int t,n,cases,ans,tot ;
int a[maxn+1],prime[maxn+1],mark[maxn+1],dp[maxn+1],pos[maxn+1];
vector <int> fac[maxn+1];
int read(){
    char c; 
    while (c=getchar(),c<'0' || c>'9'); 
    int x=c-'0';
    while (c=getchar(),c>='0' && c<='9') x=x*10+c-'0';
    return x;
}
int gcd(int x,int y) {
    if (!y) return x;
    else return gcd(y,x%y);
}
void init(){ //预处理因子 
    mark[1]=1;
    for (int i=2;i<=maxn;++i) {
        if (!mark[i]) {
            prime[++tot]=i;
            fac[i].push_back(i);
            for (int j=2;i*j<=maxn;++j){
                mark[j*i]=1;
                fac[i*j].push_back(i);
            }
        }
    }
}

int main() {
    init();
    t=read();
    while (t--){
        memset(pos,0,sizeof(pos));
        memset(dp,0,sizeof(dp));
        n=read();cases++;
        for (int i=1;i<=n;++i) a[i]=read();
        int g=a[1]; 
        printf("Case %d: ",cases);
        for (int i=1;i<=n;++i) g=gcd(g,a[i]);
        if (g>1) ans=-1; 
        else ans=n;
        printf("%d ",ans);
        ans=1;
        for (int i=1;i<=n;++i){
            int cnt=0;
            for (int j=0;j<fac[a[i]].size();++j){ //动态规划转移 
                cnt=max(cnt,pos[fac[a[i]][j]]);
                pos[fac[a[i]][j]]=i;
            }
            dp[i]=min(dp[i-1]+1,i-cnt);
            ans=max(ans,dp[i]);
        }
        if (ans==1) ans=-1;
        printf("%d\n",ans);
    }
    return 0;
}