結果

問題 No.147 試験監督(2)
ユーザー kk
提出日時 2021-04-17 23:15:58
言語 C++17(1z)
(gcc 10.2.0 + boost 1.73.0)
結果
AC  
実行時間 1,928 ms / 2,000 ms
コード長 2,410 Byte
コンパイル時間 2,240 ms
使用メモリ 3,448 KB
最終ジャッジ日時 2021-04-17 23:16:09
合計ジャッジ時間 11,497 ms
ジャッジサーバーID
(参考情報)
judge8 / judge7
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
使用メモリ
testcase_00 AC 1,928 ms
3,420 KB
testcase_01 AC 1,928 ms
3,428 KB
testcase_02 AC 1,926 ms
3,380 KB
testcase_03 AC 1 ms
3,348 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, string p, long long MOD) {
  reverse(p.begin(), p.end());

  long long ret = 1;
  long long base = x;
  for (char c: p) {
    for (int i = 0; i < (c - '0'); i++)
      ret = ret * base % MOD;
    long long t = 1;
    for (int i = 0; i < 10; i++)
      t = t * base % MOD;
    base = t;
  }
  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;
    ret *= modpow(fib(c), d, MOD);
    ret %= MOD;
  }

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