結果
| 問題 |
No.2523 Trick Flower
|
| コンテスト | |
| ユーザー |
Kude
|
| 提出日時 | 2023-10-28 00:04:29 |
| 言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 151 ms / 4,000 ms |
| コード長 | 2,087 bytes |
| コンパイル時間 | 3,381 ms |
| コンパイル使用メモリ | 281,752 KB |
| 実行使用メモリ | 34,864 KB |
| 最終ジャッジ日時 | 2024-09-25 15:40:54 |
| 合計ジャッジ時間 | 6,231 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 31 |
ソースコード
#include<bits/stdc++.h>
namespace {
#pragma GCC diagnostic ignored "-Wunused-function"
#include<atcoder/all>
#pragma GCC diagnostic warning "-Wunused-function"
using namespace std;
using namespace atcoder;
#define rep(i,n) for(int i = 0; i < (int)(n); i++)
#define rrep(i,n) for(int i = (int)(n) - 1; i >= 0; i--)
#define all(x) begin(x), end(x)
#define rall(x) rbegin(x), rend(x)
template<class T> bool chmax(T& a, const T& b) { if (a < b) { a = b; return true; } else return false; }
template<class T> bool chmin(T& a, const T& b) { if (b < a) { a = b; return true; } else return false; }
using ll = long long;
using P = pair<int,int>;
using VI = vector<int>;
using VVI = vector<VI>;
using VL = vector<ll>;
using VVL = vector<VL>;
std::vector<std::vector<int>> scc_dag(std::vector<std::vector<int>>& scc, std::vector<std::vector<int>>& graph) {
int n = graph.size();
int m = scc.size();
std::vector<int> d(n);
for(int i = 0; i < m; ++i) for(int v: scc[i]) d[v] = i;
std::vector<std::vector<int>> to(m);
std::vector<int> added(m, -1);
for(int i = 0; i < m; ++i) {
added[i] = i;
for(int u: scc[i]) for(int x: graph[u]) {
int v = d[x];
if (added[v] != i) {
added[v] = i;
to[i].push_back(v);
}
}
}
return to;
}
} int main() {
ios::sync_with_stdio(false);
cin.tie(0);
int n;
cin >> n;
VI a(n), b(n), c(n);
rep(i, n) cin >> a[i];
rep(i, n) cin >> b[i];
scc_graph g(n);
VVI to(n);
rep(i, n) {
int c;
cin >> c;
g.add_edge(c - 1, i);
to[c-1].emplace_back(i);
}
auto scc = g.scc();
int sz = scc.size();
auto toc = scc_dag(scc, to);
ll l = 0, r = accumulate(all(a), 0LL) / accumulate(all(b), 0LL) + 1;
while(r - l > 1) {
ll mid = (l + r) / 2;
VL dp(sz);
rrep(i, sz) {
for(int j : toc[i]) if (dp[j] < 0) dp[i] += dp[j], dp[j] = 0;
for(int x : scc[i]) dp[i] += a[x] - mid * b[x];
}
bool ok = all_of(all(dp), [](ll x) { return x >= 0; });
(ok ? l : r) = mid;
}
cout << l << '\n';
}
Kude