結果

問題 No.991 N×Mマス計算(構築)
ユーザー betrue12
提出日時 2020-02-14 23:14:37
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 144 ms / 2,000 ms
コード長 1,499 bytes
コンパイル時間 2,118 ms
コンパイル使用メモリ 196,928 KB
最終ジャッジ日時 2025-01-09 00:30:51
ジャッジサーバーID
(参考情報)
judge4 / judge5
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 2
other AC * 28
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;

int64_t extgcd(int64_t a, int64_t b, int64_t& x, int64_t& y){
    int64_t d = a;
    if(b != 0){
        d = extgcd(b, a%b, y, x);
        y -= (a/b) * x;
    }else{
        x = 1; y = 0;
    }
    return d;
}

int64_t inv_mod(int64_t a, int64_t MOD){
    int64_t x, y;
    extgcd(a, MOD, x, y);
    return (MOD + x%MOD) % MOD;
}

int64_t gcd(int64_t a, int64_t b){
    return b==0 ? a : gcd(b, a%b);
}

int main(){
    int64_t X;
    cin >> X;
    if(X == 0){
        cout << "2 1 2" << endl;
        cout << "* 1" << endl;
        cout << "1\n1\n";
        return 0;
    }
    int64_t N = 1e5;
    int64_t Q = X/N, R = X%N;
    int64_t K = 1e9;
    for(int t=1; ; t++){
        if(gcd(K, 2+t) > 1) continue;
        int64_t v = (K/2*t%K*R + (2+t)*(N-R-1)) % K;
        v = (X-v+K) % K;
        if(v == 0) continue;
        int64_t x = inv_mod(2+t, K) * v % K;
        if(x < K/2){
            int M = Q + 1 + t;
            vector<int> A, B;
            for(int i=0; i<R; i++) A.push_back(K/2);
            A.push_back(x);
            while(A.size() < N) A.push_back(1);
            for(int i=0; i<Q; i++) B.push_back(K);
            B.push_back(2);
            while(B.size() < M) B.push_back(1);

            cout << N << " " << M << " " << K << endl;
            cout << "*";
            for(int b : B) cout << " " << b;
            cout << endl;
            for(int a : A) cout << a << endl;
            return 0;
        }
    }
    return 0;
}
0