51nod1670 打怪兽 期望，数学

题目大意：

lyk在玩一个叫做“打怪兽”的游戏。
游戏的规则是这样的。
lyk一开始会有一个初始的能量值。每次遇到一个怪兽，若lyk的能量值>=怪兽的能量值，那么怪兽将会被打败，lyk的能量值增加1，否则lyk死亡，游戏结束。
若怪兽全部打完，游戏也将会结束。
共有n个怪兽，由于lyk比较弱，它一开始只有0点能量值。
n个怪兽排列随机，也就是说共有n!种可能，lyk想知道结束时它能量值的期望。
由于小数点比较麻烦，所以你只需要输出期望*n!关于1000000007取模后的值就可以了！

例如有两个怪兽，能量值分别为{0,1}，那么答案为2，因为游戏结束时有两种可能，lyk的能量值分别为0和2。期望为1,1*2!=2，所以答案为2。

Input

第一行一个数n(1<=n<=100000)。
接下来一行n个数ai表示怪兽的能量(0<=ai<n)。

Output

一行表示答案

Input示例

2
0 1

Output示例

2

思路：首先将问题转化，因为是期望*n！，所以结果等价于死亡时能量值的和，死亡时能量值可能为0~n，所以扫一遍，每次都记录在此时死亡的情况有多少种。如果在i位置死亡，则得分为这个i ，个数为或者走到这里的可能情况ans1 乘以（n-i-1)的阶乘。最后求和得出ans即可。若不能货到最后，在哪里结束标记一下，跳出循环后计算即可。

#include <iostream>
#include <cstdio>
#include <string>
#include <cstring>
#include <fstream>
#include <algorithm>
#include <cmath>
#include <queue>
#include <stack>
#include <vector>
#include <map>
#include <set>
#include <iomanip>

using namespace std;
#define maxn 100005
#define MOD 1000000007
#define mem(a , b) memset(a , b , sizeof(a))
#define LL long long
#define INF 1000000000
int a[maxn];
int vis[maxn];
int jc[maxn];
int n ;

void init()
{
    jc[0]= 1;
    jc[1] = 1;
    LL tmp = 1;
    for(int i = 2 ; i < maxn ; i ++)
    {
        tmp = jc[i-1];
        tmp *= i;
        tmp %= MOD;
        jc[i] = tmp;
    }
}

void solve()
{
    int tmp = a[0];
    int num = 1;
    for(int i = 1 ; i < n ; i ++)
    {
        if(tmp == a[i]) num++;
        else
        {
            vis[tmp] = num;
            num = 1;
            tmp = a[i];
        }
    }
    vis[tmp] = num;
    for(int i = 1 ; i < n ; i ++) vis[i] += vis[i-1] ;
}

int main()
{
    init();
    while(scanf("%d" , &n) != EOF)
    {
        mem(vis , 0);
        for(int i = 0 ; i < n ; i ++) scanf("%d" , &a[i]);
        sort(a , a + n);
        if(a[0] > 0)
        {
            printf("0\n");
            continue;
        }
        solve();
        LL ans1 = 1 , tmp , ans = 0;
        LL num = 1 ;int flag = 0;
        for(int i = 0 ; i < n ; i ++)
        {
            if(vis[i] - i <= 0)
            {
                flag = i;
               // cout << "YES" << flag << endl;
                break;
            }
            tmp = n - vis[i];
            ans += ( (ans1*tmp)%MOD* jc[n-i-1] % MOD*i) % MOD;
            ans % MOD;
            ans1 *= (vis[i] - i);
            ans1 %= MOD;
        }
        if(!flag)
        {
            ans1 *= n;
            ans1 %= MOD;
            ans = (ans + ans1) % MOD;
        }
        else ans += ((ans1 * jc[n-flag]) % MOD * flag ) % MOD;
        printf("%lld\n" , ans%MOD );
    }
    return 0;
}