计算日期（快速幂+打表）

题意：
今天星期六，求1¹⁺²2……N^N天后是星期几

思路：
同余与模算术，利用快速幂取模的算法，时间复杂度为O(logn)。

1.先用快速幂求出1^1 ， 2^2 +，3^3 ， … ，N^N
对7取模之后的结果，发现循环节长度为42，即
(1^1)%7=(4343)%7,
(2^2)%7=(4444)%7,
(3^3)%7=(4545)%7,
(n^n)%7=( (42+n)^(42+n) )%7
2.然后打表求出[1,42]区间每个数n的(n^n)%7,再求a数组的前缀和b数组，

sum表示一个循环节所贡献的天数，即sum=（1^1 +3^3+ … 41^41 + 42^42）%7=6;
对于每一个样例n，直接计算即可

#include <math.h>
#include <stdio.h>
char chars[7][10] = {"Saturday",  "Sunday",   "Monday", "Tuesday",
                     "Wednesday", "Thursday", "Friday"};
// 快速乘  类似于快速幂
int qmul(int a, int b, int mod) {
  int ans = 0;
  while (b) {
    if (b & 1) {
      ans = (ans + a) % mod;
    }
    b >>= 1;
    a = (a << 1) % mod;
  }
  return ans;
}
// 实现快速幂
int qpow(int a, int n, int mod) {
  int res = 1;
  while (n) {
    if (n & 1) {
      // 原型为res = res*a %mod
      res = qmul(res, a, mod) % mod;
    }
    n >>= 1;
    // 原型为a = a * a % mod
    a = qmul(a, a, mod) % mod;
  }
  return res;
}

// a[n]表示n*n%7
int a[50];
// b[n]表示 a[1]+a[2]+a[3]...a[n]
int b[50];

int main() {
  b[0] = 0;
  for (int i = 1; i <= 42; i++) {
    a[i] = qpow(i, i, 7);
    b[i] = b[i - 1] + a[i];
  }
  int t;
  scanf("%d", &t);
  while (t--) {
    int n;
    scanf("%d", &n);
    int sum = b[42] % 7;
    // 一个循环总的天数
    // 循环节个数
    int num = n / 42;
    int ans = ((sum * num) % 7 + b[n % 42] % 7) % 7;
    printf("%s", chars[ans]);
  }
  return 0;
}

另一种求sum

for (int i = 1; i <= 42; i++) //循环节为42
    {
        a[i] =qpow(i, i, 7);
        b[i] = b[i - 1] + a[i];
        sum = (sum % 7 + a[i] % 7) % 7; //重点理解
    }