洛谷 P1659【manacher】【快速幂】

#include <bits/stdc++.h>
#define int long long

using namespace std;

const int N = 1e6 + 10;
char s[N << 1];
int f[N << 1], mp[N << 1];

const int mod = 19930726;

int ksm(int a, int b, int c) {
  if (!b) return 1;
  int ans = ksm(a, b >> 1, c);
  ans = ans * ans % c;
  if (b & 1) ans = ans * a % c;
  return ans;
}

void manacher()
{
  int n, k;
  cin >> n >> k;
  scanf ("%s", s + 1);
  int l = strlen(s + 1);
  int nowmid = 1, nowr = 1;
  for (int i = 1; i <= l; i ++)
  {
    f[i] = 1;
    if (i <= nowr)
      f[i] = min(f[nowmid * 2 - i], nowr - i);
    while (i - f[i] >= 0 && i + f[i] <= l && s[i - f[i]] == s[i + f[i]])
      f[i] ++; // 进行暴力
    if (i + f[i] - 1 > nowr)
      nowmid = i, nowr = i + f[i];
    mp[f[i] * 2 - 1] ++;
  }
  int ans = 1, sum = 0;
  if (n % 2 == 0)
    n --;
  for (int i = n; i; i -= 2)
  {
    sum += mp[i];
    if (k < sum)
    {
      (ans *= ksm(i, k, mod)) %= mod;
      break;
    }
    (ans *= ksm(i, sum, mod)) %= mod; k -= sum;
  }
  cout << ans << '\n';
}

signed main()
{ manacher(); }

### 矩阵快速幂的概念 矩阵快速幂是一种高效的算法,用于计算一个给定矩阵的高次幂。该方法结合了矩阵乘法和快速幂的思想,在处理大指数运算时具有显著的时间复杂度优势[^1]。 ### C语言中的矩阵定义与初始化 为了实现矩阵快速幂,首先需要定义并初始化矩阵结构体: ```c #define N 105 // 假设最大维度不超过N*N typedef struct { long long mat[N][N]; } Matrix; ``` 此段代码声明了一个名为`Matrix`的结构体用来存储矩阵数据,并设定其大小上限为`N×N`。 ### 矩阵乘法规则 两个相同尺寸的方阵相乘遵循特定规则:新矩阵中每个元素等于原两矩阵对应行列交叉位置上的元素积之和。具体到编程层面可表示如下: ```c void multiply(Matrix* A, Matrix* B, Matrix* res) { static long long tmp[N][N]; memset(tmp, 0, sizeof(tmp)); for (int i = 0; i < N; ++i) for (int j = 0; j < N; ++j) for (int k = 0; k < N; ++k) tmp[i][j] += ((A->mat[i][k]) * (B->mat[k][j])); memcpy(res->mat, tmp, sizeof(tmp)); } ``` 这段程序实现了标准的三重循环形式的矩阵乘法操作,并将结果保存至目标矩阵对象内[^2]。 ### 快速幂逻辑应用于矩阵 接下来是核心部分——利用快速幂技巧加速矩阵自乘过程。这里给出完整的C语言版本实现方式: ```c // 定义单位矩阵I作为初始值 void initIdentityMatrix(Matrix* I){ memset(I->mat, 0, sizeof(I->mat)); for(int i=0;i<N;++i)(I->mat)[i][i]=1LL; } // 计算M^n的结果存入res中 void powerModulo(Matrix M, int n, Matrix* res){ Matrix result,tmp=M; initIdentityMatrix(&result); while(n>0){ if(n&1)multiply(&tmp,&result,&result); multiply(&tmp,&tmp,&tmp),n>>=1; } *res=result; } ``` 上述代码片段展示了如何通过迭代而非递归来完成整个矩阵快速幂的过程。其中特别注意的是当遇到奇数次幂时需额外做一次普通的矩阵乘法;而每次平方后的基底都要更新以便于下一轮判断是否继续累加当前项[^3]。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值