#include #include #include #include using namespace std; const int FACTORIAL = 1000005; const int COLOR = 100005; const int MOD = 1000000007; int mul(int a, int b) { return 1LL * a * b % MOD; } int modpow(int p, int n) { if (n == 0) return 1; int t = modpow(p, n / 2); return n & 1 ? mul(t, mul(t, p)) : mul(t, t); } int inv(int v) { return modpow(v, MOD - 2); } int gcd(int a, int b) { return b == 0 ? a : gcd(b, a % b); } int factorial[FACTORIAL]; int invFactorial[FACTORIAL]; void makeFactorial() { factorial[0] = 1; for (int i = 1; i < FACTORIAL; ++i) { factorial[i] = mul(factorial[i - 1], i); invFactorial[i] = inv(factorial[i]); } } int nC, C[COLOR]; void read() { cin >> nC; for (int i = 0; i < nC; ++i) cin >> C[i]; } void work() { // divisor2cnt[C の約数]: C の約数が現れた回数 map divisor2cnt; for (int i = 0; i < nC; ++i) { int n = C[i]; for (int div = 1; div * div <= n; ++div) { if (n % div == 0) { if (div * div == n) { ++divisor2cnt[div]; } else { ++divisor2cnt[div]; ++divisor2cnt[n / div]; } } } } // 各 C が共通で持っている約数 vector commonDivisor; for (map::iterator it = divisor2cnt.begin(); it != divisor2cnt.end(); ++it) { if (it->second == nC) { commonDivisor.push_back(it->first); } } // toMul[i]: 0 ~ sum of nC - 1 の うち、 sum of nC と gcd をとって i // となる個数 int toMul[FACTORIAL] = {}; for (int i = 0, t = accumulate(C, C + nC, 0); i < t; ++i) ++toMul[gcd(t, i)]; int ans = 0; for (int loop = 0; loop < commonDivisor.size(); ++loop) { int divisor = commonDivisor[loop]; int total = 0; int curC[COLOR]; for (int i = 0; i < nC; ++i) { total += C[i] / divisor; curC[i] = C[i] / divisor; } // toAdd: total 個の花びらがあり、同じ色の花びらは curC[i] 個存在 // 並び方が何通り存在するか int toAdd = 1; toAdd = mul(toAdd, factorial[total]); for (int i = 0; i < nC; ++i) { toAdd = mul(toAdd, invFactorial[curC[i]]); } ans = (ans + mul(toMul[total], toAdd)) % MOD; } // C の総数で割る ans = mul(ans, inv(accumulate(C, C + nC, 0))); cout << ans << endl; } int main() { makeFactorial(); read(); work(); return 0; }