FZU 1911 矩阵快速幂

题意:
给一个n,让你输出一个S(n)*S(n)的矩阵
S(n)是前n项斐波那契的和%m
这个矩阵必须满足每一行和每一列中任意两个的和不能相等
思路:
奇数肯定是无解的,因为奇数总有一行是相同的
偶数的情况,自己打个6*6和4*4的表就看出来规律了
就是把整个矩形分成三部分,左上部分和右下部分和从右上角到左下角的分界线
然后发现左上部分全是1,右下部分全是-1,分界线前n/2个是0,剩下是1就可以满足了
我们会发现每一列的和如果>0就是偶数,否则是奇数
每一列如果>0就是奇数,否则是偶数
另外:
我如果把LL n放到int main里就WA了,放到全局就AC了,这是什么鬼情况？
玄学？

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<math.h>
#include<queue>
#include<stack>
#include<string>
#include<vector>
#include<map>
#include<set>
using namespace std;
#define lowbit(x) (x&(-x))
typedef long long LL;
const int maxn = 100005;
const int inf=(1<<28)-1;
#define Matrix_Size 5
LL MOD;
int Size;
struct Matrix
{
    LL mat[Matrix_Size][Matrix_Size];
    void clear()
    {
        memset(mat,0,sizeof(mat));
    }
    void output()
    {
        for(int i = 0;i < Size;i++)
        {
            for(int j = 0;j < Size;j++)
                printf("%d ",mat[i][j]);
            printf("\n");
        }
    }
    Matrix operator *(const Matrix &b)const
    {
        Matrix ret;
        for(int i = 0;i < Size;i++)
            for(int j = 0;j < Size;j++)
            {
                ret.mat[i][j] = 0;
                for(int k = 0;k < Size;k++)
                {
                    long long tmp = (long long)mat[i][k]*b.mat[k][j]%MOD;
                    ret.mat[i][j] = (ret.mat[i][j]+tmp);
                    if(ret.mat[i][j]>=MOD)
                        ret.mat[i][j] -= MOD;
                    if(ret.mat[i][j]<0)//注意是否需要MOD 
                        ret.mat[i][j] += MOD;
                }
            }
        return ret;
    }
};
Matrix pow_M(Matrix a,long long n)
{
    Matrix ret;
    ret.clear();
    for(int i = 0;i < Size;i++)
        ret.mat[i][i] = 1;
    Matrix tmp = a;
    while(n)
    {
        if(n&1)ret = ret*tmp;
        tmp = tmp*tmp;
        n>>=1;
    }
    return ret;
}

LL get_row(LL n)
{
    Size=3;
    if(n==1) return 1%MOD;
    if(n==2) return 2%MOD;
    Matrix A,B;
    A.clear();B.clear();
    A.mat[0][0]=A.mat[0][1]=A.mat[0][2]=1;
    A.mat[1][1]=A.mat[1][2]=1;
    A.mat[2][1]=1;
    B.mat[0][0]=2;B.mat[1][0]=1;B.mat[2][0]=1;
    A=pow_M(A,n-2);
    A=A*B;
    return A.mat[0][0];
}


int Map[205][205];
LL n;
int main()
{
    int T,Case=0;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%lld%lld",&n,&MOD);
        printf("Case %d: ",++Case);
        n=get_row(n);
        //printf("****%lld\n",n);
        if(n&1||n==0)
        {
            printf("No\n");
            continue;
        }
        else
        printf("Yes\n");
        memset(Map,-1,sizeof(Map));
        for(int i=1;i<=n;++i)
        {
            for(int j=1;j<=n-i;++j)
                Map[i][j]=1;
            if(i<=n/2)
                Map[i][n-i+1]=0;
            else
                Map[i][n-i+1]=1;
        }
        for(int i=1;i<=n;++i)
        {
            for(int j=1;j<=n;++j)
            {   
                if(j>1) printf(" ");
                printf("%d",Map[i][j]);
            }
            printf("\n");
        }
    }
    return 0;
}