結果
問題 | No.2573 moving up |
ユーザー |
![]() |
提出日時 | 2023-11-30 16:13:29 |
言語 | C++23 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 6 ms / 2,000 ms |
コード長 | 2,863 bytes |
コンパイル時間 | 5,801 ms |
コンパイル使用メモリ | 328,264 KB |
実行使用メモリ | 6,548 KB |
最終ジャッジ日時 | 2023-12-06 15:48:25 |
合計ジャッジ時間 | 6,986 ms |
ジャッジサーバーID (参考情報) |
judge13 / judge12 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
6,548 KB |
testcase_01 | AC | 4 ms
6,548 KB |
testcase_02 | AC | 6 ms
6,548 KB |
testcase_03 | AC | 3 ms
6,548 KB |
testcase_04 | AC | 4 ms
6,548 KB |
testcase_05 | AC | 3 ms
6,548 KB |
testcase_06 | AC | 4 ms
6,548 KB |
testcase_07 | AC | 4 ms
6,548 KB |
testcase_08 | AC | 3 ms
6,548 KB |
testcase_09 | AC | 3 ms
6,548 KB |
testcase_10 | AC | 3 ms
6,548 KB |
testcase_11 | AC | 3 ms
6,548 KB |
testcase_12 | AC | 3 ms
6,548 KB |
testcase_13 | AC | 4 ms
6,548 KB |
testcase_14 | AC | 2 ms
6,548 KB |
testcase_15 | AC | 2 ms
6,548 KB |
testcase_16 | AC | 2 ms
6,548 KB |
testcase_17 | AC | 3 ms
6,548 KB |
testcase_18 | AC | 2 ms
6,548 KB |
testcase_19 | AC | 4 ms
6,548 KB |
testcase_20 | AC | 3 ms
6,548 KB |
testcase_21 | AC | 2 ms
6,548 KB |
testcase_22 | AC | 2 ms
6,548 KB |
testcase_23 | AC | 2 ms
6,548 KB |
testcase_24 | AC | 3 ms
6,548 KB |
testcase_25 | AC | 4 ms
6,548 KB |
testcase_26 | AC | 2 ms
6,548 KB |
testcase_27 | AC | 2 ms
6,548 KB |
testcase_28 | AC | 2 ms
6,548 KB |
testcase_29 | AC | 2 ms
6,548 KB |
testcase_30 | AC | 2 ms
6,548 KB |
ソースコード
#include <bits/stdc++.h> using namespace std; #include <atcoder/all> using namespace atcoder; template<typename T> inline bool chmax(T &a, T b) { return ((a < b) ? (a = b, true) : (false)); } template<typename T> inline bool chmin(T &a, T b) { return ((a > b) ? (a = b, true) : (false)); } #define rep(i, n) for (long long i = 0; i < (long long)(n); i++) #define rep2(i, m ,n) for (int i = (m); i < (long long)(n); i++) #define REP(i, n) for (long long i = 1; i < (long long)(n); i++) typedef long long ll; #define updiv(N,X) (N + X - 1) / X #define l(n) n.begin(),n.end() #define YesNo(Q) Q==1?cout<<"Yes":cout<<"No" using P = pair<int, int>; using mint = modint; const int MOD = 998244353LL; const ll INF = 999999999999LL; vector<long long> fact, fact_inv, inv; /* init_nCk :二項係数のための前処理 計算量:O(n) */ template <typename T> void input(vector<T> &v){ rep(i,v.size()){cin>>v[i];} return; } void init_nCk(int SIZE) { fact.resize(SIZE + 5); fact_inv.resize(SIZE + 5); inv.resize(SIZE + 5); fact[0] = fact[1] = 1; fact_inv[0] = fact_inv[1] = 1; inv[1] = 1; for (int i = 2; i < SIZE + 5; i++) { fact[i] = fact[i - 1] * i % MOD; inv[i] = MOD - inv[MOD % i] * (MOD / i) % MOD; fact_inv[i] = fact_inv[i - 1] * inv[i] % MOD; } } /* nCk :MODでの二項係数を求める(前処理 int_nCk が必要) 計算量:O(1) */ long long nCk(int n, int k) { assert(!(n < k)); assert(!(n < 0 || k < 0)); return fact[n] * (fact_inv[k] * fact_inv[n - k] % MOD) % MOD; } long long modpow(long long a, long long n, long long mod) { long long res = 1; while (n > 0) { if (n & 1) res = res * a % mod; a = a * a % mod; n >>= 1; } return res; } ll POW(ll a,ll n){ long long res = 1; while (n > 0) { if (n & 1) res = res * a; a = a * a; n >>= 1; } return res; } int main() { int h,w;cin>>h>>w; vector<ll> x(w); vector<ll> y(w); rep(i,w){cin>>x[i]>>y[i];} ll t = 2;while(t<=w){t*=2;} mcf_graph<ll,ll> g(t*2-1+w+2); // t*2+w+1 // t*2-1 is segtree,w is first,2 is end //cout<<"横幅:" << t << endl; rep(i,w){ ll l = max(0LL,y[i]-x[i])+t;ll r = min(t,y[i])+t; // cout<<l<<" "<<r<<endl; while(l<r){ if(l%2==1){ // cout << "edge:" << l << endl; g.add_edge(t*2+i,l,1LL,0LL); l++;} l /= 2; if(r%2==1){ // cout << "edge:" << r-1 << endl; g.add_edge(t*2+i,r-1,1LL,0LL); r--; } r /= 2; // cout <<"range:" <<l << " " << r << endl; } //g.add_edge(i,w+j,1LL,abs(j-clamp(j,(y[i]-1)-(x[i]-1),y[i]-1))); } rep(i,t-1){g.add_edge(t+i,t+i+1,5000LL,1LL);g.add_edge(t+1+i,t+i,5000LL,1LL);} rep(i,w){g.add_edge(0,t*2+i,1LL,x[i]-1); g.add_edge(t+i,t*2+w,1LL,0LL); } for(int i = 2;i<t*2;i++){ g.add_edge(i/2,i,5000LL,0LL); } cout<<g.flow(0,t*2+w).second<<endl; }