hdu 2065 "红色病毒"问题(快速幂)

"红色病毒"问题

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
  
  
4 1 4 20 11 3 14 24 6 0
 

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 + x2/1! + x4/2! + ....); B:  (1 + x/1! + x2/2! + x3/3! + ...); C:(1 + x2/1! + x4/2! + ....); D: (1 + x/1! + x2/2! + x3/3! + ...);

所以有: G(x) = (1 + x2/1! + x4/2! + ....)2 * (1 + x/1! + x2/2! + x3/3! + ...)2;

由于泰勒展开式:

ex = 1 + x/1! + x2/2! + x3/3! + ...

e-x = 1 - x/1! + x2/2! - x3/3! + ...

所以

G(x) = e2x + ((ex + e-x)/2)2;

 = (1/4) * (e2x + 1)2

 = (1/4) * (e4x  + 2*e2x + 1);

又因为:

e4x = 1 + (4x)/1! + (4x)2/2! + (4x)3/3! + ... + (4x)n/n!;

e2x = 1 + (2x)/1! + (2x)2/2! + (2x)3/3! + ... + (2x)n/n!;
所以:

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

即得所求的解为(1/4)(4n + 2*2n)%100 = (4n-1 + 2n-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;
}



  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值