結果

問題 No.147 試験監督（2）
コンテスト
ユーザー dnish
提出日時 2018-07-02 10:26:11
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,778 bytes
コンパイル時間 1,689 ms
コンパイル使用メモリ 175,412 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-01 01:27:11
合計ジャッジ時間 4,338 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 3 WA * 1
権限があれば一括ダウンロードができます

ソースコード

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';
            d %= (mod-1);
        }
        ans *= fast_power(F[0][0] , d);
        ans %= mod;
    }
    p(ans);
    return 0;
}
0