FZU 1911 Construct a Matrix(矩阵快速幂+构造)

思路:S(n) = S(n-1) + f(n),通过构造矩阵可以求出S(n),主要问题就是构造一个S(n) * S(n)的矩阵满足条件。

自己写一写就能找出构造方法,因为矩阵只包括0,1,-1,那么一个n*n的矩阵,每一行的sum在-n到n之间,又因为-n和n不能同时存在,所以我们可以让每一行的和在-n到n-1之间,一共2*n个数,刚好对应n行,n列,每一行每一列的sum都不相同,然后构造就行了。

#include<cstdio>
#include<cstring>
using namespace std;
const int mod = 1e9 + 7;
const int maxn = 1e6 + 10;
typedef long long LL;
struct Mutrix{
    LL a[3][3];
    Mutrix(){
        memset(a, 0, sizeof a);
    }
};
Mutrix A, B;
LL m;
void init(){
    A.a[0][0] = 1, A.a[0][1] = 1, A.a[0][2] = 1;
    B.a[0][0] = 1;
    B.a[1][2] = 1;
    B.a[2][0] = 1, B.a[2][1] = 1, B.a[2][2] = 1;
}
Mutrix mul(Mutrix p1, Mutrix p2){
    Mutrix s;
    for(int i = 0; i < 3; ++i){
        for(int k = 0; k < 3; ++k){
            for(int j = 0; j < 3; ++j){
                s.a[i][j] += p1.a[i][k] * p2.a[k][j] % m;
                s.a[i][j] %= m;
            }
        }
    }
    return s;
}
Mutrix quick_mod(Mutrix x, LL y){
    Mutrix ret;
    for(int i = 0; i < 3; ++i) ret.a[i][i] = 1;
    while(y){
        if(y & 1) ret = mul(ret, x);
        x = mul(x, x);
        y >>= 1;
    }
    return ret;
}
int a[220][220];
int main(){
    LL n;
    int kase = 0, T;
    scanf("%d", &T);
    while(T--){
        scanf("%I64d%I64d", &n, &m);
        init();
        LL w = mul(A, quick_mod(B, n - 1)).a[0][0] % m;
        printf("Case %d: ", ++kase);
        if(w == 0 || w & 1) {
            printf("No\n");
            continue;
        }
        else printf("Yes\n");
        memset(a, 0, sizeof a);
        for(int i = 0; i < w; ++i){
            for(int j = 0; j < w; ++j){
                if(i < w / 2){
                    if(j >= i) a[i][j] = 1;
                }
                else{
                    if(j < i) a[i][j] = -1;
                }
            }
        }
        for(int i = 0; i < w; ++i){
            for(int j = 0; j < w; ++j){
                printf("%d ", a[i][j]);
            }
            printf("\n");
        }
    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值