指数型母函数&&hdu 2065 "红色病毒"问题

最新推荐文章于 2020-11-04 21:29:48 发布

theArcticOcean

最新推荐文章于 2020-11-04 21:29:48 发布

阅读量843

点赞数

分类专栏： algorithm_组合数学文章标签： algorithm

本文链接：https://blog.csdn.net/theArcticOcean/article/details/46913347

版权

algorithm_组合数学专栏收录该内容

42 篇文章 0 订阅

订阅专栏

常用的指数型母函数：e^x=1+x/1!+x^2/2!+x^3/3!+…… e^(-x)=1-x/1!+x^2/2!-x^3/3!+……
指数型母函数用于解决多重集合的排列问题，g(x)=p(n,0)+p(n,1)x+p(n,2)x^2/2!+p(n,3)x^3/3!+……+P(n,n)x^n/n!=C(n,0)+C(n,1)x+C(n,2)x^2+C(n,3)x^3+……+C(n,n)x^n=(1+x)^n.所以(1+x)^n是n的各个排列的母函数。比如用red,white,black给1*n的棋盘涂色，方案数：g(x)=(1+x+x^2/2!+x^3/3!+……)(1+x+x^2/2!+x^3/3!+……)(1+x+x^2/2/3!+……)=e^(3x)=∑(0,∞)3^n*x^n/n!.-->h(n)=3^n。
如果要求red出现偶数次，总的方案数应该是：g(x)=(1+x+x^2/2!+x^4/4!+……)(1+x+x^2/2!+x^3/3!+……)(1+x+x^2/2!+x^3/3!+……)=1/2*(e^x+e^(-x))*e^x*e^(-x)=1/2*(e^3x+e^x)=1/2*(∑(0,∞)(3^n+1)*x^n/n!. --> h(n)=(3^n+1)/2.
hdu 2065 "新红色病毒"问题: http://acm.hdu.edu.cn/showproblem.php?pid=2065
A,C可以出现偶数次，B,D可以是任意次数。g(x)=(1+x+x^2/2!+x^4/4!+……)^2*(1+x+x^2/2!+x^3/3!+……)^2=[1/2*(e^x+e^-x)]^2*[e^x]^2=[1/2*(e^2x+1)]^2=1/4*(e^4x+2e^2x+1)=1/4*[ ∑(0,∞)4^n*x^n/n!+2*∑(0,∞)2^n*x^n/n!+1]=∑(0,∞)4^(n-1)*x^n/n!+∑(0,∞)2^(n-1)*x^n/n!+1/4 h(n)=4^(n-1)+2^(n-1). 最后对这个结果取模即可(用到快速幂)。

#include <iostream>
#include<cstdio>
using namespace std;
typedef long long LL;
LL powmod(LL n){
    LL q1=n-1,q2=n-1,res1=1,t1=4,res2=1,t2=2;
    while(q1){
        if(q1&1){
            res1=(res1*t1)%100;
            q1--;
        }
        q1>>=1;
        t1=t1*t1%100;
    }
    while(q2){
        if(q2&1){
            res2=(res2*t2)%100;
            q2--;
        }
        q2>>=1;
        t2=t2*t2%100;
    }
    return (res1+res2)%100;
}
int main()
{
    //freopen("cin.txt","r",stdin);
    int t;
    while(cin>>t&&t){
        LL n;
        for(int i=1;i<=t;i++){
            scanf("%I64d",&n);
            printf("Case %d: %I64d\n",i,powmod(n));
        }
        putchar('\n');
    }
    return 0;
}