#define _USE_MATH_DEFINES #include using namespace std; #define FOR(i,m,n) for(int i=(m);i<(n);++i) #define REP(i,n) FOR(i,0,n) #define ALL(v) (v).begin(),(v).end() using ll = long long; constexpr int INF = 0x3f3f3f3f; constexpr long long LINF = 0x3f3f3f3f3f3f3f3fLL; constexpr double EPS = 1e-8; constexpr int MOD = 1000000007; // constexpr int MOD = 998244353; constexpr int dy[] = {1, 0, -1, 0}, dx[] = {0, -1, 0, 1}; constexpr int dy8[] = {1, 1, 0, -1, -1, -1, 0, 1}, dx8[] = {0, -1, -1, -1, 0, 1, 1, 1}; template inline bool chmax(T &a, U b) { return a < b ? (a = b, true) : false; } template inline bool chmin(T &a, U b) { return a > b ? (a = b, true) : false; } struct IOSetup { IOSetup() { std::cin.tie(nullptr); std::ios_base::sync_with_stdio(false); std::cout << fixed << setprecision(20); } } iosetup; template struct Dinic { struct Edge { int dst, rev; T cap; Edge(int dst, T cap, int rev) : dst(dst), cap(cap), rev(rev) {} }; std::vector> graph; Dinic(int n) : graph(n), level(n), itr(n) {} void add_edge(int src, int dst, T cap) { graph[src].emplace_back(dst, cap, graph[dst].size()); graph[dst].emplace_back(src, 0, graph[src].size() - 1); } T maximum_flow(int s, int t, T limit) { T res = 0; while (true) { std::fill(level.begin(), level.end(), -1); std::queue que; level[s] = 0; que.emplace(s); while (!que.empty()) { int ver = que.front(); que.pop(); for (const Edge &e : graph[ver]) { if (level[e.dst] == -1 && e.cap > 0) { level[e.dst] = level[ver] + 1; que.emplace(e.dst); } } } if (level[t] == -1) return res; std::fill(itr.begin(), itr.end(), 0); T f; while ((f = dfs(s, t, limit)) > 0) res += f; } } private: std::vector level, itr; T dfs(int ver, int t, T flow) { if (ver == t) return flow; for (; itr[ver] < graph[ver].size(); ++itr[ver]) { Edge &e = graph[ver][itr[ver]]; if (level[ver] < level[e.dst] && e.cap > 0) { T tmp = dfs(e.dst, t, std::min(flow, e.cap)); if (tmp > 0) { e.cap -= tmp; graph[e.dst][e.rev].cap += tmp; return tmp; } } } return 0; } }; template