agc023C – Painting Machines(組合數)

自為風月馬前卒發表於2019-02-22
GCAIMac

題意

題目連結

有(n)個位置，每次你需要以(1 sim n-1)的一個排列的順序去染每一個顏色，第(i)個數可以把(i)和(i+1)位置染成黑色。一個排列的價值為最早把所有位置都染成黑色的次數。問所有排列的分數之和

Sol

神仙題Orz

不難想到我們可以列舉染色的次數(i in [lceil frac{n}{2}
ceil, n – 1])。那麼問題轉化為求有多少排列是在(i)次之後把所有位置染成黑色（需要O(1)的複雜度）

我們把這個問題具體化一下，首先第一個位置必須要染，用去一次操作，此時還有(t = n-2)個位置沒有染。對於一次染色，我們可以把兩個位置染成黑色，也可以把一個位置染成黑色。那麼問題轉化為每次可以讓(t)(-1)或者(-2)。問在(i – 1)操作後(t = 0)的方案數。

這個問題好像還是不好做，但是我們可以去求至多(i-1)次操作後(t=0)的方案。可以先把每次的(-1)算進去，這樣每次就變成了(-0/-1)那麼只要在剩下的(i-1)次操作中有((n-2)-(i-1))次選了(-1)就一定可行。

因此方案數為(C_{i-1}^{n-i-1})，最後每次計算的時候在乘上排列的係數就行了

#include<bits/stdc++.h>
#define LL long long 
using namespace std;
const int MAXN = 1e6 + 10, INF = 1e9 + 7, mod = 1000000007;
template<typename A, typename B> inline void chmax(A &x, B y) {
    x = x > y ? x : y;
}
template<typename A, typename B> inline void chmin(A &x, B y) {
    x = x < y ? x : y;
}
template<typename A, typename B> inline int mul(A x, B y) {
    return 1ll * x * y % mod;
}
template<typename A, typename B> inline void add2(A &x, B y) {
    if(x + y < 0) x = x + y + mod;
    else x = (x + y >= mod ? x + y - mod : x + y);
}
template<typename A, typename B> inline int add(A x, B y) {
    if(x + y < 0) return x + y + mod;
    else return x + y >= mod ? x + y - mod : x + y;
}
template<typename A, typename B> inline int fp(A a, B p) {
    int base = 1;
    while(p) {
        if(p & 1) base = mul(base, a);
        a = mul(a, a); p >>= 1;
    }
    return base;
}
inline int read() {
    char c = getchar(); int x = 0, f = 1;
    while(c < `0` || c > `9`) {if(c == `-`) f = -1; c = getchar();}
    while(c >= `0` && c <= `9`) x = x * 10 + c - `0`, c = getchar();
    return x * f;
}
int N, fac[MAXN], ifac[MAXN], f[MAXN];
int C(int N, int M) {
    if(M > N) return 0;
    return mul(fac[N], mul(ifac[M], ifac[N - M]));
}
int main() {
    N = read(); fac[0] = 1;
    for(int i = 1; i <= N; i++) fac[i] = mul(i, fac[i - 1]);
    ifac[N] = fp(fac[N], mod - 2);
    for(int i = N; i >= 1; i--) ifac[i - 1] = mul(ifac[i], i);
    int down = N / 2 + (N & 1);
    for(int i = down; i < N; i++) 
        f[i] = mul(mul(C(i - 1, N - i - 1), fac[i]), fac[N - i - 1]);
    int ans = 0;
    for(int i = down; i < N; i++) add2(ans, mul(i, add(f[i], -f[i - 1])));
    cout << ans;
    return 0;
}
/*
3
0 1 1
5 7 3
*/

相關文章