結果

問題 No.147 試験監督(2)
ユーザー kk
提出日時 2021-04-17 23:28:56
言語 C++17(1z)
(gcc 10.2.0 + boost 1.73.0)
結果
AC  
実行時間 599 ms / 2,000 ms
コード長 2,400 Byte
コンパイル時間 2,068 ms
使用メモリ 7,028 KB
最終ジャッジ日時 2021-04-17 23:29:01
合計ジャッジ時間 5,553 ms
ジャッジサーバーID
(参考情報)
judge10 / judge7
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
使用メモリ
testcase_00 AC 599 ms
3,352 KB
testcase_01 AC 582 ms
3,316 KB
testcase_02 AC 578 ms
3,312 KB
testcase_03 AC 1 ms
3,332 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>

using namespace std;

const long long MOD = 1e9 + 7;

typedef vector<long long> vec;
typedef vector<vec> matrix;

matrix eye(int n) {
  matrix mat(n, vec(n));
  for (int i = 0; i < n; i++) {
    mat[i][i] = 1;
  }
  return mat;
}

matrix zeros(int r, int c) {
  return matrix(r, vec(c));
}

matrix zeros(int n) {
  return zeros(n, n);
}

matrix add(const matrix &a, const matrix &b) {
  matrix c = zeros(a.size(), a[0].size());
  for (int i = 0; i < c.size(); i++) {
    for (int j = 0; j < c[0].size(); j++) {
      c[i][j] = a[i][j] + b[i][j];
    }
  }
  return c;
}

matrix mul(const matrix &a, const matrix &b) {
  matrix c = zeros(a.size(), b[0].size());
  for (int i = 0; i < a.size(); i++) {
    for (int k = 0; k < b.size(); k++) {
      for (int j = 0; j < b[0].size(); j++) {
        c[i][j] += a[i][k] * b[k][j];
      }
    }
  }
  return c;
}

matrix modmul(const matrix &a, const matrix &b, long long mod) {
  matrix c = zeros(a.size(), b[0].size());
  for (int i = 0; i < a.size(); i++) {
    for (int k = 0; k < b.size(); k++) {
      for (int j = 0; j < b[0].size(); j++) {
        c[i][j] = (c[i][j] + a[i][k] * b[k][j]) % mod;
      }
    }
  }
  return c;
}

matrix pow(const matrix &a, long long p) {
  matrix ret = eye(a.size());
  matrix x = a;
  while (p) {
    if (p & 1)
      ret = mul(ret, x);
    x = mul(x, x);
    p >>= 1;
  }
  return ret;
}

matrix modpow(const matrix &a, long long p, long long mod) {
  matrix ret = eye(a.size());
  matrix x = a;
  while (p) {
    if (p & 1)
      ret = modmul(ret, x, mod);
    x = modmul(x, x, mod);
    p >>= 1;
  }
  return ret;
}

long long fib(long long n) {
  matrix a = zeros(2);
  a[0][0] = a[0][1] = a[1][0] = 1;
  matrix b = modpow(a, n, MOD);
  return (b[0][0] + b[0][1]) % MOD;
}

long long modpow(long long x, long long p, long long mod) {
  long long ret = 1;
  while (p) {
    if (p & 1)
      ret = ret * x % mod;
    x = x * x % mod;
    p >>= 1;
  }
  return ret;
}

int main() {
  ios_base::sync_with_stdio(0);
  cin.tie(0);

  int n;
  cin >> n;

  long long ret = 1;
  for (int i = 0; i < n; i++) {
    long long c;
    string d;
    cin >> c >> d;

    long long r = 0;
    for (char x: d) {
      r = r * 10 + (x - '0');
      r %= MOD - 1;
    }

    if (!r)
      r = MOD - 1;
    
    ret *= modpow(fib(c), r, MOD);
    ret %= MOD;
  }

  cout << ret << endl;
  
  return 0;
}
0