hdu 2065 "红色病毒"问题（快速幂）

最新推荐文章于 2021-01-26 00:05:38 发布

zthgreat

最新推荐文章于 2021-01-26 00:05:38 发布

阅读量647

点赞数

分类专栏：【基础题】文章标签：快速幂母函数

本文链接：https://blog.csdn.net/u014634338/article/details/38386525

版权

【基础题】专栏收录该内容

175 篇文章 1 订阅

订阅专栏

"红色病毒"问题

Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 3898 Accepted Submission(s): 1657

Problem Description

医学界发现的新病毒因其蔓延速度和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

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

Sample Input

Sample Output

  
  
   
   Case 1: 2
Case 2: 72
Case 3: 32
Case 4: 0

Case 1: 56
Case 2: 72
Case 3: 56

　证明(转)：从题目可以知道 A: (1 + x²/1! + x⁴/2! + ....); B: (1 + x/1! + x²/2! + x³/3! + ...); C:(1 + x²/1! + x⁴/2! + ....); D: (1 + x/1! + x²/2! + x³/3! + ...);

所以有： G(x) = (1 + x²/1! + x⁴/2! + ....)² * (1 + x/1! + x²/2! + x³/3! + ...)²;

由于泰勒展开式：

e^x = 1 + x/1! + x²/2! + x³/3! + ...

e^-x = 1 - x/1! + x²/2! - x³/3! + ...

所以

G(x) = e^2x + ((e^x + e^-x)/2)²;

= (1/4) * (e^2x + 1)²

= (1/4) * (e^4x + 2*e^2x + 1);

又因为：

e^4x = 1 + (4x)/1! + (4x)²/2! + (4x)³/3! + ... + (4x)ⁿ/n!;

e^2x = 1 + (2x)/1! + (2x)²/2! + (2x)³/3! + ... + (2x)ⁿ/n!;
所以：

n次幂的排列数为 (1/4)(4ⁿ + 2*2ⁿ)

即得所求的解为(1/4)(4ⁿ + 2*2ⁿ)%100 = (4^n-1 + 2^n-1)%100;

由组合数学公式 (a + b) % c = (a%c + b%c)%c;

所以可以用快速幂取模求解。

#include <iostream>
#include <cstring>
using namespace std;
long long quick_mod(long long a,long long b,int m)
{
    long long ans=1;
    while (b)
    {
        if(b&1)
        {
            ans=ans*a%m;
            b--;
        }
        b/=2;
        a=a*a%m;
    }
    return ans;
}
int main()
{
    int ans;
    long long n;
    int t;
    while(cin>>t && t)
    {
        ans=0;
        while (t--)
        {
            cin>>n;
            cout<<"Case "<<++ans<<": ";
            cout<<(quick_mod(4, n-1, 100)+quick_mod(2, n-1, 100))%100<<endl;
        }
        cout<<endl;
    }
    return 0;
}