import java.io.*;
import java.util.*;

public class Main {
    static final long MOD = 1000000007;

    public static void main(String[] args) throws Exception {
        // 使用 BufferedReader 以提高输入效率
        BufferedReader br = new BufferedReader(new InputStreamReader(System.in));
        PrintWriter out = new PrintWriter(System.out);
        StringTokenizer st = new StringTokenizer(br.readLine());

        int n = Integer.parseInt(st.nextToken());
        // 木棍长度范围是1~5000
        int MAX = 5000;
        long[] cnt = new long[MAX + 1];

        st = new StringTokenizer(br.readLine());
        for (int i = 0; i < n; i++){
            int x = Integer.parseInt(st.nextToken());
            cnt[x]++;
        }

        // 构造前缀和数组：prefix[i] = cnt[1] + cnt[2] + ... + cnt[i]
        long[] prefix = new long[MAX + 1];
        prefix[0] = 0;
        for (int i = 1; i <= MAX; i++){
            prefix[i] = prefix[i - 1] + cnt[i];
        }

        // 预先计算逆元，注意模数 MOD 为素数
        long inv2 = modPow(2, MOD - 2, MOD);  // 2 的逆元（常数 500000004）
        long inv6 = modPow(6, MOD - 2, MOD);  // 6 的逆元（166666668）

        long invalid = 0; // 无效（三根木棍不能组成三角形）的三元组数量
        // 枚举最大边 z（注意只有 z 必须出现在木棍中，即 cnt[z] > 0）
        for (int z = 1; z <= MAX; z++) {
            if (cnt[z] == 0) continue;
            long pairCount = 0; // 统计所有从长度 < z 中选出的 (x, y) 满足 x + y <= z 的组合数
            // x 的取值只需从 1 到 floor(z/2)
            for (int x = 1; x <= z / 2; x++) {
                int jmax = z - x; // y 的最大长度（必须满足 x + y <= z）
                // 若 jmax < x 则没有符合条件的 y，但当 x <= floor(z/2) 时这一条件一定满足
                // 1. 当 y = x 时：数量为 C(cnt[x],2)
                long waysSame = 0;
                if (cnt[x] >= 2) {
                    waysSame = (cnt[x] * (cnt[x] - 1)) % MOD;
                    waysSame = (waysSame * inv2) % MOD;
                }
                // 2. 当 y > x 时：枚举 y 从 x+1 到 jmax 的个数之和（可以用前缀和快速求和）
                long waysDiff = (cnt[x] * ((prefix[jmax] - prefix[x]) % MOD)) % MOD;
                long curr = (waysSame + waysDiff) % MOD;
                pairCount = (pairCount + curr) % MOD;
            }
            // 对于每个 z，所有 (x, y) 对与 cnt[z] 中任取1根组成一个无效三角形（x + y <= z）的组合
            invalid = (invalid + (pairCount * (cnt[z] % MOD)) % MOD) % MOD;
        }

        // 总的三元组数：C(n, 3) = n*(n-1)*(n-2)/6
        long totalTriples = (((((long)n * (n - 1)) % MOD) * (n - 2)) % MOD * inv6) % MOD;
        // 求概率：invalid/totalTriples mod MOD，模意义下使用乘以分母的逆元
        long result = (invalid * modPow(totalTriples, MOD - 2, MOD)) % MOD;
        out.println(result);
        out.flush();
        out.close();
    }

    // 快速幂取模（计算 base^exp mod mod）
    static long modPow(long base, long exp, long mod) {
        long result = 1;
        base %= mod;
        while(exp > 0) {
            if ((exp & 1) == 1) {
                result = (result * base) % mod;
            }
            base = (base * base) % mod;
            exp >>= 1;
        }
        return result;
    }
}