zoukankan      html  css  js  c++  java
  • ACM数论之旅17---反演定理 第一回 二项式反演(神说要有光 于是就有了光(´・ω・`))

    终于讲到反演定理了,反演定理这种东西记一下公式就好了,反正我是证明不出来的~(~o ̄▽ ̄)~o

    首先,著名的反演公式

    我先简单的写一下o( ̄ヘ ̄*o)

    比如下面这个公式

    f(n) = g(1) + g(2) + g(3) + ... + g(n)

    如果你知道g(x),蓝后你就可以知道f(n)了

    如果我知道f(x),我想求g(n)怎么办

    这个时候,就有反演定理了

    反演定理可以轻松的把上面的公式变为

    g(n) = f(1) + f(2) + f(3) + ... + f(n)

    当然,我写的只是个形式,怎么可能这么简单。◕‿◕。

    其实每一项再乘一个未知的函数就对了,但是这个函数我们不知道(不用担心,数学家已经帮我们解决了,我们直接用就可以了)

    反演公式登场( >ω<)

    反演定理1

    c和d是两个跟n和r有关的函数

    根据用法不同,c和d是不同的

    一般数学家会先随便弄c函数

    然后经过复杂的计算和证明,得到d函数

    然后公式就可以套用了

    正片开始

    二项式反演公式

    二项式反演1

    那个括号起来的就是组合数,我记得组合数那章我有说过

    组合数4

    二项式反演也就是记住这个公式就算结束了

    然后我们开始实战(/ω\)

    容斥那章讲过的全错排(装错信封问题)

    hdu 1465

    http://acm.hdu.edu.cn/showproblem.php?pid=1465

     设g(i)表示正好有i封信装错信封

    那么全部的C(n, i)*g(i)加起来正好就是所有装信的情况,总共n!种情况

    n! = Σ C(n, i)*g(i) (i从0到n)

    那么f(n) = n!,所以f(x) = x!

    那么我们要求g(n)

    根据公式

    g(n) = Σ (-1)^(n-i) * C(n, i) * f(i)  (i从0到n)

    那么就可以计算啦~(≧▽≦)/~

    AC代码:

     1 #include<cstdio>
     2 typedef long long LL;
     3 int n, flag;
     4 LL fac[25];
     5 LL ans;
     6 void init(){
     7     fac[0] = 1;
     8     for(int i = 1; i <= 20; i ++) fac[i] = fac[i-1] * i;    
     9 }
    10 int main(){
    11     init();
    12     while(~scanf("%d", &n)){
    13         ans = 0;
    14         flag = n & 1 ? -1 : 1;//起始符号
    15         for(int i = 0; i <= n; i ++){
    16             ans += flag * fac[n] / fac[n-i];
    17             flag = -flag;  
    18         }
    19         printf("%I64d
    ", ans);
    20     }
    21 }
    View Code

    是不是很好用但是不容易想到T_T

    这也没有办法

    再来一题吧

    还是容斥那一章讲过的题目的

    UVALive 7040

    https://icpcarchive.ecs.baylor.edu/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=5052

    题意:给n盆花涂色,从m种颜色中选取k种颜色涂,保证正好用上k种颜色,你必须用上这k种颜色去涂满n个相邻的花,并且要求相邻花的颜色不同,求方案数。

     (1 ≤ n, m ≤ 1e9 , 1 ≤ k ≤ 1e6 , k ≤ n, m)

    首先,用k种颜色涂花,假如不考虑全部用上,那么总的方案数是多少

    第一盆花有k种颜色选择,之后的花因为不能跟前一盆花的颜色相同,所以有k-1种选择

    于是总方案数为k*(k-1)^(n-1)

    我们设必须用 i 种颜色两两不相邻的涂花的方案数为 g(i)

    那么

    k*(k-1)^(n-1) = Σ C(k, i)*g(i) (i从1到k)

    令f(k) = k*(k-1)^(n-1),那么f(x) = x*(x-1)^(n-1)

    二项式反演公式出现了

    所以g(k) = Σ (-1)^(k-i) * C(k, i) * f(i)  (i从1到k)

    最终的答案就是C(m, k) * g(k)

    (这里m有1e9,C(m, k)直接用for循环算,直接for循环从m*(m-1)*...*(m-k+1)再乘k的阶乘的逆元)

    AC代码:

     1 #include<cstdio>
     2 typedef long long LL;
     3 const int N = 1000000 + 5;
     4 const int MOD = (int)1e9 + 7;
     5 int F[N], Finv[N], inv[N];
     6 LL pow_mod(LL a, LL b, LL p){ 
     7     LL ret = 1;
     8     while(b){
     9         if(b & 1) ret = (ret * a) % p;
    10         a = (a * a) % p;
    11         b >>= 1;
    12     }
    13     return ret;
    14 }
    15 void init(){
    16     inv[1] = 1;
    17     for(int i = 2; i < N; i ++){
    18         inv[i] = (MOD - MOD / i) * 1ll * inv[MOD % i] % MOD;
    19     }
    20     F[0] = Finv[0] = 1;
    21     for(int i = 1; i < N; i ++){
    22         F[i] = F[i-1] * 1ll * i % MOD;
    23         Finv[i] = Finv[i-1] * 1ll * inv[i] % MOD;
    24     }
    25 }
    26 int comb(int n, int m){
    27     if(m < 0 || m > n) return 0;
    28     return F[n] * 1ll * Finv[n - m] % MOD * Finv[m] % MOD;
    29 }
    30 int main(){
    31     init();
    32     int T, n, m, k, ans, flag, temp;
    33     scanf("%d", &T);
    34     for(int cas = 1; cas <= T; cas ++){
    35         scanf("%d%d%d", &n, &m, &k);
    36         ans = 0;
    37         flag = (k & 1) ? 1 : -1;
    38         //计算g(k) 
    39         for(int i = 1; i <= k; i ++){
    40             ans = (ans + 1ll * flag * comb(k, i) * i % MOD * pow_mod(i-1, n-1, MOD) % MOD) % MOD;
    41             flag = -flag;
    42         }
    43         //接下来计算C(m, k) 
    44         temp = Finv[k];
    45         for(int i = 1; i <= k; i ++){
    46             temp = 1ll * temp * (m-k+i) % MOD;
    47         }
    48         ans = ((1ll * ans * temp) % MOD + MOD) % MOD;
    49         printf("Case #%d: %d
    ", cas, ans);
    50     }
    51 }
    View Code

    做完两题之后就感觉二项式反演变得容易了,遇到题目还是要多想( ̄▽ ̄")

    等等。。。做完两题的我突然发现二项式反演和容斥推倒出来的公式总是一样的。。。。。。(为什么有种被骗的感觉T_T)

  • 相关阅读:
    HDU3336 Count the string —— KMP next数组
    CodeForces
    51Nod 1627 瞬间移动 —— 组合数学
    51Nod 1158 全是1的最大子矩阵 —— 预处理 + 暴力枚举 or 单调栈
    51Nod 1225 余数之和 —— 分区枚举
    51Nod 1084 矩阵取数问题 V2 —— 最小费用最大流 or 多线程DP
    51Nod 机器人走方格 V3 —— 卡特兰数、Lucas定理
    51Nod XOR key —— 区间最大异或值 可持久化字典树
    HDU4825 Xor Sum —— Trie树
    51Nod 1515 明辨是非 —— 并查集 + 启发式合并
  • 原文地址:https://www.cnblogs.com/linyujun/p/5210650.html
Copyright © 2011-2022 走看看