LOJ #6432. 「PKUSC2018」真实排名

题目在这里......

对于这道题，现场我写炸了......谁跟我说组合数O（n）的求是最快的？（~！@#￥￥%……&

#include <cstdio>
#include <algorithm>
using namespace std;

typedef long long ll;
;
;

ll ksm(ll x, ll n) {
    ll ans = 1ll; //printf("%lld ", x);
    while( n) {
        ) ans = ans * x % mod; x = x * x % mod;
        n >>= ;
    }
//    printf("ans = %lld, %lld\n", ans, x);
    return ans;
}
ll zt[N];
ll lucas(int m, int n) {
    ll x = zt[m] * zt[n-m] % mod; //printf("%lld\n", ksm(x, mod-2));
    ) % mod;
}
int data[N], a[N], n, k;
ll query1(int x) {
    , data++n, a[x]) - data;
    , data++n, (a[x] +  >> ) - ) - data - ;
    if( l + n - r < k) return 0ll;
    return lucas(k, l + n - r);
}

ll query2(int x) {
    , data++n, a[x]) - data;
    , data++n, a[x] << ) - data;
    if( l == r) r ++;
    if( r - l > k) return 0ll;
    return lucas(k-r+l, l + n - r);
}

int main() {
    zt[] = zt[] = 1ll;
    ; i <= N-; i ++) zt[i] = zt[i-] * i % mod;
//    for( int i = 1; i <= 10; i ++) printf("%d ", zt[i]); puts("");
    scanf("%d%d", &n, &k);
    ; i <= n; i ++) scanf("%d", a+i), data[i] = a[i];
    sort(data+, data++n);
    ; i <= n; i ++)
        printf("%lld\n", (query1(i) + query2(i)) % mod);
}

---

我们先复制一份排个序

我们考虑2种情况：

　　1：当前节点不变的情况，则大于他的节点和小于他的一半的节点都可以改变（显然），注意一下0，然后我们可以求得可以改变的数的个数t1。

　　2：当前节点要改变，那么想要维持他的排名不变，那么大于等于他且小于他的2倍的数都要改变（显然），然后我们可以求得可以改变的数t2。

然后就是求组合数（心态炸了QAQ）......

预处理一下，然后因为mod为质数，所以可以用费马小定理在log的时间内求得逆元......

每次logN的求t1， t2,。logmod的时间求逆元，所以复杂度是O(nlogn+logmod)。

费马小定理 ： a^(p-1) % p == 1 (p为质数);

组合数 ： 

（当时我脑抽用这个公式每次O(n)的求组合数只有45分QAQ，应该预处理阶乘）

---