中国石油大学ACM俱乐部开放训练赛 问题 F: 求和 等比矩阵求和

等比数列是指从第二项起，每一项与它的前一项的比值等于同一个常数的一种数列。对于一个等比数列an=a1 * q^n-1，它的前n项的和Sn=a1 * (1-qⁿ)/(1-q)(q≠1)。现在已知A为n*n的矩阵，S=A+A²+A³+…+A^m，你能否正确求出S，并且输出S中的每一个元素对1000000007取模后的值。

输入

输入包括n+1行，第一行包括两个正整数n, m，分别代表矩阵A的大小和S中的项数，其中1≤n≤30, 1≤m≤10⁹。接下来n行，每行n个元素，相应地代表A中的元素x，其中0≤x≤10⁶。

输出

输出包括n行，每行n个元素，相应地代表S中的每一个元素对1000000007取模后的值。

样例输入 Copy

1 2019
1

样例输出 Copy

2019

第一反应就是分治，因为等比可分治。
纪念lz学长，带老。
不过题解不是分治写的。
这里就都记一下啦，分治公式以后用到了自己推。。

分治写法：

#include<iostream>
#include<cstring>
#include<string>
#include<algorithm>
#include<cstdio>
#include<cmath>
#include<map>
#include<queue>
#define ll long long
using namespace std;
const int maxn=35;
const int mod=1000000007;
struct node
{
   
    ll a[maxn][maxn];
}res;
 
int n,m;
 
node operator * (const node &a,const node &b)
{
   
    node c;
    memset(c.a,0,sizeof