#include<bits/stdc++.h>
using namespace std;
const int N=1e6+1,mod=998244353;
int n,a[N],inv[N],w[N];
int main()
{
	inv[1]=1;
	for(int i=2;i<N;i++)
		inv[i]=(-1ll*inv[mod%i]*(mod/i)%mod+mod)%mod;
	scanf("%d",&n);
	int ans=0;
	for(int i=1;i<=n;i++)
	{
		int tmp;
		scanf("%d",&tmp);
		a[tmp]+=tmp;
		a[tmp]%=mod;
		ans-=tmp;
		ans%=mod;
	}
	for(int i=1;i<N;i++)
	{
		int tot=0;
		(w[i]+=inv[i])%=mod;
		for(int j=i+i;j<N;j+=i)
		{
			w[j]-=w[i],w[j]%=mod;
			tot+=a[j],tot%=mod;
		}
		(tot+=a[i])%=mod;
		(ans+=1ll*tot*tot%mod*w[i]%mod)%=mod;
	}
	ans=(1ll*ans*inv[2]%mod+mod)%mod;
	cout<<ans;
}
