#include <bits/stdc++.h>
using namespace std;

#define fopen                                   \
    freopen("E:/vscode/oi/in.txt", "r", stdin); \
    freopen("E:/vscode/oi/out.txt", "w", stdout);
#define ios                  \
    ios::sync_with_stdio(0); \
    cin.tie(0);

#define i64 long long
#define ull unsigned i64
#define pii pair<int, int>
#define pdd pair<double, double>
#define ld long double
#define ls(x) (x << 1)
#define rs(x) (x << 1 | 1)
#define lowbit(x) (x & -x)
#define de(x) cout << #x << " = " << x << '\n'
#define MAXP 20

const int N = 2e3 + 10, M = 3e3 + 10, mod = 998244353;
const double eps = 1e-12;
const double PI = acos(-1);

bool Mst;

#define int i64

const int pri[16] = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53};
vector<array<int, 2>> vec[M];
int qsm(int a, int b)
{
    int res = 1;
    while (b)
    {
        if (b & 1)
            res = res * a % mod;
        a = a * a % mod;
        b >>= 1;
    }
    return res;
}

int calc(int p)
{
    return (p - 1) * qsm(p, mod - 2) % mod;
}

void solve()
{
    int n;
    cin >> n;
    for (int id = 0; id < n; id++)
    {
        int x;
        cin >> x;
        int mask = 0, val = x;
        for (int i = 0; i < 16; i++)
        {
            int p = pri[i];
            if (x % p == 0)
            {
                while (x % p == 0)
                    x /= p;
                mask |= 1 << i;
            }
        }
        vec[x].push_back({mask, val});
    }
    vector<int> dp(1 << 17, 0);
    vector<int> f(1 << 16), g(16);
    for (int i = 0; i < 16; i++)
        g[i] = calc(pri[i]);

    for (int x = 0; x < (1 << 16); x++)
    {
        f[x] = 1;
        for (int i = 0; i < 16; i++)
            if (x >> i & 1)
                f[x] = f[x] * g[i] % mod;
    }

    dp[0] = 1;
    for (int p = 1; p < M; p++)
    {
        if (vec[p].empty())
            continue;
        int maxp = calc(p);
        bool isfir = true;
        for (auto [mask, val] : vec[p])
        {
            auto ndp = dp;
            if (p > 1)
                mask |= 1 << 16;
            for (int i = 0; i < (1 << 17); i++)
            {
                int j = i | mask, k = j ^ i;
                int res = dp[i] * val % mod;
                if (isfir && p > 1)
                    k |= 1 << 16;
                if (k >> 16 & 1)
                {
                    res = res * maxp % mod;
                    k ^= 1 << 16;
                }
                assert(k < (1 << 16));
                res = res * f[k] % mod;
                ndp[j] = (ndp[j] + res) % mod;
            }
            dp = ndp;
            isfir = false;
        }
    }
    int ans = 0;
    for (auto x : dp)
        ans = (ans + x) % mod;
    cout << ans << '\n';
}

bool Med;

signed main()
{
    ios;
    // fopen;
    cerr << (&Med - &Mst) / 1024 / 1024;

    int t = 1;
    srand(time(0));
    cout << fixed << setprecision(12);

    // cin >> t;
    for (int i = 1; i <= t; i++)
    {
        solve();
    }

    return 0;
}