#include <float.h>
#include <limits.h>
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <time.h>

// 内部定数
#define D_MOD			1000000007								// 除数(10の9乗+7)
#define D_ARRAY_MAX		10000									// 最大配列数
#define D_VAL_MAX		1000000									// 最大値
#define D_PRM_MAX		100										// 最大素数数
#define D_SEGT_CNT		1048576									// セグメントツリーデータ数(2の20乗)

// 内部構造体 - 素数情報
typedef struct Prm {
	int miVal;													// 値
	int miCnt;													// 個数
} Prm;

// 内部変数
static FILE *szpFpI;											// 入力
static int siACnt;												// 配列数
static int si1Array[D_ARRAY_MAX];								// 配列
static int si1Lcm[D_VAL_MAX];									// 公倍数
static Prm sz1Prm[D_PRM_MAX];									// 素数
static int siPCnt;												// 素数数
static int si1SegT[D_SEGT_CNT * 2];								// セグメントツリー
static int siSCNo;												// セグメントツリー - 子の開始位置

// 内部変数 - テスト用
#ifdef D_TEST
	static int siRes;
	static FILE *szpFpA;
	static int siTNo;
#endif

// １行出力
int
fOutLine(
	char *pcpLine				// <I> １行
)
{
	char lc1Buf[1024];

#ifdef D_TEST
	lc1Buf[0] = '\0';
	fgets(lc1Buf, sizeof(lc1Buf), szpFpA);
	if (strcmp(lc1Buf, pcpLine)) {
		siRes = -1;
	}
#else
	printf("%s", pcpLine);
#endif

	return 0;
}

// 最大値 - 取得
int
fGetMax(
	int piVal1					// <I> 値１
	, int piVal2				// <I> 値２
)
{
	if (piVal1 > piVal2) {
		return piVal1;
	}
	else {
		return piVal2;
	}
}

// 和 - 取得
int
fGetModSum(
	int piVal1					// <I> 値１
	, int piVal2				// <I> 値２
)
{
	return (piVal1 + piVal2) % D_MOD;
}

// 積 - 取得
int
fGetModMulti(
	int piVal1					// <I> 値１
	, int piVal2				// <I> 値２
)
{
	return (int)((long long)piVal1 * (long long)piVal2 % D_MOD);
}

// 最大公約数 - 取得
int
fGetGcd(
	int piVal1					// <I> 値１
	, int piVal2				// <I> 値２
)
{
	int liWork;

	// ユークリッドの互除法
	if (piVal1 < piVal2) {
		liWork = piVal1;
		piVal1 = piVal2;
		piVal2 = liWork;
	}
	while (piVal2 > 0) {
		liWork = piVal1 % piVal2;
		piVal1 = piVal2;
		piVal2 = liWork;
	}

	return piVal1;
}

// 素因数分解
int
fDivPrm(
	int piVal					// <I> 対象値
)
{
	int i;

	// 初期化
	memset(sz1Prm, 0, sizeof(sz1Prm));
	siPCnt = 0;

	// 上限
	int liLimit = (int)sqrt((double)piVal);

	// 素因数分解
	for (i = 2; i <= liLimit; i++) {

		// 分解終了
		if (piVal == 1) {
			break;
		}

		// 分解
		while (piVal % i == 0) {
			piVal /= i;
			sz1Prm[siPCnt].miCnt++;
		}

		// 分解あり
		if (sz1Prm[siPCnt].miCnt > 0) {
			sz1Prm[siPCnt].miVal = i;
			siPCnt++;
		}
	}

	// 上限を超える素数
	if (piVal != 1) {
		sz1Prm[siPCnt].miCnt = 1;
		sz1Prm[siPCnt].miVal = piVal;
		siPCnt++;
	}

	return 0;
}

// セグメントツリー - 子の開始位置 - セット
int
fSegTSetCSNo(
	int piCnt					// <I> 子データ数
)
{
	siSCNo = 1;
	while (siSCNo < piCnt) {
		siSCNo *= 2;
	}

	return 0;
}

// セグメントツリー - ２つの値の処理 - 積
int
fSegTVal2(
	int piVal1					// <I> 値１
	, int piVal2				// <I> 値２
)
{
	return fGetModMulti(piVal1, piVal2);
}

// セグメントツリー - 親データセット
int
fSegTSetPData(
)
{
	int i;

	// 最初の親番号範囲
	int liPNo1 = siSCNo / 2;
	int liPNo2 = siSCNo - 1;

	// 作成
	while (liPNo1 > 0) {
		for (i = liPNo1; i <= liPNo2; i++) {

			// セット
			si1SegT[i] = fSegTVal2(si1SegT[i * 2], si1SegT[i * 2 + 1]);
		}

		// 次の親番号範囲
		liPNo2 = liPNo1 - 1;
		liPNo1 /= 2;
	}

	return 0;
}

