【组合数学】 红色病毒问题 详解（指数型生成函数+泰勒展开）

红色病毒问题：

题目链接：https://vjudge.net/problem/HDU-2065

医学界发现的新病毒因其蔓延速度和Internet上传播的"红色病毒"不相上下,被称为"红色病毒",经研究发现,该病毒及其变种的DNA的一条单链中,胞嘧啶,腺嘧啶均是成对出现的。
现在有一长度为N的字符串,满足一下条件:
(1) 字符串仅由A,B,C,D四个字母组成;
(2) A出现偶数次(也可以不出现);
(3) C出现偶数次(也可以不出现);
计算满足条件的字符串个数.
当N=2时,所有满足条件的字符串有如下6个:BB,BD,DB,DD,AA,CC.
由于这个数据肯能非常庞大,你只要给出最后两位数字即可.

Input

每组输入的第一行是一个整数T,表示测试实例的个数,下面是T行数据,每行一个整数N(1<=N<2^64),当T=0时结束.

Output

对于每个测试实例,输出字符串个数的最后两位,每组输出后跟一个空行.

分析：

转换为多重集排列问题：

有四种物品，其中B，D能取1~N个，A，C能取1~N中的偶数个，求其排列个数。

 指数型母函数解多重集排列问题原理：指数型母函数

 构造母函数：

 

前一个多项式代表B和D，后一个多项式代表A和C，因为只能取偶数个，所以指数只包含偶数。

但是这题N太大了，没法直接模拟多项式相乘来做，所以要借助泰勒展开：

泰勒展开：

 

也就是说e^x近似等于上述的式子

那我们就将第一个多项式化为了

而第二个多项式就是去掉了第一个多项式去掉了奇数项

因此第二个多项式可以化为

我们的构造的母函数就化为了

而我们要求的是多项式中的系数

我们将和再次展开为多项式

取的系数即

快速幂即可

所以这题是考数学吧（摔桌！！！）

AC代码：

#include<iostream>
#include<cstdio>
#include<string.h>
using namespace std;

long long quickpower(long long base,long long power){
   long long ret=1;
   while(power){
     if(power%2)
        ret=ret*base%100;
     base=base*base%100;
     power/=2;
   }
   return ret;
}

int main()
{
   int T;
   while(~scanf("%d",&T)&&T){
      for(int i=1;i<=T;i++){
          long long N;
          scanf("%lld",&N);
          printf("Case %d: %lld\n",i,(quickpower(4,N-1)+quickpower(2,N-1))%100);
      }
      printf("\n");
   }

   return 0;
}