結果
問題 | 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 4000000000000000001 int 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; }