矩(方)阵运算

打完一个矩阵快速幂的板子以后心血来潮写了个namespace……不过因为是快速幂的板子所以也就能支持长宽一致的矩阵之间的运算，别的之后再改吧…
这篇Blog可能随时就会改掉，毕竟之后还要往里面塞一些乱七八糟的东西什么的，比如至少先支持真正的矩阵之间的加减乘运算的说…
列个表格大致说一下变量的作用吧。

---

变量名称	变量作用
n	矩(方)阵的长&&宽
q	进行的运算代号,1为加,2为减,3为对一个矩阵的快速幂(之后总会加别的运算的)
k	幂运算的指数 (仅幂运算)
a	输入的矩阵1
b	输入的矩阵2 (仅加减运算)
NULL	预留位置

---

为了增 jiang 强 di 代码可读性，就选择了重载运算符的方法瞎搞这个程序……到底也是搞出来了还扔到了namespace里面……不说别的了上代码:

#include<cstdio>
#define N 101
#define mod 1000000007
#define int long long
int n,k,q;
namespace Square{
    struct square{
        int n[N][N];
    };
    void reset(square &s){
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                s.n[i][j]=0;
    }
    square operator +(square &a,square &b){
        square c;reset(c);
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                c.n[i][j]=a.n[i][j]+b.n[i][j];
        return c;
    }
    square operator -(square &a,square &b){
        square c;reset(c);
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                c.n[i][j]=a.n[i][j]-b.n[i][j];
        return c;
    }
    square operator *(square &a,square &b){
        square c;reset(c);
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                for(int k=1;k<=n;k++)
                    c.n[i][j]=(c.n[i][j]+a.n[i][k]*b.n[k][j]%mod)%mod;
        return c;
    }
    square operator ^(square a,int k){
        square rst=a,r=a;k--;
        while(k){
            if(k&1)rst=rst*r;
            r=r*r;k>>=1;
        }return rst;
    }
}using namespace Square;
square a,b,ans;
main(){
    scanf("%lld%lld",&n,&q);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            scanf("%lld",&a.n[i][j]),a.n[i][j]%=mod;
    if(q!=3)
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                scanf("%lld",&b.n[i][j]),b.n[i][j]%=mod;
    else scanf("%lld",&k);
    switch(q){
        case 1:ans=a+b;break;
        case 2:ans=a-b;break;
        case 3:ans=a^k;break;
    }
    for(int i=1;i<=n;i++,putchar('\n'))
        for(int j=1;j<=n;j++,putchar(' '))
            printf("%lld",ans.n[i][j]);
}

至于为什么把int给define成long long？嗯……那个模板题数据范围比较鬼，而我比较懒不想一个一个改，就那么搞呗……
希望这不是一个弃坑的Blog就好了