循环矩阵

#include <stdio.h>
#include <string.h>
#define ll long long
const int maxn = 505;
ll n,m,d,k,a[maxn];
struct mat{ // 结构体循环矩阵，原矩阵式方阵
    ll v[maxn];
    mat() {memset(v,0,sizeof(v));} // 构造函数
    mat operator*(mat c){ // 重载乘法运算符
        mat ans;
        // 只需计算出循环矩阵的第一行即可
        for(int i = 0 ; i < n ; i++)
        for(int j = 0 ; j < n ; j++)
            ans.v[i] = ((ans.v[i] + v[j] * c.v[(i - j + n) % n]) % m + m) % m;
          return ans;
     }
};
// 此处不用担心溢出，况且考虑到不好写单位矩阵，就用的方式啦
mat pow_mod(mat x,ll k){
    mat y; y.v[0] = 1;
    while(k){
        if(k & 1) y = y * x;
        x = x * x;
        k >>= 1;
    }
    return y;
    // 递归写法
    //if(k == 1) return x;
    //mat sb = x * x;
    //mat ans = pow_mod(sb,k >> 1);
    //if(k & 1) ans = ans * x;
    //return ans;
}