結果

問題 No.147 試験監督（2）
ユーザー HachimoriHachimori
提出日時 2015-02-09 01:31:51
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 935 ms / 2,000 ms
コード長 2,267 bytes
コンパイル時間 573 ms
コンパイル使用メモリ 58,188 KB
実行使用メモリ 27,776 KB
最終ジャッジ日時 2024-06-23 11:22:17
合計ジャッジ時間 5,039 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 935 ms
27,724 KB
testcase_01 AC 912 ms
27,776 KB
testcase_02 AC 919 ms
27,776 KB
testcase_03 AC 2 ms
6,940 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
const int NUM_DATA = 30005;
const int MOD = 1000000007;
const int LENG = 205;

class Matrix {
public:
    int v[2][2];

    Matrix(){
        memset(v, 0, sizeof(v));
    }
};


int N;
long long C[NUM_DATA];
int lengD[NUM_DATA], D[NUM_DATA][LENG];

void read() {
    cin >> N;
    for (int i = 0; i < N; ++i) {
        string DStr;
        cin >> C[i] >> DStr;
        lengD[i] = (DStr.size() + 3) / 4;
        for (int j = DStr.size() - 1, len = lengD[i] - 1; j >= 0; --len, j -= 4) {
            for (int k = max(j - 3, 0); k <= j; ++k)
                D[i][len] = D[i][len] * 10 + DStr[k] - '0';
        }
    }
}


Matrix matmul(const Matrix &a, const Matrix &b) {
    Matrix c;
    for (int i = 0; i < 2; ++i)
        for (int j = 0; j < 2; ++j) {
            int toPush = 0;
            for (int k = 0; k < 2; ++k)
                toPush = (toPush + 1LL * a.v[i][k] * b.v[k][j]) % MOD;
            c.v[i][j] = toPush;
        }
    return c;
}


Matrix matpow(Matrix &p, long long n) {
    if (n == 0) {
        Matrix I;
        I.v[0][0] = 1; I.v[0][1] = 0;
        I.v[1][0] = 0; I.v[1][1] = 1;
        return I;
    }
    
    Matrix t = matpow(p, n / 2);
    return n & 1 ? matmul(matmul(t, t), p) : matmul(t, t);
}


int mypow(int p, int leng, int n[LENG]) {
    int bgn = 0;
    int ret = 1;
    
    while (1) {
        if (n[leng - 1] % 2 == 1) {
            ret = 1LL * ret * p % MOD;
        }
        p = 1LL * p * p % MOD;
        
        
        while (bgn < leng && n[bgn] == 0) ++bgn;
        if (bgn == leng) break;
        
        // n /= 2
        int carry = 0;
        for (int i = bgn; i < leng; ++i) {
            int nextCarry = (n[i] + carry * 10000) & 1;
            n[i] = (n[i] + carry * 10000) >> 1;
            carry = nextCarry;
        }
    }
    
    return ret;
}


void work() {
    // 1, 1
    // 1, 0
    Matrix mat;
    mat.v[0][0] = 1; mat.v[0][1] = 1;
    mat.v[1][0] = 1; mat.v[1][1] = 0;
    
    int ans = 1;
    for (int i = 0; i < N; ++i) {
        Matrix t = matpow(mat, C[i] + 1);
        ans = 1LL * ans * mypow(t.v[0][0], lengD[i], D[i]) % MOD;
    }
    cout << ans << endl;
}


int main() {
    read();
    work();
    return 0;
}
0