枚举保留的那些数,可以先枚举前缀的gcd,然后因为一个数与log个和它不同的数取gcd后就会变成1,所以可以二分,这样是
nlog2n
的,为什么网上有人说两个log能过……
不过因为不同前缀gcd个数也只用log个,所以判一下当前枚举到的前缀gcd与上一位相同的话特殊处理就可以了……
T了两发就放弃了……不应该啊
#include <cstdio>
#include <iostream>
#include <algorithm>
using namespace std;
const int N=100010,M=30,P=998244353;
int a[N],n;
int _gcd[N][M],_log[N];
inline char nc(){
static char buf[100000],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline void rea(int &x){
char c=nc(); x=0;
for(;c>'9'||c<'0';c=nc());for(;c>='0'&&c<='9';x=x*10+c-'0',c=nc());
}
int gcd(int x,int y){
return y?gcd(y,x%y):x;
}
inline void Pre(){
for(int i=1;i<=n;i++) _log[i]=_log[i-1]+((1<<_log[i-1]+1)==i);
int t=_log[n];
for(int i=1;i<=n;i++) _gcd[i][0]=a[i];
for(int k=1;k<=t;k++)
for(int i=1;i+(1<<k)-1<=n;i++)
_gcd[i][k]=gcd(_gcd[i][k-1],_gcd[i+(1<<k-1)][k-1]);
}
inline int check(int l,int r){
int t=_log[r-l+1];
return gcd(_gcd[l][t],_gcd[r-(1<<t)+1][t]);
}
inline int query(int l,int r,int x){
int t=_log[r-l+1];
return (_gcd[l][t]%x==0)&&(_gcd[r-(1<<t)+1][t]%x==0);
}
int ans,subans[N],pre[N],suf[N];
int main(){
rea(n);
for(int i=1;i<=n;i++)
rea(a[i]),pre[i]=gcd(pre[i-1],a[i]);
for(int i=n;i;i--)
suf[i]=gcd(suf[i+1],a[i]);
Pre();
for(int i=0;i<n-1;i++)
if(i<3||pre[i]!=pre[i-1]){
int cur,pos;
if(i==0) cur=gcd(a[n],a[n-1]),pos=n-1;
else if(i==1) cur=gcd(a[1],a[n]),pos=n;
else cur=check(1,i),pos=n+1;
while(cur^1&&pos>i+1){
int l=i+2,r=pos-1,nxt=pos,mid;
while(l<=r){
mid=l+r>>1;
if(!query(mid,pos-1,cur)) l=mid+1;
else r=(nxt=mid)-1;
}
subans[i]=(1LL*(pos-nxt+1)*cur%P+subans[i])%P;
pos=nxt-1; cur=gcd(cur,a[pos]);
}
subans[i]=(1LL*(pos-i-1)*cur%P+subans[i])%P;
}
else subans[i]=(1LL*subans[i-1]+P-gcd(pre[i],suf[i+1]))%P;
for(int i=0;i<n-1;i++) ans=(ans+subans[i])%P;
printf("%d\n",(ans+check(1,n-1))%P);
return 0;
}