tjut 2842

#include <iostream>  
#include <cstring>  
#include <cstdio>  
#include <algorithm>  
#include <cmath>  
using namespace std;  
typedef long long ll;  
const int maxn = 10;  
const int mod = 200907;  
  
int cnt;  
struct Matrix {  
    int v[maxn][maxn];  
    Matrix() {}  
    Matrix(int x) {  
        init();  
        for (int i = 0; i < maxn; i++)   
            v[i][i] = x;  
    }  
    void init() {  
        memset(v, 0, sizeof(v));  
    }  
    Matrix operator *(Matrix const &b) const {  
        Matrix c;  
        c.init();  
        for (int i = 0; i < cnt; i++)  
            for (int j = 0; j < cnt; j++)  
                for (int k = 0; k < cnt; k++)  
                    c.v[i][j] = (c.v[i][j] + (ll)v[i][k]*b.v[k][j]) % mod;  
        return c;  
    }  
    Matrix operator ^(int b) {  
        Matrix a = *this, res(1);  
        while (b) {  
            if (b & 1)  
                res = res * a;  
            a = a * a;  
            b >>= 1;  
        }  
        return res;  
    }  
} a, b, tmp;  
  
int main() {  
    int n;  
    while (scanf("%d", &n) != EOF && n != 0) {  
        if (n < 3) {  
            printf("%d\n", n);  
            continue;  
        }  
        a.init();  
        cnt = 3;  
        a.v[0][0] = a.v[0][2] = a.v[1][0] = a.v[2][2] = 1;  
        a.v[0][1] = 2;  
        tmp = a^(n-2);  
        printf("%d\n", (tmp.v[0][0]*2 + tmp.v[0][1] + tmp.v[0][2]) % mod);  
    }  
    return 0;  
}  

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值