#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 1e6+10;
const ll mod = 998244353;
int prime[N],mu[N],tot;
bool vis[N];
ll ksm(ll a,ll p)
{
	ll res=1;
	while(p)
	{
		if(p&1)
		res=res*a%mod;
		a=a*a%mod;
		p>>=1;
	}
	return res;
}
void sieve()
{
	memset(vis,1,sizeof vis);
	vis[0]=vis[1]=0;
	mu[1]=1;
	for(int i=2;i<N;i++)
	{
		if(vis[i])
		{
			prime[tot++]=i;
			mu[i]=-1;
		}
		for(int j=0;j<tot&&1ll*i*prime[j]<N;j++)
		{
			int x=prime[j];
			vis[i*x]=0;
			if(i%x==0)
			{
				mu[i*x]=0;
				break;
			}
			else
			mu[i*x]=-mu[i];
		}
	}
}
ll n,a[N],f[N],cnt[N],ma;
int main()
{
	sieve();
	scanf("%lld",&n);
	for(int i=1;i<=n;i++)
	{
		scanf("%lld",&a[i]);
		cnt[a[i]]++;
		ma=max(ma,a[i]);
	}
	for(int i=1;i<=ma;i++)
	{
		ll sum=0;
		for(int j=i;j<=ma;j+=i)
		{
			f[i]=(f[i]+cnt[j]*(cnt[j]-1)/2*j%mod*j%mod+sum*j*cnt[j]%mod)%mod;
			sum=(sum+j*cnt[j]%mod)%mod;
		}
	}
	ll ans=0;
	for(int i=1;i<=ma;i++)
	{
		ll now=0;
		for(int j=i;j<=ma;j+=i)
		now=(f[j]*mu[j/i]%mod+now)%mod;
		ans=(ans+now*ksm(i,mod-2))%mod;
	}
	printf("%lld\n",ans);
}