結果

問題 No.147 試験監督（2）
コンテスト
ユーザー tnakao0123
提出日時 2016-03-22 17:53:23
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 185 ms / 2,000 ms
コード長 1,872 bytes
コンパイル時間 882 ms
コンパイル使用メモリ 86,720 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-01 13:02:25
合計ジャッジ時間 2,434 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 4
権限があれば一括ダウンロードができます

ソースコード

diff # 

/* -*- coding: utf-8 -*-
 *
 * 147.cc: No.147 試験監督(2) - yukicoder
 */

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<iostream>
#include<string>
#include<vector>
#include<map>
#include<set>
#include<stack>
#include<list>
#include<queue>
#include<deque>
#include<algorithm>
#include<numeric>
#include<utility>
#include<complex>
#include<functional>
 
using namespace std;

/* constant */

typedef long long ll;
const ll MOD = 1000000007;

/* typedef */

typedef ll mat[2][2];

/* global variables */

/* subroutines */

void mulmat(const mat &a, const mat &b, mat &c) { // a * b -> c
  c[0][0] = (a[0][0] * b[0][0] + a[0][1] * b[1][0]) % MOD;
  c[0][1] = (a[0][0] * b[0][1] + a[0][1] * b[1][1]) % MOD;
  c[1][0] = (a[1][0] * b[0][0] + a[1][1] * b[1][0]) % MOD;
  c[1][1] = (a[1][0] * b[0][1] + a[1][1] * b[1][1]) % MOD;
}

void powmat(const mat &a, ll b, mat &c) { // a^b -> c
  c[0][0] = c[1][1] = 1, c[0][1] = c[1][0] = 0;
  mat t, a0;
  memcpy(a0, a, sizeof(a));
  
  while (b > 0) {
    if ((b & 1) != 0) {
      mulmat(a0, c, t);
      memcpy(c, t, sizeof(t));
    }

    mulmat(a0, a0, t);
    memcpy(a0, t, sizeof(t));
    b >>= 1;
  }
}

ll powmod(ll a, ll n) {  // a^n % MOD
  ll pm = 1;
  while (n > 0) {
    if ((n & 1) != 0) pm = (pm * a) % MOD;
    a = (a * a) % MOD;
    n >>= 1;
  }
  return pm;
}

/* main */

int main() {
  int n;
  cin >> n;

  ll mul = 1;

  while (n--) {
    ll ci;
    string ds;
    cin >> ci >> ds;

    ll di = 0;
    int dn = ds.size();
    for (int i = 0; i < dn; i++) di = (di * 10 + (ds[i] - '0')) % (MOD - 1);

    mat a = {{0, 1}, {1, 1}}, c;
    powmat(a, ci, c);

    ll k = (c[0][0] * 1 + c[0][1] * 2) % MOD;
    ll kd = (k == 0) ? 0 : powmod(k, di);
    mul = (mul * kd) % MOD;

    //printf("%lld %lld -> %lld -> %lld\n", ci, di, k, kd);
  }

  printf("%lld\n", mul);
  return 0;
}
0