結果

問題 No.991 N×Mマス計算(構築)
ユーザー betrue12betrue12
提出日時 2020-02-14 23:14:37
言語 C++17
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 134 ms / 2,000 ms
コード長 1,499 bytes
コンパイル時間 2,112 ms
コンパイル使用メモリ 198,216 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-04-16 02:52:28
合計ジャッジ時間 7,538 ms
ジャッジサーバーID
(参考情報)
judge4 / judge2
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
5,248 KB
testcase_01 AC 128 ms
5,376 KB
testcase_02 AC 128 ms
5,376 KB
testcase_03 AC 128 ms
5,376 KB
testcase_04 AC 128 ms
5,376 KB
testcase_05 AC 130 ms
5,376 KB
testcase_06 AC 130 ms
5,376 KB
testcase_07 AC 130 ms
5,376 KB
testcase_08 AC 132 ms
5,376 KB
testcase_09 AC 131 ms
5,376 KB
testcase_10 AC 129 ms
5,376 KB
testcase_11 AC 134 ms
5,376 KB
testcase_12 AC 131 ms
5,376 KB
testcase_13 AC 129 ms
5,376 KB
testcase_14 AC 130 ms
5,376 KB
testcase_15 AC 127 ms
5,376 KB
testcase_16 AC 130 ms
5,376 KB
testcase_17 AC 129 ms
5,376 KB
testcase_18 AC 132 ms
5,376 KB
testcase_19 AC 132 ms
5,376 KB
testcase_20 AC 128 ms
5,376 KB
testcase_21 AC 130 ms
5,376 KB
testcase_22 AC 130 ms
5,376 KB
testcase_23 AC 133 ms
5,376 KB
testcase_24 AC 131 ms
5,376 KB
testcase_25 AC 133 ms
5,376 KB
testcase_26 AC 129 ms
5,376 KB
testcase_27 AC 129 ms
5,376 KB
testcase_28 AC 134 ms
5,376 KB
testcase_29 AC 134 ms
5,376 KB
権限があれば一括ダウンロードができます

ソースコード

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