結果

問題 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
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#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;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0