結果

問題 No.147 試験監督（2）
ユーザー dnishdnish
提出日時 2018-07-02 10:40:18
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 511 ms / 2,000 ms
コード長 1,790 bytes
コンパイル時間 1,592 ms
コンパイル使用メモリ 173,724 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-09-13 16:51:57
合計ジャッジ時間 4,443 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 511 ms
4,376 KB
testcase_01 AC 509 ms
4,380 KB
testcase_02 AC 510 ms
4,376 KB
testcase_03 AC 1 ms
4,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include "bits/stdc++.h"
#define REP(i, n, N) for(ll i=(n); i<(N); i++)
#define RREP(i, n, N) for(ll i=(N-1); i>=n; i--)
#define CK(n, a, b) ((a)<=(n)&&(n)<(b))
#define ALL(v) (v).begin(),(v).end()
#define MCP(a,b) memcpy(b,a,sizeof(b))
#define p(s) cout<<(s)<<endl
#define p2(a, b) cout<<(a)<<" "<<(b)<<endl
#define v2(T) vector<vector<T>>
typedef long long ll;
using namespace std;
const ll mod = 1e9 + 7;
const ll inf = 1e18;
typedef vector<ll> vec;
typedef vector<vec> mat;

ll N,C;
//A*B
mat mul(mat &A, mat &B){
    mat C(A.size(), vec(B[0].size(),0));
    REP(i, 0, A.size())
        REP(k, 0, B.size())
            REP(j, 0, B[0].size())
                (C[i][j] += A[i][k] * B[k][j]) %= mod;
    return C;
}

//A^n
mat mat_pow(mat A, long long n){
    mat B(2, vec(2,0));
    REP(i, 0, 2) B[i][i] = 1;
    while(n > 0){
        if(n & 1) B = mul(B, A);
        A = mul(A, A);
        n >>= 1;
    }
    return B;
}

long long fast_power(long long a, long long b) {
    long long ret = 1;
    // aが負の場合
    if (a < 0) {
        a = -a;
        if (b % 2 == 1) ret = -1;
    }
    // 累乗計算
    while (b>0) {
        if (b & 1) {
            ret *= a;
        }
        a *= a;
        a %= mod;
        ret %= mod;
        b >>= 1;
    }
    return ret;
}

string D;
int main(){
    cin>>N;
    ll ans=1;
    mat A(2, vec(2));
    A[0][0] = A[1][0] = A[0][1] = 1;
    A[1][1] = 0;
    REP(i,0,N){
        cin>>C>>D;
        mat F = mat_pow(A, C+1);
        mat B(2, vec(1));
        B[0][0] = 1, B[1][0] = 0;
        F = mul(F, B);
        ll d=0;
        REP(i,0,D.size()){
            d *= 10;
            d += D[i] - '0';
            if(d>mod-1) d %= (mod-1);
        }
        ans *= fast_power(F[0][0] , d);
        ans %= mod;
    }
    p(ans);
    return 0;
}
0