HDU 2065 "红色病毒"问题（指数生成函数+快速幂）

题目描述

传送门

题目大意：有A，B，C，D四种字母，问组成长度为N的字符串的方案数，其中A，C只能选偶数个（可不选，N<2^64）。

---

题解

指数型生成函数果果题。

不知道是谁说的：组合用常生成函数，排列用指数生成函数。

根据肥肥图的博客，答案就是$4^{n-1}+2^{n-1}$，快速幂一下。

将每个字母取法的生成函数相乘，得到一堆：

$$F(x)=(1+x+{x^2\over 2!}+{x^3\over 3!}+...)^2(1+{x^2\over2!}+{x^4\over4!}+{x^6\over6!}...)^2$$

而$e^x$的泰勒展开是$e^x=\sum_{i=0}^{+\infty}{x^i\over i!}$，$e^{-x}$就是$e^{-x}=\sum_{i=0}^{+\infty}{(-1)^i}{x^i\over i!}$，
所以$F(x)=e^{2x}({e^x+e^{-x}\over2})^2={e^{4x}+2e^{2x}+1\over4}$。

到这里直接求$n$阶导就得到第$n$项的系数：${4^n+2*2^n\over4}=4^{n-1}+2^{n-1}$。

也可以直接强拆，第$n$项
$A(n)={1\over4}·({(4x)^n\over n!}+{2*(2x)^n\over n!})=(4^{n-1}+2^{n-1} ){x^n\over n!}$。

又到了感性理解其正确性的时候了：对于取$n$个的排列方案，在考虑数量的前提下，乘上了在$n$个位置里选$i$个给一种放，另外$n-i$个位置则放另一种的方案数，即$A_i*{x^i\over i!}*B_{n-i}*{x^{n-i}\over{(n-i)!}}=A_i*B_{n-i}*{n!\over {i!*(n-i)!}}*{x^n\over n!}$，就是在卷积的时候乘上组合数。

换种形式表述的话，常生成函数就是在搞$C_m^n=\sum_{i=0}^nC_L^i*C_{m-L}^{n-i}$，而指数生成函数就是在搞$A_m^n=\sum_{i=0}^nC_n^i*A_L^i*A_{m-L}^{n-i}$。

所以这两种看似很神奇的生成函数，其本质也是基于排列组合。

---

代码

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <cmath>
#include <algorithm>
#define MOD 100

using namespace std;

typedef unsigned long long ULL;

int T;
ULL N;

int Pow(int x, ULL y){
    int res = 1;
    for(; y; x = x * x % MOD, y >>= 1)
        if(y & 1)  res = res * x % MOD;
    return res;
}

int main(){

    while(~ scanf("%d", &T) && T){

        for(int i = 1; i <= T; i++){
            scanf("%llu", &N);
            printf("Case %d: %d\n", i, (Pow(4, N-1) + Pow(2, N-1)) % MOD);
        }

        puts("");
    }

    return 0;
}

---