\(LCM\) 只和每个质数的最大次数有关，所以考虑如果每个数都 \(\le 16\) 我们就可以状压每个质数出现的次数

\(f[i2][i3][i5][i7][i11][i13]\) 为 这 \(6\) 个质数的出现次数 这个状态的方案数，转移的话每一项都取 \(max\) 即可。

现在每个数都 \(\le 200\) ,然后需要用到一个套路，建议先做[NOI2015]寿司晚宴

一个数的 \(> \sqrt{200}\) 的质因数次数最多是 \(1\) ，且最多有一个\(> \sqrt{200}\) 的质因数,以下简称\(> \sqrt{200}\) 的质因数为 \(mx\)

我们设 \(dp[i2][i3][i5][i7][i11][i13]\) 为 这 \(6\) 个质数的出现次数 这个状态所有集合的 \(mx\) 的积的和

将数按照 \(> \sqrt{n}\) 的质因数 分类后，按照上面的题的套路，因为 \(mx\) 最多会算一次，所以我们用 \(f\) 表示先不管 \(mx\) 的转移，在最后再算出有贡献的一部分（容斥）后再乘上 \(mx\)

最后统计答案的时候直接枚举就行了，注意最后要 \(-1\)，因为不包括空集。

代码是真的毒瘤。。。

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#define LL long long
using namespace std;
int n, x;
LL ans;
const int N = 505, mod = 1e9 + 7;
int zhi[10] = {0, 2, 3, 5, 7, 11, 13}; //6
LL dp[8][5][4][3][3][3], f[8][5][4][3][3][3];
struct shu
{
	int mx, s[10];
	void Insert(int x)
	{
		for (int i = 1; i <= 6; ++i)
			if (!(x % zhi[i]))
				while (!(x % zhi[i]))++s[i], x /= zhi[i];
		mx = x;
	}
	friend bool operator <(const shu &a, const shu &b) {return a.mx > b.mx;}
} e[N];
signed main()
{
	cin >> n;
	for (int i = 1; i <= n; ++i)scanf("%d", &x), e[i].Insert(x);
	sort(e + 1, e + 1 + n); dp[0][0][0][0][0][0] = 1;
	for (int i = 1; i <= n; ++i)
	{
		if (e[i - 1].mx != e[i].mx || e[i].mx == 1)memcpy(f, dp, sizeof(dp));
		for (int i2 = 7; i2 >= 0; --i2)
			for (int i3 = 4; i3 >= 0; --i3)
				for (int i5 = 3; i5 >= 0; --i5)
					for (int i7 = 2; i7 >= 0; --i7)
						for (int i11 = 2; i11 >= 0; --i11)
							for (int i13 = 2; i13 >= 0; --i13)
								(f[max(i2, e[i].s[1])][max(i3, e[i].s[2])][max(i5, e[i].s[3])][max(i7, e[i].s[4])][max(i11, e[i].s[5])][max(i13, e[i].s[6])] += f[i2][i3][i5][i7][i11][i13]) %= mod;
		if (e[i + 1].mx != e[i].mx || e[i].mx == 1 || i == n)
			for (int i2 = 0; i2 <= 7; ++i2)
				for (int i3 = 0; i3 <= 4; ++i3)
					for (int i5 = 0; i5 <= 3; ++i5)
						for (int i7 = 0; i7 <= 2; ++i7)
							for (int i11 = 0; i11 <= 2; ++i11)
								for (int i13 = 0; i13 <= 2; ++i13)
									dp[i2][i3][i5][i7][i11][i13] += (f[i2][i3][i5][i7][i11][i13] - dp[i2][i3][i5][i7][i11][i13]) * e[i].mx % mod;
	}
	for (int i2 = 0, t2 = 1; i2 <= 7; ++i2, t2 *= 2)
		for (int i3 = 0, t3 = 1; i3 <= 4; ++i3, t3 *= 3)
			for (int i5 = 0, t5 = 1; i5 <= 3; ++i5, t5 *= 5)
				for (int i7 = 0, t7 = 1; i7 <= 2; ++i7, t7 *= 7)
					for (int i11 = 0, t11 = 1; i11 <= 2; ++i11, t11 *= 11)
						for (int i13 = 0, t13 = 1; i13 <= 2; ++i13, t13 *= 13)
							(ans += (LL)t2 * t3 % mod * t5 % mod * t7 % mod * t11 % mod * t13 % mod * dp[i2][i3][i5][i7][i11][i13]) %= mod;
	cout << ((ans - 1) % mod + mod) % mod;
	return 0;
}