結果
問題 | No.572 妖精の演奏 |
ユーザー |
![]() |
提出日時 | 2020-05-02 20:08:29 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 6 ms / 2,000 ms |
コード長 | 2,417 bytes |
コンパイル時間 | 3,202 ms |
コンパイル使用メモリ | 203,580 KB |
最終ジャッジ日時 | 2025-01-10 06:06:36 |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 20 |
ソースコード
//#include <tourist>#include <bits/stdc++.h>using namespace std;typedef long long ll;typedef pair<ll, ll> p;const int INF = 1e9;const ll LINF = ll(1e18);const int MOD = 1000000007;const int dx[4] = {0, 1, 0, -1}, dy[4] = {-1, 0, 1, 0};const int Dx[8] = {0, 1, 1, 1, 0, -1, -1, -1}, Dy[8] = {-1, -1, 0, 1, 1, 1, 0, -1};#define yes cout << "Yes" << endl#define YES cout << "YES" << endl#define no cout << "No" << endl#define NO cout << "NO" << endl#define rep(i, n) for (int i = 0; i < n; i++)#define ALL(v) v.begin(), v.end()#define debug(v) \cout << #v << ":"; \for (auto x : v) \{ \cout << x << ' '; \} \cout << endl;template <class T>bool chmax(T &a, const T &b){if (a < b){a = b;return 1;}return 0;}template <class T>bool chmin(T &a, const T &b){if (b < a){a = b;return 1;}return 0;}//cout<<fixed<<setprecision(15);有効数字15桁//-std=c++14//-std=gnu++17ll gcd(ll a, ll b) { return b ? gcd(b, a % b) : a; }ll lcm(ll a, ll b) { return a / gcd(a, b) * b; }ll n, m;vector<vector<ll>> a;int main(){cin.tie(0);ios::sync_with_stdio(false);cin >> n >> m;a.resize(m);rep(i, m){a[i].resize(m);rep(j, m){cin >> a[i][j];}}vector<vector<vector<ll>>> dp(40, vector<vector<ll>>(m, vector<ll>(m,0)));rep(i, m){rep(j, m) dp[0][j][i] = a[i][j];}vector<vector<ll>>ans(m,vector<ll>(m,0));for (int i = 1; i < 40; i++){rep(j, m){rep(k, m){rep(z, m){chmax(dp[i][j][k], dp[i - 1][j][z]+dp[i-1][z][k]);}}}}n--;//1回目は遷移ではないですint d=0;while(n){if(n%2){vector<vector<ll>> ans2(m, vector<ll>(m, 0));rep(i,m){//jからiの遷移rep(j,m){rep(k,m){chmax(ans2[i][j],ans[i][k]+dp[d][k][j]);}}}swap(ans,ans2);}n/=2;d++;}ll ansa=0;rep(i,m){rep(j,m){chmax(ansa,ans[i][j]);}}cout<<ansa<<"\n";}