Description有 n n n 只变色龙,一开始都是蓝色。现在你喂了 k k k 次球,每次指定一只变色龙吃下你指定颜色的球。
一只变色龙从蓝色变成红色当且仅当它吃的红球比蓝球多; 一只变色龙从红色变成蓝色当且仅当它吃的蓝球比红球多。
求最后能使所有变色龙都变成红色的方案数。
两个方案不同当且仅当至少一次喂的球颜色不同(而不是喂的变色龙不同 )。
注意:存在一次喂的变色龙不同的两个方案可能是相同的方案。
1 ≤ n , k ≤ 5 × 1 0 5 1\leq n,k\leq 5\times 10^5 1 ≤ n , k ≤ 5 × 1 0 5 。
Solution考虑对于一只变色龙,怎样给他喂球才能变成红色。
设喂了 a a a 个红球,b b b 个蓝球。
a > b a>b a > b :为红色。a < b a<b a < b :为蓝色。a = b a=b a = b :为不同于最后喂的颜色的另一个颜色,即最后的颜色异或 1 1 1 。首先如果 a − b ≥ 2 a-b\geq 2 a − b ≥ 2 ,那么把多出来的 a − b − 1 a-b-1 a − b − 1 个球给别人一定不劣,所以这里只需要满足 a = b + 1 a=b+1 a = b + 1 或 a = b a=b a = b 且最后一次选了蓝球。
不妨设一个方案总共放了 R R R 个红球,B B B 个蓝球,则要满足条件一定要满足 R ≥ B R\geq B R ≥ B 。
容易发现 a = b + 1 a=b+1 a = b + 1 的变色龙总共 R − B R-B R − B 个,a = b a=b a = b 的变色龙共 n − ( R − B ) n-(R-B) n − ( R − B ) 个,不妨设 c n t = n − ( R − B ) cnt=n-(R-B) c n t = n − ( R − B ) 。
所以判断一个方案是否合法等价于能否在其中找到 c n t cnt c n t 个形如 RB 的匹配。
但是这个判定还不够简单。
注意到这些匹配一定是前 c n t cnt c n t 个 R 和后 c n t cnt c n t 个 B 按顺序进行匹配,那么对于倒数第 i i i 个 B,它前面一定有至少 c n t − i + 1 cnt-i+1 c n t − i + 1 个 R。
所以对于每个前缀,他的 R 的个数 ≥ \geq ≥ B 的个数 − ( B − c n t ) -(B-cnt) − ( B − c n t ) 。
那么把 R 看作 + 1 +1 + 1 ,B 看作 − 1 -1 − 1 做一遍前缀和,只要满足每个前缀和 ≥ R − n \geq R-n ≥ R − n 即可。
这就转化为了一个类似卡特兰数的东西,答案就是:
( R + B R ) − ( R + B n + B − R − 1 ) \binom{R+B}{R}-\binom{R+B}{n+B-R-1} ( R R + B ) − ( n + B − R − 1 R + B )
这里对于 R = B R=B R = B 要特判,因为最后一步必须是 B B B ,所以只要把 B − 1 B-1 B − 1 再做上面那个式子即可。
时间复杂度:O ( k ) O(k) O ( k ) 。
Code1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 #include <bits/stdc++.h> const int kMaxN = 5e5 + 5 , kMod = 998244353 ;int n, k;int inv[kMaxN], fac[kMaxN], ifac[kMaxN];int add (int x, int y) { return (x + y) >= kMod ? (x + y - kMod) : (x + y); }int sub (int x, int y) { return (x < y) ? (x - y + kMod) : (x - y); }int C (int m, int n) { if (m < n || m < 0 || n < 0 ) return 0 ; return 1ll * fac[m] * ifac[n] % kMod * ifac[m - n] % kMod; }void prework () { fac[0 ] = ifac[0 ] = fac[1 ] = ifac[1 ] = inv[1 ] = 1 ; for (int i = 2 ; i <= k; ++i) { inv[i] = 1ll * (kMod - kMod / i) * inv[kMod % i] % kMod; fac[i] = 1ll * i * fac[i - 1 ] % kMod; ifac[i] = 1ll * inv[i] * ifac[i - 1 ] % kMod; } }void dickdreamer () { std::cin >> n >> k; if (k < n) return void (std::cout << "0\n" ); prework (); int ans = 0 ; for (int i = (k + 1 ) / 2 ; i <= k; ++i) { int j = k - i; if (i == j) -- j; ans = add (ans, sub (C (i + j, i), C (i + j, n + j - i - 1 ))); } std::cout << ans << '\n' ; }int32_t main () {#ifdef ORZXKR freopen ("in.txt" , "r" , stdin); freopen ("out.txt" , "w" , stdout);#endif std::ios::sync_with_stdio (0 ), std::cin.tie (0 ), std::cout.tie (0 ); int T = 1 ; while (T--) dickdreamer (); return 0 ; }