結果
問題 | No.2523 Trick Flower |
ユーザー |
![]() |
提出日時 | 2023-10-27 23:20:32 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 320 ms / 4,000 ms |
コード長 | 1,222 bytes |
コンパイル時間 | 4,287 ms |
コンパイル使用メモリ | 259,004 KB |
最終ジャッジ日時 | 2025-02-17 15:59:59 |
ジャッジサーバーID (参考情報) |
judge2 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 31 |
ソースコード
#include <stdio.h>#include <atcoder/all>#include <bits/stdc++.h>using namespace std;using namespace atcoder;using mint = modint998244353;#define rep(i,n) for (int i = 0; i < (n); ++i)#define Inf32 1000000001#define Inf64 4000000000000000001int main(){int n;cin>>n;vector<long long> a(n),b(n);rep(i,n)cin>>a[i];rep(i,n)cin>>b[i];vector<int> c(n);scc_graph S(n);rep(i,n){cin>>c[i];c[i]--;S.add_edge(c[i],i);}auto s = S.scc();vector<int> pos(n);vector<long long> aa(s.size()),bb(s.size());vector<vector<int>> E(s.size());rep(i,s.size()){rep(j,s[i].size()){pos[s[i][j]] = i;aa[i] += a[s[i][j]];bb[i] += b[s[i][j]];}}rep(i,n){if(pos[i]==pos[c[i]])continue;E[pos[i]].push_back(pos[c[i]]);}long long ok = 0,ng = 2000000000000000;while(ng-ok>1LL){long long mid =(ok+ng)/2;bool f = true;auto ta = aa,tb = bb;vector<__int128> need(s.size());for(int i=s.size()-1;i>=0;i--){__int128 t = bb[i];t *= mid;need[i] += t;if(ta[i]>=need[i]){continue;}if(E[i].size()==0){f = false;break;}need[E[i][0]] += need[i] - ta[i];}if(f)ok = mid;else ng = mid;}cout<<ok<<endl;return 0;}