HDU-2065 "红色病毒"问题（矩阵快速幂）

https://blog.csdn.net/idealism_xxm/article/details/51612355
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
对于每个测试实例,输出字符串个数的最后两位,每组输出后跟一个空行.
设dp[i][0~3]表示当前状态的方案数
i表示长度为i的字符串
0代表含有‘A’和’C’的个数均为偶数
1代表含有’A’的个数为奇数，含有‘C’的个数为偶数
2代表含有’C’的个数为奇数，含有‘A’的个数为偶数
3代表含有‘A’和’C’的个数均为奇数
然后枚举第i+1个字符，得到状态转移方程
dp[i+1][0]=dp[i][0]*2+dp[i][1]+dp[i][2];
dp[i+1][1]=dp[i][1]*2+dp[i][0]+dp[i][3];
dp[i+1][2]=dp[i][2]*2+dp[i][0]+dp[i][3];
dp[i+1][3]=dp[i][3]*2+dp[i][1]+dp[i][2];
则可得矩阵乘法

设其为：dp[i]*P=dp[i+1]
根据矩阵乘法的结合律可得：dp[0]*P^n=dp[n]
对P^n进行快速幂即可

import java.util.Scanner;
public class D2065 {
 static int[][] p = { { 2, 1, 1, 0 }, { 1, 2, 0, 1 }, { 1, 0, 2, 1 }, { 0, 1, 1, 2 } };
 static int[][] e = { { 1, 0, 0, 0 }, { 0, 1, 0, 0 }, { 0, 0, 1, 0 }, { 0, 0, 0, 1 } };
    static int  MOD=100;
 public static void main(String[] args) {
  Scanner input = new Scanner(System.in);
  while (input.hasNext()) {
   long n = input.nextLong();
   if(n==0)return;
   for (int i = 0; i < n; i++) {
    long x = input.nextLong();
                  System.out.println("Case "+(i+1)+": "+xx(x)[0][0]);
   }
              System.out.println();
  }
 }
static int[][] chen(int [][]a,int [][]b) {
 int [][]c=new int [4][4];
 for (int i = 0; i < c.length; i++) {
  for (int j = 0; j < c.length; j++) {
   c[i][j]=0;
   for (int k = 0; k < c.length; k++) {
    c[i][j]=(c[i][j]+a[i][k]*b[k][j])%MOD;
   }
  }
 }
 return c;
}
static int[][] xx(long x) {
   int [][]res=e;
   int [][]a=p;
   while(x>0) {
    if((x&1)==1) {
     res=chen(res,a);
    }
    a=chen(a,a);
    x>>=1; 
   }
    return res;
 }
}