// セグメントツリー - 取得
int
fSegTGetRngMain(
	int piGetS					// <I> 取得範囲 - 開始 0～
	, int piGetE				// <I> 取得範囲 - 終了 0～
	, int piNNo					// <I> 現在位置 1～
	, int piNowS				// <I> 現在範囲 - 開始 0～
	, int piNowE				// <I> 現在範囲 - 終了 0～
)
{
	int liRet, liRet2;

	// 内包チェック
	if (piGetS <= piNowS && piNowE <= piGetE) {
		return si1SegT[piNNo];
	}

	// 中間位置
	int liCenter = (piNowS + piNowE) / 2;

	// 初期値
	liRet = 1;

	// 左側
	if (piGetS <= liCenter) {
		liRet2 = fSegTGetRngMain(piGetS, piGetE, piNNo * 2, piNowS, liCenter);
		liRet = fSegTVal2(liRet, liRet2);
	}

	// 右側
	if (piGetE >= liCenter + 1) {
		liRet2 = fSegTGetRngMain(piGetS, piGetE, piNNo * 2 + 1, liCenter + 1, piNowE);
		liRet = fSegTVal2(liRet, liRet2);
	}

	return liRet;
}
int
fSegTGetRng(
	int piGetS					// <I> 取得範囲 - 開始 0～
	, int piGetE				// <I> 取得範囲 - 終了 0～
)
{
	return fSegTGetRngMain(piGetS, piGetE, 1, 0, siSCNo - 1);
}

// 実行メイン
int
fMain(
)
{
	int i, j, k;
	char lc1Buf[1024];

	// データ - 初期化
	memset(si1Lcm, 0, sizeof(si1Lcm));
	memset(si1SegT, 0, sizeof(si1SegT));

	// 配列数 - 取得
	fgets(lc1Buf, sizeof(lc1Buf), szpFpI);
	sscanf(lc1Buf, "%d", &siACnt);

	// 配列 - 取得
	for (i = 0; i < siACnt; i++) {
		fscanf(szpFpI, "%d", &si1Array[i]);
	}
	fgets(lc1Buf, sizeof(lc1Buf), szpFpI);

	// 最大公約数 - 除算
	int liGcd = si1Array[0];
	for (i = 1; i < siACnt; i++) {
		liGcd = fGetGcd(liGcd, si1Array[i]);
	}
	for (i = 0; i < siACnt; i++) {
		si1Array[i] /= liGcd;
	}

	// 公倍数 - セット
	for (i = 0; i < siACnt; i++) {
		fDivPrm(si1Array[i]);
		for (j = 0; j < siPCnt; j++) {
			si1Lcm[sz1Prm[j].miVal] = fGetMax(si1Lcm[sz1Prm[j].miVal], sz1Prm[j].miCnt);
		}
	}

	// セグメントツリー - 子の開始位置 - セット
	fSegTSetCSNo(D_VAL_MAX);

	// セグメントツリー - 子データセット
	for (i = 0; i < D_VAL_MAX; i++) {
		si1SegT[siSCNo + i] = 1;
		for (j = 0; j < si1Lcm[i]; j++) {
			si1SegT[siSCNo + i] *= i;
		}
	}

	// セグメントツリー - 親データセット
	fSegTSetPData();

	// 和 - 取得
	int liSum = 0;
	for (i = 0; i < siACnt; i++) {
		fDivPrm(si1Array[i]);

		// 番兵追加
		sz1Prm[siPCnt].miVal = D_VAL_MAX - 1;
		sz1Prm[siPCnt].miCnt = 0;
		siPCnt++;

		// 積 - 取得
		int liNo = 0;
		int liMulti = 1;
		for (j = 0; j < siPCnt; j++) {
			int liPrm = sz1Prm[j].miVal;

			// 手前まで
			if (liNo < liPrm) {
				int liVal = fSegTGetRng(liNo, liPrm - 1);
				liMulti = fGetModMulti(liMulti, liVal);
			}

			// 自身
			int liCnt = si1Lcm[liPrm] - sz1Prm[j].miCnt;
			for (k = 0; k < liCnt; k++) {
				liMulti = fGetModMulti(liMulti, liPrm);
			}

			// 次の開始位置
			liNo = liPrm + 1;
		}

		// 和 - 加算
		liSum = fGetModSum(liSum, liMulti);
	}

	return liSum;
}

// １回実行
int
fOne(
)
{
	int liRet;
	char lc1Buf[1024];

	// 入力 - セット
#ifdef D_TEST
	sprintf(lc1Buf, ".\\Test\\T%d.txt", siTNo);
	szpFpI = fopen(lc1Buf, "r");
	sprintf(lc1Buf, ".\\Test\\A%d.txt", siTNo);
	szpFpA = fopen(lc1Buf, "r");
	siRes = 0;
#else
	szpFpI = stdin;
#endif

	// 実行メイン
	liRet = fMain();

	// 出力
	sprintf(lc1Buf, "%d\n", liRet);
	fOutLine(lc1Buf);

	// 残データ有無
#ifdef D_TEST
	lc1Buf[0] = '\0';
	fgets(lc1Buf, sizeof(lc1Buf), szpFpA);
	if (strcmp(lc1Buf, "")) {
		siRes = -1;
	}
#endif

	// テストファイルクローズ
#ifdef D_TEST
	fclose(szpFpI);
	fclose(szpFpA);
#endif

	// テスト結果
#ifdef D_TEST
	if (siRes == 0) {
		printf("OK %d\n", siTNo);
	}
	else {
		printf("NG %d\n", siTNo);
	}
#endif

	return 0;
}

// プログラム開始
int
main()
{

#ifdef D_TEST
	int i;
	for (i = D_TEST_SNO; i <= D_TEST_ENO; i++) {
		siTNo = i;
		fOne();
	}
#else
	fOne();
#endif

	return 0;
}



Measured memory_bytes: 12288