Codeforces 938E组合数学逆元

题目链接：http://codeforces.com/contest/938/problem/E

题意：给你n个数字，求所有排列的fa的和：

对于一个排列来说，一开始fa为1，m=1，若a[m]<a[i] fa=fa+a[m]。

题解：我们考虑一个数的贡献，若这个数能产生贡献，则比这个数大的数都在这个数后面（包括相等），那么我们就可以得到如下能得到贡献的数目。

意思为：从n个位置中挑出n-i+1（大于等于这个数的个数）个位置，这n-i+1个位置中这个数必须是第一位，其他数可以任意排列。然后把剩下的小于他的数插入剩下的空位。

化简之后得到：

这样就可以通过逆元O（logn）计算当前这个数的答案了，对于相同的数，只要乘上其个数就可以了。

#include <bits/stdc++.h>

using namespace std;
using ll = long long ;
using ld = long double;

#define pb push_back
#define SZ(X) ((int)X.size())
#define mp make_pair
#define Fi first
#define Se second
#define pii pair<int,int>
#define pll pair<ll,ll>
#define pli pair<ll,int>
#define pil pair<int,ll>
#define ALL(X) X.begin(),X.end()
#define RALL(X) X.rbegin(),X.rend()
#define rep(i,j,k) for(int i = j;i <= k;i ++)
#define per(i,j,k) for(int i = j;i >= k;i --)
#define mem(a,p) memset(a,p,sizeof(a))


const int N = 1E6 + 7;
const int MOD = 1E9 + 7;
ll a[N];
ll inv[N];
ll M;

ll qmod(ll n, ll k) // n ^k
{
    ll res = 1;
    while(k) {
        if(k & 1) res = res * n % MOD;
        n = n * n % MOD;
        k >>= 1;
    }
    return res;
}
int main()
{
    inv[1] = 1;
    for(int i = 2;i < N;i ++) {
        inv[i] = (MOD - MOD / i) * inv[MOD % i] % MOD;
    }
//    for(int i = 1;i <= 100;i ++) {
//        printf("%lld %lld\n",inv[i],qmod(i,MOD-2));
//    }
    int n;
    scanf("%d", &n);
    for(int i = 1;i <= n;i ++) scanf("%lld", &a[i]);
    M = 1;
    for(int i = 1;i <= n;i ++) M = M * i % MOD;
    sort(a + 1, a + 1 + n);
    ll res = 0;
    int j;
    for(int i = 1;i < n;i = j) {
        int num = 0;
        for(j = i;a[j] == a[i];j ++) num ++;
        if(a[i] == a[n]) break;
        res = (res + a[i] * M % MOD * num % MOD * inv[n - i + 1] % MOD);
        if(res >= MOD)  res %= MOD;
    }
    printf("%lld\n",res);
    return 0;
}