結果
問題 |
No.3230 Mutual Corresponding System
|
ユーザー |
|
提出日時 | 2025-08-08 22:28:13 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 150 ms / 2,000 ms |
コード長 | 2,419 bytes |
コンパイル時間 | 2,354 ms |
コンパイル使用メモリ | 199,708 KB |
実行使用メモリ | 7,716 KB |
最終ジャッジ日時 | 2025-08-08 22:28:17 |
合計ジャッジ時間 | 4,181 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 23 |
ソースコード
#include<bits/stdc++.h> using namespace std; //#pragma GCC optimize("Ofast") #define rep(i,n) for(ll i=0;i<n;i++) #define repl(i,l,r) for(ll i=(l);i<(r);i++) #define per(i,n) for(ll i=(n)-1;i>=0;i--) #define perl(i,r,l) for(ll i=r-1;i>=l;i--) #define fi first #define se second #define pb push_back #define ins insert #define pqueue(x) priority_queue<x,vector<x>,greater<x>> #define all(x) (x).begin(),(x).end() #define CST(x) cout<<fixed<<setprecision(x) #define vtpl(x,y,z) vector<tuple<x,y,z>> #define rev(x) reverse(x); using ll=long long; using vl=vector<ll>; using vvl=vector<vector<ll>>; using pl=pair<ll,ll>; using vpl=vector<pl>; using vvpl=vector<vpl>; const ll MOD=1000000007; const ll MOD9=998244353; const int inf=1e9+10; const ll INF=4e18; const ll dy[9]={1,0,-1,0,1,1,-1,-1,0}; const ll dx[9]={0,1,0,-1,1,-1,1,-1,0}; template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return true; } return false; } template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return true; } return false; } //O(N^3logk) //https://yukicoder.me/submissions/733933 //https://atcoder.jp/contests/abc236/submissions/28778130 using vec=vector<ll>; using matrix=vector<vector<ll>>; matrix matmul(matrix &left,matrix &right){ int n=left.size(); matrix ret(n,vec(n,-INF));//+の単位元で埋める for(int i=0;i<n;i++){ for(int k=0;k<n;k++){//ループ順序の最適化 for(int j=0;j<n;j++){ chmax(ret[i][j],left[i][k]+right[k][j]); } } } return ret; } matrix matpow(matrix &mat,ll k){ int n=mat.size(); matrix ret(n,vec(n,-INF));//+の単位元で埋める for(int i=0;i<n;i++)ret[i][i]=0;//*の単位元で埋める while(k){ if(k&1)ret=matmul(mat,ret); k>>=1; if(k)mat=matmul(mat,mat); } return ret; } vec update(matrix &mat,vec &v){ int n=mat.size(); vec ret(n);//+の単位元で埋める for(int i=0;i<n;i++){ for(int j=0;j<n;j++){ chmax(ret[i],mat[i][j]+v[j]); } } return ret; } int main(){ ll n,m;cin >> n >> m; vl s(n); rep(i,n)cin >> s[i]; vvl mat(n,vl(n)); rep(i,n)rep(j,n)cin >> mat[j][i]; mat=matpow(mat,m); /*rep(i,n){ rep(j,n)cout << mat[i][j] <<" ";cout << endl; }*/ s=update(mat,s); rep(i,n)cout << s[i] <<" ";cout << endl; }