結果

問題 No.2405 Minimal Matrix Decomposition
ユーザー 👑 獅子座じゃない人獅子座じゃない人
提出日時 2023-07-28 13:42:09
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
(最新)
AC  
(最初)
実行時間 -
コード長 2,568 bytes
コンパイル時間 1,825 ms
コンパイル使用メモリ 177,004 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-08 04:53:34
合計ジャッジ時間 6,875 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,248 KB
testcase_02 AC 2 ms
5,248 KB
testcase_03 WA -
testcase_04 AC 41 ms
5,248 KB
testcase_05 AC 170 ms
5,248 KB
testcase_06 AC 2 ms
5,248 KB
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 AC 84 ms
5,248 KB
testcase_11 AC 107 ms
5,248 KB
testcase_12 AC 4 ms
5,248 KB
testcase_13 AC 29 ms
5,248 KB
testcase_14 AC 26 ms
5,248 KB
testcase_15 AC 11 ms
5,248 KB
testcase_16 AC 11 ms
5,248 KB
testcase_17 AC 48 ms
5,248 KB
testcase_18 AC 29 ms
5,248 KB
testcase_19 AC 28 ms
5,248 KB
testcase_20 AC 5 ms
5,248 KB
testcase_21 AC 8 ms
5,248 KB
testcase_22 AC 87 ms
5,248 KB
testcase_23 AC 3 ms
5,248 KB
testcase_24 AC 4 ms
5,248 KB
testcase_25 AC 6 ms
5,248 KB
testcase_26 AC 61 ms
5,248 KB
testcase_27 AC 23 ms
5,248 KB
testcase_28 AC 51 ms
5,248 KB
testcase_29 AC 16 ms
5,248 KB
testcase_30 AC 71 ms
5,248 KB
testcase_31 AC 32 ms
5,248 KB
testcase_32 AC 8 ms
5,248 KB
testcase_33 AC 57 ms
5,248 KB
testcase_34 AC 16 ms
5,248 KB
testcase_35 AC 2 ms
5,248 KB
testcase_36 AC 28 ms
5,248 KB
testcase_37 AC 38 ms
5,248 KB
testcase_38 AC 37 ms
5,248 KB
testcase_39 AC 61 ms
5,248 KB
testcase_40 AC 37 ms
5,248 KB
testcase_41 AC 6 ms
5,248 KB
testcase_42 AC 77 ms
5,248 KB
testcase_43 AC 23 ms
5,248 KB
testcase_44 AC 18 ms
5,248 KB
testcase_45 AC 8 ms
5,248 KB
testcase_46 AC 16 ms
5,248 KB
testcase_47 AC 20 ms
5,248 KB
testcase_48 AC 9 ms
5,248 KB
testcase_49 AC 4 ms
5,248 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

#include <atcoder/modint>
using namespace atcoder;
using mint=modint;

int main(void)
{
    int p;
    cin >> p;
    int n,m;
    cin >> n >> m;
    mint::set_mod(p);
    vector<vector<mint>> a(n,vector<mint>(m));
    for(int i=0;i<n;++i){
        for(int j=0;j<m;++j){
            int atmp;
            cin >> atmp;
            a[i][j]=atmp;
        }
    }
    vector<vector<mint>> a_copy=a;
    int rank=0;
    for(int i=0;i<m;++i){
        if(rank==n){
            break;
        }
        if(a_copy[rank][i]==0){
            bool not_pivot=true;
            for(int j=rank+1;j<n;++j){
                if(a_copy[j][i]!=0){
                    not_pivot=false;
                    for(int k=i;k<m;++k){
                        swap(a_copy[rank][k],a_copy[j][k]);
                    }
                    break;
                }
            }
            if(not_pivot){
                continue;
            }
        }
        for(int j=m-1;j>=i;--j){
            a_copy[rank][j]/=a_copy[rank][i];
        }
        for(int j=0;j<n;++j){
            if(j==rank){
                continue;
            }
            for(int k=m-1;k>=i;--k){
                a_copy[j][k]-=a_copy[rank][k]*a_copy[j][i];
            }
        }
        ++rank;
    }
    if(n*m<=n*rank+rank*m){
        cout << 1 << endl;
        cout << n << " " << m << endl;
        for(int i=0;i<n;++i){
            for(int j=0;j<m;++j){
                cout << a[i][j].val() << " ";
            }
            cout << endl;
        }
    } else {
        vector<vector<mint>> c(n,vector<mint>(rank));
        int pivot=0;
        for(int i=0;i<m;++i){
            if(pivot==n){
                break;
            }
            if(a_copy[pivot][i]==0){
                continue;
            }
            for(int j=0;j<n;++j){
                c[j][pivot]=a[j][i];
            }
            ++pivot;
        }
        vector<vector<mint>> f(rank,vector<mint>(m));
        for(int i=0;i<rank;++i){
            for(int j=0;j<m;++j){
                f[i][j]=a_copy[i][j];
            }
        }
        cout << 2 << endl;
        cout << n << " " << rank << endl;
        for(int i=0;i<n;++i){
            for(int j=0;j<rank;++j){
                cout << c[i][j].val() << " ";
            }
            cout << endl;
        }
        cout << rank << " " << m << endl;
        for(int i=0;i<rank;++i){
            for(int j=0;j<m;++j){
                cout << f[i][j].val() << " ";
            }
            cout << endl;
        }
    }
    return 0;
}
0