結果

問題 No.891 隣接3項間の漸化式
コンテスト
ユーザー pia019
提出日時 2026-10-08 16:07:55
言語 C++23
(gcc 15.3.0 + boost 1.92.0 + ACL)
コンパイル:
g++-15 -O2 -lm -std=c++23 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 1 ms / 2,000 ms
+ 382µs
コード長 1,901 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 5,319 ms
コンパイル使用メモリ 417,888 KB
実行使用メモリ 9,908 KB
最終ジャッジ日時 2026-10-08 16:08:05
合計ジャッジ時間 7,801 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge2_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 39
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
typedef long long ll;
const int INF = 1<<30;
const ll INFLL = 1LL<<60;
const ll MOD = 998244353;
const double INFD = 1.0E10;
const int dx[4] = {1, 0, -1, 0};
const int dy[4] = {0, -1, 0, 1};
map<int, char> dm = {{0, 'D'}, {1, 'L'}, {2, 'U'}, {3, 'R'}};
// const int dx[8] = {1, 1, 0, -1, -1, -1, 0, 1};
// const int dy[8] = {0, 1, 1, 1, 0, -1, -1, -1};
using Pair = pair<ll, ll>;
// using mint = atcoder::modint998244353;
using mint = atcoder::modint1000000007;

vector<mint> convolution_naive(const vector<mint>& a, const vector<mint>& b) {
    vector<mint> c(a.size() + b.size() - 1);
    for (int i = 0; i < (int)a.size(); i++) {
        for (int j = 0; j < (int)b.size(); j++) {
            c[i + j] += a[i] * b[j];
        }
    }

    return c;
}


mint bostan_mori(vector<mint> P, vector<mint> Q, long long n){
    while (n > 0){
        vector<mint> Q_neg(Q.size());
        for (int i = 0; i < (int)Q.size(); i++) Q_neg[i] = (i % 2) ? -Q[i] : Q[i];
        
        // vector<mint> U = atcoder::convolution(P, Q_neg);
        // vector<mint> V = atcoder::convolution(Q, Q_neg);
        vector<mint> U = convolution_naive(P, Q_neg);
        vector<mint> V = convolution_naive(Q, Q_neg);
        
        P.clear();
        Q.clear();
        
        for (int i = n & 1; i < (int)U.size(); i += 2) P.push_back(U[i]);
        for (int i = 0; i < (int)V.size(); i += 2) Q.push_back(V[i]);
        
        n >>= 1;
    }
    
    return P[0] / Q[0];
}


int main(){
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cout << fixed << setprecision(15);
    
    ll a, b, n; cin >> a >> b >> n;
    vector<mint> p = {0, 1};
    vector<mint> q = {1, -a, -b};
    
    mint ans = bostan_mori(p, q, n);
    
    cout << ans.val() << endl;
    
    return 0;
}
0