hdu 2065 "红色病毒"问题

最新推荐文章于 2020-03-08 18:32:47 发布

Tgotp

最新推荐文章于 2020-03-08 18:32:47 发布

阅读量439

点赞数

分类专栏： hdu 文章标签：生成函数

本文链接：https://blog.youkuaiyun.com/Tgotp/article/details/80085017

版权

hdu 专栏收录该内容

1 篇文章

订阅专栏

本文探讨了通过生成函数解析组合数列的问题，利用指数函数的幂级数展开形式进行推导，最终给出了一种计算特定数列通项公式的高效方法，并附带提供了C++实现代码。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

考虑生成函数
令 $G(x) = (1 + \frac{x}{1!} + \frac{x^2}{2!} + \frac{x^3}{3!} + \frac{x^4}{4!} \cdots)^2 * (1 + \frac{x^2}{2!} + \frac{x^4}{4!} \cdots)^2$
考虑 $e^x = 1 + \frac{x}{1!} + \frac{x^2}{2!} + \frac{x^3}{3!} + \frac{x^4}{4!} \cdots$
$e^{-x} = 1 - \frac{x}{1!} + \frac{x^2}{2!} - \frac{x^3}{3!} + \frac{x^4}{4!} \cdots$
则 $G(x) = e^{2x} + (\frac{e^x-e^{-x}}{2})^2$
化简可得 $\begin{align*} G(x) &=e^{2x}* \frac{e^{2x} + e^{-2x}+ 2}{4} \\ &= \frac{e^{4x} + 2*e^{2x}+ 1}{4}\\ &= \frac{(1 + \frac{4*x}{1!} + \frac{16*x^2}{2!} + \frac{64*x^3}{3!} + \frac{256*x^4}{4!} \cdots) + 2* (1 + \frac{2*x}{1!} + \frac{4*x^2}{2!} + \frac{8*x^3}{3!} + \frac{16*x^4}{4!} \cdots) + 1}{4}\\ &= (\frac{1}{4} + \frac{x}{1!} + \frac{4*x^2}{2!} + \frac{16*x^3}{3!} + \frac{64*x^4}{4!} \cdots) + (\frac{1}{2} + \frac{x}{1!} + \frac{2*x^2}{2!} + \frac{4*x^3}{3!} + \frac{8*x^4}{4!} \cdots) + \frac{1}{4} \end{align*}$
所以综上第 $n$ 项答案为 $4^{n-1}+2^{n-1} (mod\ 100)$
c++代码如下：

#include<bits/stdc++.h>
#define rep(i,x,y) for(register int i = x ; i <= y; ++ i)
#define repd(i,x,y) for(register int i = x ; i >= y; -- i)
using namespace std;
typedef long long ll;
template<typename T>inline void read(T&x)
{
    char c;int sign = 1;x = 0;
    do { c = getchar(); if( c == '-' ) sign = -1; }while(!isdigit(c));
    do { x = x * 10 + c - '0'; c = getchar(); }while(isdigit(c));
    x *= sign;
}

const int mod = 100;
int ksm(int x,ll y)
{
    int ans = 1;
    while(y)
    {
        if(y&1) ans = ans * x %mod;
        x = x*x %mod;
        y >>= 1;
    }
    return ans;
}

int main()
{
    int t; ll n;

    while(true){
    read(t);
    if(t == 0) break;
    rep(i,1,t)
    {
        read(n);
        printf("Case %d: %d\n",i,(ksm(2,n-1)+ksm(4,n-1))%mod );
    }
    puts("");
    }
    return 0;    
}