結果

問題 No.147 試験監督（2）
ユーザー kk
提出日時 2021-04-17 23:24:41
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 2,364 bytes
コンパイル時間 2,239 ms
コンパイル使用メモリ 201,492 KB
最終ジャッジ日時 2025-01-20 21:26:42
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 3 WA * 1
権限があれば一括ダウンロードができます

ソースコード

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;
    }

    ret *= modpow(fib(c), r, MOD);
    ret %= MOD;
  }

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