結果

問題 No.2122 黄金比で擬似乱数生成
コンテスト
ユーザー 259_Momone
提出日時 2022-11-04 22:52:26
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 20 ms / 2,000 ms
コード長 1,160 bytes
コンパイル時間 2,729 ms
コンパイル使用メモリ 253,284 KB
最終ジャッジ日時 2025-02-08 18:14:52
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/extc++.h>

int main() {
    using namespace std;
    constexpr unsigned long mod{10000};
    unsigned long S, M, L;
    {
        string s;
        cin >> s;
        while(size(s) > 1 && s[0] == '0')s.erase(begin(s));
        S = stoi(s);
    }
    cin >> M >> L;
    const auto m_shuffle{[&M](unsigned long x){
        unsigned long a{}, b{1}, c{x}, d{1}, e{}, f{1};
        unsigned long n{M};
        while(n){
            if(n & 1){
                tie(d, e, f) = make_tuple(a * d + b * e, b * d + c * e, b * e + c * f);
                d %= mod;
                e %= mod;
                f %= mod;
            }
            tie(a, b, c) = make_tuple(a * a + b * b, b * (a + c), b * b + c * c);
            a %= mod;
            b %= mod;
            c %= mod;
            n /= 2;
        }
        if(M & 1)(e += 9999) %= mod;
        return e;
    }};
    vector<optional<unsigned long>> num(10000);
    for(unsigned long i{}; i < L; ++i){
        if(num[S]){
            const auto C{i - num[S].value()};
            i += (L - i) / C * C;
        }else num[S] = i;
        S = m_shuffle(S);
    }
    printf("%04lu\n", S);
    return 0;
}
0