#include<bits/stdc++.h>

using namespace std;

#define    pf               printf
#define    PF(a)            printf("%d\n",(a))
#define    PF2(a,b)         printf("%d %d\n",(a), (b))
#define    PF3(a,b,c)       printf("%d %d %d\n",(a), (b), (c))
#define    PFL(a)           printf("%lld\n",(a))
#define    PFL2(a,b)        printf("%lld %lld\n",(a), (b))
#define    PFL3(a,b,c)      printf("%lld %lld %lld\n",(a), (b), (c))
#define    get(a)           cin>>a
#define    sf               scanf
#define    SF(a)            scanf("%d",&a)
#define    SF2(a,b)         scanf("%d %d",&a, &b)
#define    SF3(a,b,c)       scanf("%d %d %d",&a, &b, &c)
#define    SFL(a)           scanf("%lld",&a)
#define    SFL2(a,b)        scanf("%lld %lld",&a, &b)
#define    SFL3(a,b,c)      scanf("%lld %lld %lld",&a, &b, &c)
#define    gl(a)            getline(cin,(a))
#define    gc()             getchar()
#define    pb               push_back
#define    m_p              make_pair
#define    pc()             printf("Case %d: ",tt++)
#define    tc()             cout<<"Case "<<tt++<<": "
#define    SZ(a)            (int)(a).size()
#define    all(a)           a.begin(),a.end()
#define    FF               first
#define    SS               second
#define    MAX              2134567891
#define    MOD              1000000007
#define    MM               100005
#define    PI               2*acos(0.0)
#define    cond(n,m)        (xx>=0 && xx<(n) && yy>=0 && yy<(m))

#define    mem(a)           memset((a),0,sizeof (a))
#define    SET(a)           memset((a),-1,sizeof (a))
#define    output           freopen("output.txt","w",stdout);
#define    input            freopen("input.txt","r",stdin);
#define    I_O              ios_base::sync_with_stdio(0); cin.tie(0)
#define    rep(a)           for(int i=0;i<(a);i++)
#define    FOR(a,b)         for(int i=(a);i<(b);i++)
#define    REP(a)           for(int j=0;j<(a);j++)
#define    rev(a)           for(int i=(a);i>=0;i--)
#define    fr(i,n)          for(i=0;i<n;i++)
#define    fr1(i,n,m)       for(i=n;i<m;i++)

int dx[]  = {0, 0, 1, -1};
int dy[]  = {1, -1, 0, 0};
int dx8[] = {0, 0, 1, 1, 1, -1, -1, -1};
int dy8[] = {1, -1, -1, 0, 1, -1, 0, 1};
int kx8[] = {1, 1, 2, 2, -1, -1, -2, -2};
int ky8[] = {2, -2, 1, -1, 2, -2, 1, -1};

typedef long long ll;
typedef unsigned long long llu;
typedef priority_queue < int > PQ;
typedef pair < int , int > pi;
typedef pair < int , pi > pii;
typedef vector < int > vi;
typedef vector < ll > vl;
typedef vector < string > vs;

ll bigmod(ll a, ll b, ll c)
{
    if(b==0)return 1%c;ll x=bigmod(a,b/2,c);x=x*x%c;
    if(b%2==1)x=x*a%c;return x;
}
ll mod_inverse(ll a, ll mod){return bigmod(a,mod-2,mod);}
/*
//FUNCTION//

ll poww(ll a, ll b)
{
    if(b==0)return 1;ll x=poww(a,b/2);x=x*x;if(b%2==1)x=x*a;return x;
}


ll LCM(ll a, ll b){ return a*b/ __gcd(a,b);}

int pr = 1000000;
vi primelist;
bool a[MM*10];
void seive( )
{
    int i,j,k=sqrt(pr);
    a[1]=1;
    primelist.pb(2);
    for(i=4;i<=pr;i+=2)a[i]=1;
    for(i=3;i<=k;i+=2)if(!a[i])for(j=i*i;j<=pr;j+=2*i)a[j]=1;
    for(i=3;i<=pr;i+=2)if(!a[i])primelist.pb(i);
}

ll fact_divs( ll n, ll p)
{
    ll cnt=0;while(p<=n){cnt += n/p;n /= p;}return cnt;
}
struct point{double x, y;};
double area( point A, point B, point C){
return 0.5*(A.x*(B.y-C.y) - A.y*(B.x-C.x) + B.x*C.y - C.x*B.y);
}


ll Set(ll N,ll pos){return N=N | (1LL<<pos);}
ll reset(ll N,ll pos){return N= N & ~(1LL<<pos);}
bool check(ll N,ll pos){return (bool)(N & (1LL<<pos));}

//*/
const int maxn = 1000005;

ll inv_md[maxn+2];
void inverse_precal(ll m)
{
    //m is prime;
    inv_md[1] = 1;
    for(ll i = 2; i < maxn; ++i)
    inv_md[i] = (m - (m/i) * inv_md[m%i] % m) % m;
}

ll inv(ll a, ll b)
{
    return inv_md[a];
}

vector<int> F[maxn];
ll p[maxn];
ll Count[maxn];

void sieve(){
    for(int i=1; i<maxn; i++){
        for(int j=i; j<maxn; j+=i){
            F[j].push_back(i);
        }
    }
}

int a[maxn];
ll mod = 998244353;
int main(){
    sieve();
    inverse_precal(mod);
    int n;
    scanf("%d",&n);
    ll tot = 0;
    for(int i=1; i<=n; i++){
        scanf("%d",&a[i]);
        tot += a[i];
        tot %= mod;
        for(int x : F[a[i]])
        {
            p[x] += a[i];
            if(p[x]>=mod)p[x] -= mod;
        }
    }

    for(int i=maxn; i>=1; i--){
        if(p[i]==0)continue;
        Count[i] += (p[i] * p[i])%mod;
        Count[i] %= mod;
        if(Count[i]<0)Count[i] += mod;
        ll md_inv = inv(i,mod);
        Count[i] = (Count[i]*md_inv)%mod;
        for(auto j:F[i])
        {
            if(i==j)break;
            Count[j] -= (Count[i]*i)%mod;
            if(Count[j]<0) Count[j] += mod;
        }
    }

    ll sum = 0;
    for(int i=1; i<maxn; i++)
    {
        sum += Count[i];
        sum %= mod;
    }
    sum -= tot;
    if(sum<0)sum += mod;
    ll md_inv = inv(2,mod);
    sum *= md_inv;
    sum %= mod;
    PFL(sum);

    return 0;
}


Measured memory_bytes: 140928