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