日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

[BZOJ 2839]集合计数

發布時間:2024/4/17 编程问答 26 豆豆
生活随笔 收集整理的這篇文章主要介紹了 [BZOJ 2839]集合计数 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

Description

題庫鏈接

\(2^n\) 個集合,每個集合只包含 \([1,n]\) ,且這些集合兩兩不同。問有多少種選擇方法(至少選一個),使得這些集合交集大小為 \(k\)

\(0\leq k\leq n\leq 1000000\)

Solution

\(f(n)\) 為交集元素大于 \(k\) 的方案數,設 \(g(n)\) 為交集元素等于 \(k\) 的方案數。

容易得到

\[f(k)=\sum_{i=k}^n{i\choose k}g(i)\Rightarrow g(k)=\sum_{i=k}^n(-1)^{i-k}{i\choose k}f(i)\]

并且 \(f(i)={n\choose i}2^{2^{n-i}}\)

直接求就好了。

Code

#include <bits/stdc++.h> using namespace std; const int N = 1000000+5, yzh = 1000000007;int n, k, ifac[N], fac[N], ans;int quick_pow(int a, int b, int p) {int ans = 1;while (b) {if (b&1) ans = 1ll*ans*a%p;b >>= 1, a = 1ll*a*a%p;}return ans; } int C(int n, int m) {return 1ll*fac[n]*ifac[m]%yzh*ifac[n-m]%yzh; } void work() {scanf("%d%d", &n, &k);fac[0] = fac[1] = ifac[0] = ifac[1] = 1;for (int i = 2; i <= n; i++) ifac[i] = -1ll*yzh/i*ifac[yzh%i]%yzh;for (int i = 2; i <= n; i++)fac[i] = 1ll*fac[i-1]*i%yzh, ifac[i] = 1ll*ifac[i]*ifac[i-1]%yzh;for (int i = k; i <= n; i++)if ((i-k)&1) (ans -= 1ll*C(i, k)*C(n, i)%yzh*quick_pow(2, quick_pow(2, n-i, yzh-1), yzh)%yzh) %= yzh;else (ans += 1ll*C(i, k)*C(n, i)%yzh*quick_pow(2, quick_pow(2, n-i, yzh-1), yzh)%yzh) %= yzh;printf("%d\n", (ans+yzh)%yzh); } int main() {work(); return 0; }

轉載于:https://www.cnblogs.com/NaVi-Awson/p/9245549.html

總結

以上是生活随笔為你收集整理的[BZOJ 2839]集合计数的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。