hdu 2065 红色病毒(指数型母函数)

求A B C D 在规定条件下n个元素的排列个数,先写出指数型母函数
G(X) = ( 1 + x + x^2/2! + x^3/3! +... )^2 * ( 1+ x + x^2/2! + x^4/! + .. )^2
前者表示：B, D出现方式不限制；后者表示：A, C 只能出现偶数或者不出现情况    
又知：
e^x=1+x/1!+x^2/2!+x^3/3!+...
e^(-x)=1-x/1!+x^2/2!-x^3/3!+...
化简得：
G(x)  = e^(2x) * ((e^x+e^(-x))/2)^2
 = (1/4) * e^(2x) * (e^(2x) + 2 + e^(-2x))
 = (1/4) * (e^(4x) + 2*e^(2x) +1)
　　 = (1/4) * ( (1+4x/1!+(4x)^2/2!+(4x)^3/3!+...+(4x)^n/n!) + 2*(1+2x/1!+(2x)^2/2!+(2x)^3/3!+...+(2x)^n/n!) +1)
 
得:  x^n 项系数 

a(n)  = (1/4) * ((4x)^n/n! + 2*(2x)^n/n!)
 = (1/4) * ( 4^n*x^n/n! + 2^(n+1)*x^n/n!)
 = (4^(n-1) + 2^(n-1)) * x^n/n!

即所求 F(n) = (4^(n-1) + 2^(n-1)) % 100

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int mod=100;
int fun(int a,long long m)
{
    int ans=1;
    while(m)
    {
        if(m&1) ans=(ans*a)%mod;
        a=a*a%mod;
        m>>=1;
    }
    return ans;
}
int main()
{
    int n,i,t;
    __int64 m;
    while(~scanf("%d",&n),n)
    {
        for(t=i=1;i<=n;i++)
        {
            scanf("%I64d",&m);
            printf("Case %d: %d\n",t++,(fun(2,m-1)+fun(4,m-1))%mod);
        }
        printf("\n");
    }
    return 0;
}