P1349 广义斐波那契数列(矩阵的应用)

#include<iostream>
#include<cstdio>
#include<string.h>
#define mem(a,b) memset(a,b,sizeof(a));
using namespace std;
long  m;
long p,q,a1,a2,n;
struct mat
{
    long  a[2][2];
};
mat mat_mul(mat x,mat y)
{
    mat ans;
    mem(ans.a,0);
    for(int i=0;i<2;i++)
    {
        for(int j=0;j<2;j++)
        {
            for(int k=0;k<2;k++)
            {
                ans.a[i][j]+=x.a[i][k]*y.a[k][j];
            }
            ans.a[i][j]%=m;
        }
    }
    ret