[51nod 1309]Value of all Permutations

题目大意

长度为N的整数数组A，有Q个查询，每个查询包含一个数M，对A的所有不同排列，执行find函数（需用到查询中的M），你来计算find函数的返回值的和。由于结果很大，输出Mod 1000000007的结果。

find(int permutation_A[], int M)
{
    x = Length(permutation_A)
    sum = 0
    for(i = 0; i < x; i++) {
        if (permutation_A[i] <= M)
            sum = sum + permutation_A[i]
        else
            break
    }
    return sum
}

做法

先不考虑相同元素，认为元素两两不同，共有n!种排列。
考虑计算ai对答案的贡献，假设>M的有k个，那么ai必须在这k个前面。
方案数$k!*(n-k-1)!*C_n^{k+1}=\frac{n!}{k+1}$
考虑算重的，对于每一种本质不同的排列，如果认为元素两两不同，就会算重每种元素数量的阶乘的积次，答案除以这个即可。
我们可以计算出本质不同排列数P，那么对于所有<=M的ai的和S，答案为$S*\frac{P}{k+1}$
通过前缀和与二分得到S和k。
注意P不一定整除k+1。

#include<cstdio>
#include<algorithm>
#include<cmath>
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
typedef long long ll;
const int maxn=50000+10,mo=1000000007;
int sum[maxn],a[maxn];
int i,j,k,l,t,n,m,p,ans;
int quicksortmi(int x,int y){
    if (!y) return 1;
    int t=quicksortmi(x,y/2);
    t=(ll)t*t%mo;
    if (y%2) t=(ll)t*x%mo;
    return t;
}
int main(){
    scanf("%d%d",&n,&m);
    fo(i,1,n) scanf("%d",&a[i]);
    sort(a+1,a+n+1);
    fo(i,1,n) sum[i]=(sum[i-1]+a[i])%mo;
    p=1;
    fo(i,1,n) p=(ll)p*i%mo;
    i=1;
    while (i<=n){
        j=i;
        while (j<=n&&a[j+1]==a[i]) j++;
        l=1;
        fo(k,1,j-i+1) l=(ll)l*k%mo;
        p=(ll)p*quicksortmi(l,mo-2)%mo;
        i=j+1;
    }
    while (m--){
        scanf("%d",&t);
        i=upper_bound(a+1,a+n+1,t)-a;
        k=n+1-i;
        ans=(ll)sum[i-1]*p%mo*quicksortmi(k+1,mo-2)%mo;
        (ans+=mo)%=mo;
        printf("%d\n",ans);
    }
}